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C0MPF2~A  PROGRAM  FOR  CALCULATING 
POST-FLASHOVER  FIRE  TEMPERATURES 

Vytenis  Babrauskas 


C0MPF2  is  a  computer  program  for  calculating  the 
characteristics  of  a  post-f lashover  fire  in  a  single 
building  compartment,  based  on  fire-induced  ventilation 
through  a  single  door  or  window.   It  is  intended  both 
for  performing  design  calculations  and  for  the  analysis 
of  experimental  burn  data.   Wood,  thermoplastic,  and 
liquid  fuels  can  be  treated.   In  addition  to  the  capa- 
bility of  performing  calculations  for  compartments  with 
completely  determined  properties,  routines  are  included 
for  calculating  fire  behavior  by  an  innovative  variable 
abstraction  method.   A  comprehensive  output  format  is 
provided  which  gives  gas  temperatures,  heat  flow  terms, 
and  flow  variables.   The  documentation  includes  input 
instructions,  sample  problems,  and  a  listing  of  the 
program.   The  program  is  written  in  Fortran  and  consti- 
tutes an  improved  version  of  an  earlier  program,  COMPF . 

Key  words:  Computer  programs--f ire  protection;  fire 
protection;  fire  resistance;  fire  tests;  fire  walls; 
safety  engineering--f ires. 


1.   INTRODUCTION 

With  increasing  efforts  [1-4]   towards  rational  methods  of  providing 
fire  endurance  for  structural  building  components,  it  becomes  highly  desir- 
able for  both  the  designer  and  the  researcher  to  have  available  computer 
programs  for  calculating  expected  fire  temperatures  and  heat  transfer  through 
the  building  components.   A  fire  is  not  considered  as  becoming  a  threat  to  a 
structure  and  its  fire  barriers  until  it  reaches  the  flashover  stage.   Flash- 
over  of  a  room  is  defined  as  that  fire  stage  when  the  bulk  of  the  room  volume 
becomes  involved  in  flames.   Operationally,  this  roughly  coincides  with  flames 
coming  out  the  door  or  window,  or  an  upper  gas  space  temperature  of  around 
600°  C,  or  a  radiant  heat  flux  at  floor  level  of  about  20  kW/m^ .   For  the 
purpose  of  designing  for  fire  endurance,  then,  only  post-f lashover  fires  are 
considered.   The  present  report  describes  a  computer  program  for  calculating 
the  expected  temperatures,  heat  and  mass  flows  and  other  variables  in 
post-f lashover  building  fires.   Different  routines  are  incorporated  for 
producing  design  time-temperature  curves  and  for  permitting  comparative 
theoretical  curves  to  be  generated  based  on  experimental  mass  loss  rates. 

2.   HISTORY  OF  DEVELOPMENT 

The  first  computer  program  for  calculating  post-f lashover  fire  tempera- 
tures was  developed  by  Kawagoe  [5] ,  in  conjunction  with  his  pioneering  studies 
leading  to  a  theoretical  room  fire  model.   This  model  was  an  adaptation  of  an 
earlier  graphical  technique.   The  main  limitations  of  both  the  computer  program 


Numbers  in  brackets  refer  to  the  literature  references  listed  at  the  end  of 
this  paper. 
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and  the  theoretical  model  was  the  restriction  to  ventilation-limited  fires. 
Fuel-limited  fires  could  only  be  expressed  in  terms  of  an  empirical  tempera- 
ture change  rate.   Magnusson  and  Thelandersson  [6]  studied  heat  release  rates 
in  more  detail  and  produced  a  model.   An  unpublished  computer  program  was 
used  to  implement  that  model.   The  normalized  shape  of  the  fire  time-temperature 
curves  was  an  input  variable  in  this  model;  the  shape  was  based  on  sets  of 
typical  measurements.   Based  on  Magnusson  and  Thelandersson' s  theory,  Fedock 

[7]  published  a  similar  computer  program  with  emphasis  on  prestressed  concrete 
structures.  The  first  program  to  provide  for  theoretically  based  calculations 
of  both  ventilation-limited  and  fuel-limited  burning  was  written  by  Tsuchiya 

[8] .   It  was  restricted  to  fires  starting  in  ventilation  control  and  to  fuel 
consisting  of  sparsely-packed  wood  sticks. 

The  predecessor  to  the  present  program,  COMPF  [9] ,  was  issued  in  1975 
and  incorporated  several  new  advances,  including  the  ability  to  treat  entirely 
fuel-limited  fires,  to  allow  for  temperature-dependent  wall  properties,  to 
permit  the  optional  use  of  numerical  input  fuel  weight  loss  rates,  and  to 
perform  certain  variable  abstraction  ("pessimization")  calculations  as  an  aid 
to  design.   (These  techniques  enable  an  input  variable  to  be  treated  non- 
deterministically. )   Program  C0MPF2  is  intended  to  replace  program  COMPF  and 
differs  from  it  in  the  following  main  ways: 

1.  A  subroutine  has  been  added  to  allow  treatment  of  fires  where 
thermoplastic  or  liquid  fuel  exists  in  the  form  of  a  pool  on  the 
floor.   The  routine  implements  the  theory  discussed  in  reference  [10] 
and  outlined  in  section  4.1;  examples  of  calculations  are  also  given 
and  discussed  in  that  reference. 

2.  The  deterministic  wood  fuel  burning  model  has  been  extended  to  include 
the  possibility  of  densely-packed  cribs. 

3.  Both  pool  fire  and  densely-packed  crib  options  have  been  incorporated 
into  the  pessimization  routines. 

'4.   In  addition  to  performing  transient  calculations,  the  program  can 

now  also  treat  steady-state  solutions,  for  both  lossy  and  adiabatic 
walls. 

5.  The  program  is  now  in  S.I.  units  throughout. 

6.  Certain  corrections  and  improvements  have  been  incorporated  in  the 
calculation  routines.   The  method  for  the  iterative  solution  of  the 
heat  balance  equation  has  especially  been  improved. 

3 .   THEORY 

The  post-f lashover  compartment  fire  theory  has  been  given  in  some  detail 
in  reference  [11],  thus,  only  a  brief  summary  will  be  given  here.   The  main 
assumptions  are: 

•  The  compartment  represents  a  well-stirred  reactor,  i.e.,  spatial 
temperature  variations  in  the  hot  fire  gases  are  ignored. 

•  The  model  is  quasi-steady.   Time  variations  in  fuel  release  rate  and 
in  conduction  losses  are  fully  included.   However,  time  rate  of  change 
terms  in  gas  phase  mass  and  energy  balance  are  dropped. 

•  Air  supply  and  gas  outflow  is  through  a  single  window  in  a  vertical  wall 
and  is  the  result  of  fire-induced  convection. 

The  thermal  discontinuity  away  from  the  window  region  is  at  a  level  below 
the  bottom  of  the  window.   The  volume  below  the  discontinuity  is  occupied 


by  cold  incoming  air.   In  a  flashed-over  fire  this  discontinuity  is  close 
to  the  floor.   Its  exact  location  below  the  window  bottom  is  immaterial  [12] 

Burning  is  limited  by  rates  of  air  or  fuel  supply  rather  than  by  gas 
phase  chemical  kinetics. 

Walls  (including  the  ceiling)  are  modeled  as  portions  of  a  homogeneous 
solid  of  finite  thickness.   Temperature-dependent  material  properties 
are  allowed  for. 

The  heat  balance  equation  is : 


^c  -  "^f  ^^T. 
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where  h  denotes  enthalpy  and  the  definition  of  symbols  is  given  in  the 
Nomenclature  section.   The  subscripts  on  the  enthalpy  terms  denote  the  tempera- 
ture at  which  they  are  evaluated.   The  window  radiation  loss  is,  simply 
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The  wall  loss  term  has  a  radiative  and  a  convective  component, 

-1  k  k 

Q   =  A   [o  -^, =--i-7 =-   (T^  -  T  )  +  h  (T^  -  T  )]        (3) 

w     w      1/e^  +  1/e  -If     w         f     w         ^ 
r       w 

The  convective  coefficient  h  above  is  not  well  known  since  the  exact  flow 
conditions  at  the  wall  and  ceiling  surfaces  in  a  post-f lashover  fire  are  not 
known  in  detail.   The  convective  fraction  is  much  less  than  the  radiative 
fraction,  permitting  a  rather  simplified  treatment.   For  turbulent-free  con- 
vection flow  over  flat  plates  the  value  for  h  should  depend  [13]  on 
(T_  -  T  )  ^^3^  p^   value  of 
f     w 
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1     w 
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was  selected  as  being  in  reasonable  agreement  with  data. 

The  analysis  in  [11]  shows  that  for  compartments  greater  than  about  2  m 
on  a  side,  a  flame  emissivity  of  e^=0.9  may  be  used. 


of 


The    enthalpy   evolved    from   combustion,    h    ,    must   be    evaluated   as    the    lesser 


m     A  h      b 
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When  equation  (5)  is  limiting,  the  combustion  is  known  as  "fuel-limited", 
while  if  equation  (6)  is  smaller,  combustion  is  "ventilation-limited".   Here 
bp  represents  the  maximum  combustion  efficiency  and  is  a  largely  unknown 
number.   Since  it  represents,  effectively,  the  "unmixedness"  of  the  combustion, 
there  may  be  a  scale  effect,  with  smaller  compartment  spaces  showing  less 
complete  mixing.   Experimental  data  can  generally  be  correlated  within  a  range 

of  about  0.7  <   b   <  0.9. 
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The  value  of  m  .   is  obtained  from  the  Bernoulli  equation  at  the  window 
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The  discharge  coefficient  has  been  determined  by 
to  be  0.68  for  normal-shaped  windows.   This  value  does 
where  the  window  takes  up  almost  an  entire  wall.   For 
patterns  have  not  been  studied,  but  data  can  be  correl 
about  one-half  its  normal  value.   The  molecular  weight 
is  not  exactly  known  since  the  composition  of  the  gase 
unburned  fuel  gases  (excess  pyrolysates)  is  generally 
their  molecular  weight  has  been  assumed  equal  to  that 
contribution  of  carbon  monoxide  and  other  minor  combus 
ignored.   The  main  dependence  of  m^ij-  is  on  the  window 
reasonable  values  of  temperature  the  whole  expression 
equal  to 


Prahl  and  Emmons  [14] 
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But  this  approximation  has  not  been  employed  here. 

Ventilation  through  multiple  openings  has  not  been  provided  for  in  this 
program.   An  approach  for  treating  such  problems  is  given  in  [11] . 

The  heat  of  combustion.  Ah  ,  is  taken  as  the  net  value  since  the  hot  gas 
outflow  is  above  100°  c-   The  stoichiometric  ratio,  r,  is  a  constant  for  a  pure 
material;  a  tabulation  of  values  is  given  in  [15]. 

The  rate  of  pyrolysis,  Ap ,  is  one  of  , the  hardest  quantities  to  determine. 
A  discussion  of  available  values  is  given  in  the  next  section. 
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The  outflow  mass  rate,  mf ,  is  by  mass  conservation  the  sum  of  A^ir 
The  enthalpy  of  the  outflow  products,  hm   and  h     " 
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and 
is  evaluated  on  the 


assumption  that  the  combusted  fuel  goes  to  CO2  and  H„0.   No  account  is 
taken  of  CO  for  two  reasons:   because  the  effect  on  a  mass  basis  would  be 
very  small,  and  because  it  was  considered  advisable  not  to  introduce  any 
reaction  kinetics.   Also,  only  elements  C,  H,  0,  and  N  have  been  considered 
for  the  fuel  composition. 


The  excess  pyrolysate  term  Qgp/  is  the  heat  required  to  vaporize  the 

excess  pyrolysates.   Note  that  with  the  conventional  definition  of  heat  of 

combustion,  the  loss  for  vaporization  of  combusted  pyrolysates  is  already 

included  in  A  h  . 
c 

The  second  major  equation  to  be  solved  is  for  heat  conduction  through 
the  wall. 


3T     ^     3T        ,  ,  , 


The  wall  is  initially  at  ambient  temperature,  T  ,  and  is  subjected  to 
boundary  conditions  at  the  fire  side  of:  ° 


w 
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_k  __  =  h  [T^  -  T^(o)]  +  ea     [T^  -  T^(o)]  (10) 

and  on  the  unexposed  side  (x  =  L) , 

ci  T 

-^  ^  =  ^  f^w(^)  -  ^o^  ^  ^^  f^w  (L)  -  T^  (11) 

For  the  fire  side  the  convective  coefficient  has  been  given  above.   For 
the  unexposed  side  a  value  of 

h  =  1.87   [T   -  T   (L) ] 1^3  (12) 

o     w  ^   ' 

was  taken. 

4.   PYROLYSIS  RATES 

4.1   Liquid  or  Thermoplastic  Pools 

There  is  currently  only  one  fuel  arrangement  where  the  pyrolysis  rate  may 
adequately  be  predicted  from  theory.   It  consists  of  a  pool  of  thermoplastic 
or  liquid  fuel  on  the  floor.   The  fuel  is  pyrolyzed  solely  by  radiant  flux 
and  "sees"  the  compartment  with  a  view  factor  of  1.0  and  itself  with  a  view 
factor  of  zero.   In  addition,  the  fuel  must  pyrolyze  at  a  known  surface 
temperature,  T,  ,  and  with  a  known  heat  of  pyrolysis.  Ah  .   Then: 

m   =  A.  ^^  (^f  -  ^b)  (13) 
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Tewarson  and  Pion  [16]  have  measured  heats  of  pyrolysis  for  numerous 
thermoplastic  materials. 

The  above  simple  model  is  fully  adequate  for  steady-state  solutions. 
At  the  start  of  the  fire,  however,  the  radiation  feedback  is  small  from  the 
hot  gas  volume  but  may  be  larger  from  the  local  plume  above  the  pool  itself. 
Thus,  a  plume  term  should  be  added  in  to  model  the  starting  transient.   Very 
limited  experimental  data  by  Burgess  [17]  and  J^y  Modak  [18]  can  be  used  to 
derive  an  empirical  relationship  for  the  plume  pyrolysis  rate  as: 


m  =  A   0.0014  ^^c   (kg/s)  (14) 

Ah 
P 

This  relationship  does  not  take  into  account  differences  in  flame 
emissivities  for  various  materials;  as  a  result,  it  only  provides  a  crude 
measure.   In  the  present  application,  however,  the  contribution  of  this  term 
is  minor;  therefore,  an  approximate  expression  is  adequate.   Also,  as  the 
room  radiation  increases,  the  effect  of  plume  radiation  on  pyrolyzing  the 
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fuel  decreases.   For  a  radiatively  black  room,  at  high  temperature,  the  plume 
term  should  properly  be  negligible.   This  interaction  is  crudely  modeled  by 
multiplying  the  plume  term  by  a  proportionality  factor  before  adding  to  the 
far-field  term.   The  proportionality  factor,  X/  has  been  set  equal  to 


X  =  1.0 


4  4 

T   -  T 

f  b  (15) 

IVOO'*  -  Tj^ 


With  x>0. 

4. '2   Solid  Fuels 

Empirical  data  are  available  for  the  mass  loss  rates  of  wood  planks  in 
flashed-over  fires.   Because  of  the  nature  of  wood  combustion,  these  rates  are 
not  especially  sensitive  to  room  radiation  and  can  be  specified  [11]  using 
a  regression  velocity  of  7-15  ym/s.   This  relationship  is  adequate  to  describe 
the  burning  of  large,  isolated  wood  panels.   For  pieces  thin  in  two  or  three 
dimensions,  yet  still  widely  spaced,  the  following  expression  is  suitable: 


1-1/^  (16) 

Here  M   is  the  original  mass,  m  is  the  mass  at  a  given  time,  and  F  is  a 
constant  equal  to  2  for  cylinders  or  rectangular  sticks  and  equal  to  3  for 


m 

P   . 

.    F     .    m 
'    C     ^    M 

M 

o 

spheres  or  cubes.   C  is  given  by 


c  =  ° 


^^p 


with  D  being  the  smallest  fuel  dimension  and  v   the  regression  velocity.   For 
thin  fuels  v   is  approximately 

_6   _0 . 6 
v   -  1.7  X  10    D     (m/s)  (17) 

The  final  arrangement  for  wood  fuel  for  which  data  are  available  is  a 
crib,  or  a  regular  stacked  array.   From  the  data  of  Nilsson  [19]  and  Yamashika 
[20] ,  a  set  of  simplified  relationships  has  been  evolved  for  the  three  crib 
burning  regimes. 

Fuel  Surface  Control 

%  =  ^  ^p  (  S   )  ^  ^o  (^^5/^)  (^^) 

^  ^    o 

_6   _0. 6 
V   =  1.7  X  10    D 
P 


Crib 

Porosity  Control 

*p  = 

4.4  X  10-'  (Vh  )  %                                                                                                     (,9) 
c   D 

\ 

=  ratio  of  stick  clear  spacing  to  crib  height 

Room 

Ventilation  Control 

m   = 
P 

0.12  A^  /h^                                                 (20) 

In  calculations,  each  of  the  three  rates  above  are  determined  and  the 

lowest  rate  taken  as  governing. 

5.   DETAILS  OF  SUBROUTINES 

The  program  routines  are  written  in  Fortran  language.    A  complete 

listing  is  given  in  appendix  B.   The  following  are  brief  descriptions  of  the 

operation 

of  each  subroutine. 

5.1   C0MPF2 

C0MPF2  is  the  main  program.   It  calls  most  of  the  calculational  routines. 

A  flow  chart  of  C0MPF2  is  given  in  figure  1.   The  program  starts  with  the 

initialization  of  certain  constants  and  default  values.   The  input  title  and 
namelist  are  then  read  in.   If  tabular  data  are  specified,  subroutine  INC  is 
called.   ICONDS  is  then  called  in  to  set  initial  starting  values.   The  input 
data  are  echoed  in  ECHOID.   After  that,  the  appropriate  computational  routine 
is  called  in.   If  no  iteration  failure  has  occurred  the  program  then  loops 
back  to  the  start  and  goes  to  the  next  problem.   In  case  of  iteration  failure, 
the  program  returns  to  the  same  problem,  this  time  printing  out  additional 
intermediate  calculation  values.   This  intermediate  output  can  also  be  forced 
to  appear  by  specifying  KTRACE=1. 

5.2  CRIB 

Subroutine  CRIB  calculates  the  burning  of  wood  crib  fires.   A  trial  gas 
temperature  value  is  assumed  for  the  first  time  step.   This  value  is  preset, 
but  may  be  overridden  by  specifying  a  value  of  TINPT.   The  flow  quantities 
are  computed,  then  the  wall  losses  are  determined  by  calling  DESOLV.   The  heat 
balance  is  then  determined.   If  the  normalized  residue  is  greater  than  0.002, 
the  iteration  continues.   The  new  temperature  is  normally  determined  by  the 
Newton  method.   If  divergence  results,  a  scanning  technique  is  used  initially 
and  a  splitting  of  differences  once  a  bounded  oscillation  results.   After 
successful  convergence  a  new  wall  temperature  profile  is  established  by  calling 
RSTA.   The  calculation  then  proceeds  to  the  next  time  step.   Computation  is 
terminated  at  the  end  of  time  MTIME,  or  when  gas  temperature  drops  to  353  K, 
or  if  errors  or  convergence  failure  is  detected. 

5.3  DEQNS 

Subroutine  DEQNS  computes  wall  heat  conduction  using  the  Crank-Nicolson 
method  [21].   DEQNS  has  two  entry  points:   DESOLV  and  RSTA.   The  radiation 
boundary  condition  is  linearized;  updating  every  iteration  rather  than  every 
time  step  ensures  minimal  error.   An  additional  within-loop  iteration  is  also 
used.   DEQNS  calls  TRIDGF  to  solve  the  equation  matrix. 


5.4  ECHOID 

Subroutine  ECHOID  echoes  the  input  data.   The  complete  data  set  is  given 
for  each  run,  rather  than  just  the  changed  values.   Care  has  been  taken  to 
give  physical  meaning  for  the  variables  printed. 

5.5  ICONDS 

Subroutine  ICONDS  initializes  starting  values  and  does  some  preliminary 
calculations  on  the  input  data.   It  also  makes  a  few  checks  on  the  validity 
of  the  input  data.   The  user,  however,  is  cautioned  that  this  checking  is 
very  rudimentary  and  in  case  of  error  exit  or  iteration  failure  the  input 
data  must  be  carefully  examined. 

5.6  INC 

Subroutine  INC  is  called  in  when  tabular  input  data  are  to  be  read. 

5.7  OUTPUT 

Subroutine  OUTPUT  is  the  primary  output  routine.   It  writes  at  each 
time  step  a  large  number  of  variables  to  output  files  (logical  units)  2  and 
3.   The  temperatures,  burning  rates,  and  other  primary  variables  are  put  on 
file  2,  while  the  heat  balance  values  and  the  mass  fractions  are  v/ritten  on 
file  3.   OUTPUT  also  converts  temperatures  from  Kelvin  to  degrees  Celsius 
before  printing  them  out. 

5.8  PFLFIX 

PFLFIX  is  a  pessimization  design  routine.   Fuel  pyrolysis  rate  is 
calculated  according  to  governing  equations,  but  the  ventilation  is  pessimized 
by  instantaneously  adjusting  the  window  width  to  give  the  highest  possible 
temperatures.   Wood  stick  or  wood  crib  fuel  is  assumed  unless  PLFUEL=T, 
in  which  case  a  pool  fire  is  used.   The  window  width  is  not  allowed  to  exceed 
a  maximum,  as  set  by  AWDOW/HWDOW.   Calculations  stop  when  the  fuel,  as 
specified  by  FLOAD,  is  exhausted,  since  the  window  width  would  be  undefined 
beyond  that  point.   Calculational  procedures  are  similar  to  those  in  CRIB. 

5 . 9  POOL 

POOL  is  a  pool  fire  burning  routine.   Computational  details  are  similar 
to  those  as  in  CRIB.   The  pyrolysis  rate  is  based  on  equations  13,  14,  and  15. 
Three  modes  of  subroutine  operation  are  possible.   If  STOICH=T,  the  steady- 
state  temperatures  and  pool  area  are  determined  for  stoichiometric  burning. 
If  EISCAN=T,  the  steady-state  solution  is  found  for  a  given  pool  area  greater 
than  stoichiometric.   The  pool  area  is  specified  by  use  of  the  parameter  EITA, 
defined  as  [10] 


A   /h 

V     V 


n  = 


A   /h 

V     V 


(21) 


stoich 


For  constant  window  size,  this  becomes  simply  a  ratio  of  pool  areas.  No 

solutions  are  possible  for  n^l-   Finally,  a  transient  calculation  can  be  made, 

which  proceeds  similarly  as  in  the  other  transient  calculations.   The  user 

must  make  sure  that  the  pool  size  given  is  sufficiently  large  so  that  n^l. 


5.10   PP 

Subroutine  PP  is  a  plotting  routine.   Details  are  not  given  since  plotting 
routines  are  dependent  on  the  hardware  used. 

5.11  PVTFIX 

PVTFIX  is  a  pessimization  routine,  and  is  effectively  the  inverse  of 
PFLFIX.   In  PVTFIX  a  fixed  ventilation  opening  is  specified.   The  fuel  release 
rate  is  instantaneously  varied  to  always  result  in  the  highest  possible  burning 
temperature.   Temperatures  drop  after  the  fuel  load  is  consumed.   Computational 
details  are  similar  to  those  in  PFLFIX. 

5.12  RPFIX 

For  comparison  of  measured  data  against  numerical  predictions  a  routine 
is  needed  which  can  accept  m   rates  as  an  input  tabular  function  of  time. 
RPFIX  provides  for  this  type  of  checking  calculation.   The  case  of  measured 
combustion  rate  input  (as  provided,  for  instance,  by  oxygen  depletion  measure- 
ments in  the  window  outflow)  can  also  be  treated  by  dividing  the  measured  rate 
by  A h   (net)  and  setting  b   =  1.0. 

5.13  STFLOW 

Subroutine  STFLOW  is  a  wall  heat  conduction  routine.   It  is  similar  to 
DEQNS,  except  that  only  the  steady-state  temperatures  are  determined. 

5.14  TLU 

Function  TLU  is  a  tabular  data  interpolating  function  used  in  several 
subroutines.   If  the  independent  variable  entered  is  smaller  than  the  smallest 
data  point  or  larger  than  the  largest  data  point,  the  output  is  set  equal  to 
the  smallest,  or  largest  dependent  value,  respectively. 

5.15  TRIDGF 

Subroutine  TRIDGF  uses  a  Gauss  elimination  procedure  to  solve  a  set  of 
tri-diagonal  matrix  equations. 

6 .   AGREEMENT  WITH  EXPERIMENT 

A  comparison  of  numerical  predictions  with  experimental  results  has  been 
given  in  [22]  for  the  program  COMPF.   Similar  agreement  should  hold  for  C0MPF2, 
since  C0MPF2  is  improved  mainly  in  operational  features,  especially  increased 
versatility,  while  retaining  the  same  theoretical  model  as  in  COMPF.   For  pool 
fires  useable  full-scale  experimental  data  are  not  available. 

7.   INPUT  INSTRUCTIONS 

7.1   Deck  Set-up 

The  input  is  assigned  to  file  1.   Each  problem  run  consists  of  two  or 
three  card  groups,  as  follows: 

1.  Title  card  (20A4).   One  card  only.   Card  must  be  present.   The 
identifying  information  from  the  title  card  is  printed  at  the  head 
of  the  output. 

2.  Namelist  card(s).   One  or  more  cards.   Details  are  given  in  the  next 
section. 


3.   Tabular  input  (optional).   This  input  group  is  contained  only  for 
the  first  run  and  for  those  ensuing  runs  where  NEWPRP=T.   If  no 
tabular  input  is  present,  then  blank  cards  must  not  be  inserted.   If 
tabular  input  is  present,  it  is  arranged  as  follows: 

First  card:   NCN,  NCP,  NEM,  NR,  NQG  (1013) .   These  are  the  number  of 
points  for  the  wall  thermal  conductivity,  wall  heat  capacity,  wall 
emissivity,  mass  pyrolysis  rate,  and  wall  heat  generation  rate, 
respectively.   The  number  of  points  may  be  0,  1,  or  greater  than  1. 
If  N=0,  then  the  previous  run  value  is  unchanged.   If  N=l,  then  it 
is  assumed  the  value  is  a  constant,  independent  of  temperature  or 
time.   If  N>1,  then  an  array  is  inputted. 

Ensuing  card(s):   These  are  in  the  format  (8F10.0)  and  arranged  in 
pairs  (independent,  dependent).   For  wall  thermal  properties,  tempera- 
ture is  the  independent  variable,  while  for  mass  pyrolysis  rate  it  is 
time.   The  order  is:   CNDA,  CPW,  EMSA,  RPX,  QGEN .   First  all  the  points 
(if  any)  for  CNDA  are  read  in,  four  pairs  per  card.   Then  a  new  card 
is  started  even  if  the  last  card  is  part-full,  and  CPW  is  read  in. 
The  process  is  continued  for  EMSA,  RPX,  QGEN.   No  blank  cards  may  be 
inserted.   If  N=l  for  an  array,  then  the  constant  value  is  entered 
in  columns  11-20. 

After  cards  for  one  run  are  finished,  the  cards  for  the  next  one  are 
stacked,  again  with  no  blank  cards. 

7.2   Name list  VARS 

For  all  non-tabular  data,  the  namelist  format  was  adopted.   This 
undeservedly  obscure  Fortran  feature  is  highly  advantageous  for  the  present 
application.   Its  features  include: 

•  Semi  free- format  input 

•  Variables  may  be  in  any  order 

•  Unneeded  variable  values  need  not  be  specified 

•  Variables  needed,  but  not  specified  in  current  run  are  automatically 
set  equal  to  the  prior  given  value. 

The  namelist  card(s)  must  contain  the  following  information:   the  first 

card  must  start  with  $VARS  in  columns  2-6,  then  a  space,  then  the  desired 

values,  separated  by  commas.  Input  may  be  continued  on  continuation  cards, 
each  of  which  must  have  columns  1-2  blank;  The  stream  is  terminated  by  a  $ 
after  the  last  variable. 

The  user  is  cautioned  to  check  the  input  carefully,  since  namelist 
format  provides  for  only  rudimentary  error  messages.   The  namelist  VARS 
values  are  written  to  file  5  when  read  in.   In  normal  operation  file  5  can 
be  rewound  or  discarded.   If  error  failures  occur,  however,  the  VARS  listing 
on  file  5  may  be  useful  in  determining  input  errors. 

Table  1  lists  all  the  variables  inputted  in  namelist  VARS. 

7.3   Modes  of  Program  Operation 

Time 

Three  possibilities  are  available:   complete  time-temperature  curve 
calculation,  calculation  of  steady-state  temperature  for  a  given  wall,  or  the 
calculation  of  a  steady-state  temperature  for  adiabatic  walls.   To  select 
adiabatic  walls,  set  ADIA=TRUE.   To  select  steady-state  solution  for  real  walls, 
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set  STEADY=TRUE.   To  obtain  complete  time-temperature  curve,  set  ADIA  and 
STEADY  both  FALSE.   Note  that  for  some  fuel  pyrolysis  conditions  below  not 
all  three  possibilities  are  available. 

Fuel  Pyrolysis 

The  following  modes  of  operations  are  available: 

1)  Pool  fire 

Set  PLFUEL=TRUE. 

a.  Time- temperature  curve,  for  given  ventilation  and  pool  area. 
Specify  SIZE,   Set  STOICH  and  EISCAN  both  FALSE. 

b.  Burning  conditions  at  steady  state  for  stoichiometric  pool  size, 
that  is,  determine  values  for  EITA=1.   Set  STOICH=TRUE.   Do  not 
input  SIZE.   Do  not  set  EISCAN=TRUE. 

c.  Burning  conditions  for  any  other  EITA.   Set  EISCAN=TRUE.   Specify 
EITA.   This  option  must  be  preceded  by  the  stoichiometric  problem 
(option  lb,  above) .   SIZE  input  is  not  used;  if  given,  the  value 
is  disregarded. 

2)  Wood  crib  fire.   This  is  the  default  option.   Set  FLSPEC,  PLFUEL, 
RPSPEC,  and  VTSPEC  all  FALSE. 

a.  Simple  stick  burning.   Must  specify  a  value  for  REGRES  greater 
than  zero. 

b.  Nilsson's  crib  formulas  for  crib  burning  in  three  possible 
regimes.   Specify  REGRES=0 .  (default).   Also  specify  SH. 

3)  Checking  option  when  tabular  input  pyrolysis  rates  are  given. 

Set  RPSPEC=TRUE.   Also  must  set  NEWPRP=TRUE  and  give  an  appropriate 
array  of  RPX. 

4)  Pessimization  over  ventilation.   Set  FLSPEC=--TRUE .   Window  width  is 
automatically  adjusted,  but  is  no  greater  than  determined  by  the 
inputted  value  of  AWDOW/HWDOW.   Program  stops  when  fuel  is  exhausted. 

a.  Simple  stick  burning.   Must  specify  a  value  for  REGRES  greater 
than  zero  and  set  PLFUEL=FALSE. 

b.  Nilsson's  crib  formulas  for  crib  burning  in  three  possible 
regimes.   Set  PFLUEL=FALSE  and  REGRES=0.   Also  specify  SH . 

c.  Pool  burning.   Set  PLFUEL=TRUE. 

5)  Pessimization  over  fuel  pyrolysis  rate.   Set  VTSPEC=TRUE.   Fuel 
pyrolysis  rate  is  automatically  adjusted  for  pessimal  burning 
conditions . 

8.   FILES  USED 

The  Fortran  file  logical  units  must  be  declared  as  follows: 

File  1  --  Input 

File  2  --  Output  (echoed  input  and  main  calculated  variables) 

File  3  --  Output  (heat  balance  and  mass  fractions) 

File  4  --  Output  (intermediate  tracing  output  -  used  only  if  KTRACE=1) 

File  5  --  Output  (listing  of  namelist  VARS  contents) . 
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File  5  can  be  arranged  to  be  rewound  after  each  problem  so  that  it  will 
contain  data  only  in  case  of  error  failure. 

9.   IMPLEMENTATION 

Program  C0MPF2  has  been  successfully  implemented  on  a  UNIVAC  110  8  computer. 
The  predecessor  program,  COMPF,  was  run  on  a  CDC  6400  computer.   The  program 
uses,  as  much  as  possible,  only  standard  Fortran  expressions.   Minor  unavoidable 
implementation  differences  exist,  however,  in  commands  associated  with  file 
usage. 

10.   LIST  OF  VARIABLES 
Table  2  gives  a  list  of  all  the  major  problem  variables. 
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Table  1.   Variables  specified  in  the  input  Namelist  VARS 


Default  Values 


First 
run 


Following 
run 


Information 


ADIA 


FALSE 


pv 


AFLOOR 

none 

pv 

AWALL 

none 

pv 

AWDOW 

none 

pv 

BPF 

none 

pv 

CD 

0.68 

pv 

CFLPC 

44.4 

pv 

CPPYR(2) 

CPN2 

pv 

CVGROS 

none 

pv 

DENSW 

none 

pv 

DHP 

none 

pv 

DTI  ME 

none 

pv 

EF 

0.9 

pv 

EI SCAN 

FALSE 

FALSE 

EITA 


none 


pv 


FLOAD 

none 

pv 

FLSPEC 

FALSE 

FALSE 

HFLPC 

5.4 

pv 

HWDOW 

none 

pv 

I  RUN 

1 

sequential 

IX 

10 

pv 

KTRACE 

0 

0 

MTIME 

360. 

pv 

MWPYR 

28.97 

pv 

NEWPLT 

FALSE 

FALSE 

NEWPRP 

TRUE 

FALSE 

NFLPC 

0. 

pv 

OFLPC 

0. 

pv 

PLFUEL 

FALSE 

FALSE 

PLOT 

FALSE 

pv 

PNCH 

FALSE 

pv 

PRNT 

none 

pv 

REGRES 

0. 

pv 

RPSPEC 

FALSE 

FALSE 

SH 

0 

pv 

SHAPE 

2. 

pv 

SIZE 

none 

pv 

STEADY 

FALSE 

FALSE 

STOICH 

FALSE 

FALSE 

if  true,  walls  area  adiabatic  and  only 

steady-state  solution  is  sought 
area  of  floor  (m^) 

gross  area  of  walls  and  ceiling  (m^) 
area  of  window  (m^) 
maximum  fraction  of  pyrolyzed  fuel  burned 

to  be  ^1.0 
discharge  coefficient 

percent,  by  weight,  of  carbon  in  fuel 
heat  capacity  of  pyrolysis  gases  (J/kg-K) 
upper  calorific  value  for  dry  fuel  (J/kg) 
wall  density  (kg/in-^) 

total  heat  of  gasification  for  fuel  (J/kg) 
increment  of  time  step  (s) 
gas  emissivity,  assumed  gray 
if  true,  solve  steady-state  problem  in 

POOL  for  a  given  EITA 
normalized  air-fuel  parameter  for  pool 

burning 
fuel  load  (kg/m^  floor  area) 
if  true,  pessimize  ventilation  for  a 

specified  pyrolysis  rate 
percent  of  hydrogen,  by  weight,  in  fuel 
window  height  (m) 
run  problem  number 

number  of  wall  slices,  to  be  _<  10 
print  intermediate  output  if  =1 
maximum  fire  time  (s) 
molecular  weight  of  pyrolysis  gases  (g/g- 

mole) 
if  true,  start  new  plot  frame 
if  true,  new  data  arrays  will  be  given 
percent  of  nitrogen,  by  weight,  in  fuel 
percent  of  oxygen,  by  weight,  in  fuel 
if  true,  fuel  is  a  pool  fire 
if  true,  plot  time-temperature  curve 
if  true,  punch  time-temperature  curve 
interval  at  which  results  are  to  be 

printed  (s) 
rate  of  fuel  regression  (m/s) 
if  true,  use  tabular  input  fuel  pyrolysis 
ratio  of  clear  spacing  between  sticks/crib 

height  for  crib 
shape  factor  in  pyrolysis  equation  for 

wood  sticks 
for  cribs:   smallest  dimension  of 

stick  (m) 
for  pools:   pool  area  (m^) 
if  true,  only  steady-state  solution  is  to 

be  sought 
if  true,  EITA=1  solution  is  sought  in 

POOL 
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Table  1 

TBOILC 

0. 

pv 

THICKW 

none 

pv 

TINPT 

0. 

0. 

VTSPEC 

FALSE 

FALSE 

WFLPC 

0. 

pv 

Note: 

(continued) 

fuel  vaporization  temperature  for  pools  (C) 

wall  thickness  (m) 

optional  input  iteration  gas  temperature 

(K) 
if  true,  pessimize  pyrolysis  rate  for  a 

specified  ventilation 
percent  of  water,  by  weight,  in  fuel 


pv  =  previous  value 
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Table  2.   List  of  variables 


AD  I A 

AFLOOR 

AWALL 

AWALLN 

AWDOW 

BPF 

BWDOW 

BWORST 

C 

CD 

CFLDC 

CND 

CNDA 

CNG 

CNL 

CNV 

CPA 

CPCO 

CPC02 

CPH2 

CPH20 

CPN2 

CP02 

CPPYR 

CPW 

CVGROS 

CVNET 

DENF 

DENSA 

DENSW 

DENU 

DERIVl 

DERIV2 

DIF 

DTGAS 

DTIME 

DX 

EF 

EI SCAN 

EITA 

EMS 

EMSA 

FC 

FLOAD 

FLREM 

FLSPEC 

FUELPC 

Fl 

F2 

G 

II 

HCP 

HF 

HFLPC 

HU 

HIN 

HRATIO 

HWDOW 

ILINE 


true  if  walls  are  adiabatic 

area  of  floor  (m^) 

gross  area  of  walls  and  ceiling  (m^) 

AWALL  minus  window  area 

area  of  window  (m^) 

maximum  fraction  of  pyrolyzed  fuel  burned 

width  of  window  (m) 

window  width  (.LE. BWDOW)  which  maximizes  gas  temperatures  (m) 

moles  of  carbon  in  fuel  (mole/kg  fuel) 

discharge  coefficient 

percent  of  carbon,  by  weight,  in  fuel 

conductivity  of  a  given  wall  slice  (W/m-K) 

conductivity  of  the  wall,  as  a  function  of  temperature  (W/m-K) 

average  conductivity,  next  to  higher  numbered  slice 

average  conductivity,  next  to  lower  numbered  slice 

numerical  factor  in  heat  transfer  coefficient 

heat  capacity  of  ambient  air  (J/kg-K) 

heat  capacity  of  CO,  as  a  function  of  temperature  (J/kg-K) 

heat  capacity  of  CO2/  as  a  function  of  temperature  (J/kg-K) 

heat  capacity  of  H2,  as  a  function  of  temperature  (J/kg-K) 

heat  capacity  of  H2O,  as  a  function  of  temperature  (J/kg-K) 

heat  capacity  of  N2,  as  a  function  of  temperature  (J/kg-K) 

heat  capacity  of  O2/  as  a  function  of  temperature  (J/kg-K) 

heat  capacity  of  pyrolysis  gases,  as  a  function  of  temperature 

(J/kg-K) 
wall  heat  capacity,  as  a  function  of  temperature  (J/kg-K) 
upper  calorific  value  for  dry  fuel  (J/kg) 
lower  calorific  value  for  moist  fuel  (J/kg) 
Biot  Number/2--f ire  side 
ambient  air  density  (kg/m^) 
wall  density  (kg/m^) 
Biot  Number/2--unexposed  side 

current  derivative  of  heat  balance  remainder  (W/K) 
previous  derivative  of  heat  balance  remainder  (W/K) 
temperature  error  in  iteration  (K) 
increment  in  gas  temperature  (K) 
time  increment  (s) 
wall  thickness  increment  (m) 
effective  flame  grey  body  emissivity 
true  if  seeking  contant  EITA^^l  solution 
dimensionless  air/fuel  parameter  for  pool  burning 
computed  wall  emissivity  for  parallel  plane  problem 
wall  emissivity,  as  a  function  of  temperature 
true  if  in  fuel  control 
fuel  load  (kg/m^  of  floor  area) 
mass  of  fuel  remaining  at  a  given  time  (kg) 

true  if  fuel  pyrolysis  rate  is  fixed  and  ventilation  pessimized 
percent  of  original  fuel  supply  still  remaining 
current  heat  balance  error  (W) 
previous  heat  balance  error  (W) 
acceleration  of  gravity  (m/s^) 
moles  of  hydrogen  in  fuel  (mole/kg  fuel) 
variable  in  solving  differential  equation 
effective  heat  transfer  coefficient,  fire  side  (W/m^-K) 
percent  of  hydrogen,  by  weight,  in  fuel 

effective  heat  transfer  coefficient,  unexposed  side  (W/m^-K) 
height  of  neutral  plane  (m) 

fractional  height  of  neutral  plan  above  window  bottom 
height  of  window  (r") 
line  number 
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Table  2, 


(continued) 


IPG 
I  RUN 
IX 
IXC 

IXL 

J 

JM 

JPRINT 

K 

KD 

KITER 

KNTRL 

KTRACE 

MTIME 

MWIN 

MWOUT 

MWPYR 

N 

NCND 

NCPW 

NEMS 

NEWPLT 

NEWPRP 

NFLPC 

NQGEN 

NRP 

0 

OFLPC 

OPENF 

PLFUEL 

PLOT 

PNCH 

PRNT 

QCONW 

QFIRE 

QFLOW 

QFUEL 

QGEN 

QRADO 

QFLOW 

QFUEL 

QGEN 

QRADO 

QRADW 

QWLSUM 

R 

RO 

RC 

REGRES 

RMA 

RMF 

RP 

RPSPEC 

RPX 

SCAN 

SH 

SHAPE 

SIGMA 

SIZE 

SIZEl 


page  number 

run  number 

number  of  wall  slices 

number  of  middle  slice 

number  of  penultimate  slice 

number  of  current  time  step 

maximum  number  of  time  steps 

output  to  be  printed  every  JPRINT  time  steps 

number  of  trial  iterations  at  any  given  time  step 

number  of  iterations  to  converge  differential  equation 

equals  0  for  normal  operation,  equals  1  for  convergence  failure 

parameter  indicating  exit  status 

print  intermediate  tracing  output  if  KTRACE-1 

maximum  time  for  fire  simulation  (s) 

molecular  weight  of  ambient  air  (g/g-mole) 

molecular  weight  of  exhaust  gases  (g/g-mole) 

molecular  weight  of  pyrolysis  gases  (g/g-mole) 

moles  of  nitrogen  in  fuel  (mole/kg  fuel) 

number  of  points  in  CNDA  table 

number  of  points  in  CPW  table 

number  of  points  in  EMSA  table 

true  if  start  new  plot  frame  (not  overlay  previous  one) 

true  if  read  in  new  set  of  tabular  data 

percent  of  nitrogen,  by  weight,  in  fuel 

number  of  points  in  QGEN  table 

number  of  points  in  RPX  table 

moles  of  oxygen  in  fuel  (mole/kg  fuel) 

percent  of  oxygen  by  weight,  in  fuel 

opening  factor  ration  (m    ) 

true  if  pool  fire  configuration 

true  if  plot  time-temperature  curve 

true  if  punch  time-temperature  curve 

number  of  times  per  second  output  is  to  be  printed 

heat  transfered  to  walls  by  convection  (W) 

heat  generated  by  combustion  (W) 

net  flow  enthalpy  (exhaust  minus  inflow)  (W) 

heat  lost  in  heating  up  unburned  fuel  fraction  (W) 

wall  heat  generation,  as  a  function  of  temp  (W/m^) 

heat  radiated  out  the  window  (W) 

new  flow  enthalpy  (exhaust  minus  inflow)  (W) 

heat  lost  in  heating  up  unburned  fuel  fraction  (W) 

wall  heat  generation,  as  a  function  of  temp  (W/m^) 

heat  radiated  out  the  window  (W) 

heat  transfered  to  walls  by  radiation  (W) 

total  heat  removed  from  compartment  and  passing  into  the  walls  (J) 

stoichiometric  air/fuel  mass  ratio 

stoichiometric  oxygen/fuel  mass  ratio 

rate  of  burning  (kg/s) 

rate  of  fuel  surface  regression  (m/s) 

mass  inflow  rate  of  air  (kg/s) 

mass  outflow  rate  of  hot  gases  (kg/s) 

rate  of  pyrolysis  (kg/s) 

true  if  rate  of  pyrolysis  is  prescribed  as  input 

rate  of  pyrolysis,  as  a  function  of  time  (kg/s) 

true  if  search  for  solution  by  scanning  temperatures 

ratio  of  clear  spacing  between  sticks  to  crib  height 

contant  indicating  shape  of  fuel  sticks 

Stef an-Boltzmann  constant  (w/m^-K**) 

thickness  of  crib  sticks  (m) 

area  of  pool  (m^) 

pool  area  for  EITA=1  condition  (m^) 
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Table  2.   (continued) 

STEADY  true  if  only  steady-state  calculation  to  be  made 

STOICH  true  if  pool  fire  and  EITA=1 

TAMB  ambient  temperature  (K) 

TGAS  gas  temperature  (K) 

TGASl  previous  value  of  TGAS  (K) 

TGAS2  previous  value  of  TGASl  (K) 

TGASC  gas  temperature  (C) 

TGASN  closest  gas  temperature,  lower  than  true  (K) 

TGASP  closest  gas  temperature,  higher  than  true  (K) 

TGOLD  value  of  TGAS  from  prior  time  step  (K) 

THICKW  wall  thickness  (m) 

TINPT  input  trial  starting  gas  temperature  (K) 

TITLE  title  of  this  run 

TSF  wall  surface  temperature,  fire  side  (K) 

TSU  wall  surface  temperature,  unexposed  side  (K) 

TTIME  total  time  (s) 

Tl  old  wall  temperature  profile  (K) 

T2  new  wall  temperature  profile  (K) 

T2C  wall  temperature  profile  (C) 

VAVGIN  average  inflow  velocity  (m/s) 

VTSPEC  true  if  ventilation  is  fixed  and  pyrolysis  rate  pessimized 

W  moles  of  water  in  fuel  (mole/kg  fuel) 

WA  format  constant 

WB  format  constant 

WFLPC  percent  of  water,  by  weight,  in  fuel 

WTFUEL  initial  total  mass  of  fuel  (kg) 

YC02  mass  fraction  of  CO2  in  outflow 

YH20  mass  fraction  of  H2O  in  outflow 

YN2  mass  fraction  of  N2  in  outflow 

Y02  mass  fraction  of  O2  in  outflow 

YPYR  mass  fraction  of  pyrolysates  in  outflow 
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Figure  1.   Flow  chart  for  main  program  C0MPF2 
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APPENDIX  A  —  SAMPLE  PROBLEMS 

Given  below  is  a  set  of  ten  concatenated  input  problems.   Each  problem 
is  intended  to  test  out  a  subroutine  or  other  feature  of  the  program.   The 
output  for  the  problems  is  given  following  the  input. 


TEST  PROGRAM  FOR  POOL  FIRE.  STEADY  STATE*  EITA=1.0 
SVARS  AFLOOR*20«  .  AWALL=eO«  .  AWDGIK=4#  .  BPF  =  0  .7  .CD=0  .  68.  CFLPC=85  .  7  . 

CVGR0S=46«5E6.DENSW=790« .DHP=2.4E6.EITA=1 .O.FLOAD=20e .HFLPC=14o3. 
HI*DOW=1.5.0FLPC=O.O.P1_FUEL=T.STOICH=T.TBOILC=390«. 
THICK W=0« 038. WFLPC=0«0$ 
001000001 

0.17 
0*5 
TEST  PROGRAM  FOR  POOL  FIRE.  STEADY  STATE.  EITA=0«Oi 

SVARS  EISCAN=T.EITA=0«01 .PLFUEL=T$ 
TEST  PROGRAM  WITH  DELIBERATE  ERROR  TO  CHECK  KTRACE  OPERATION 

SVARS  £ISCAN=T. PLFUEL=T.TBOILC=2000» $ 
TEST  PROGRAM  FOR  POOL  FIRE.  TRANSIENT  CASE.  SIZE=7,5  M2 
SVARS  OTIME=60« .MTIME=3600«.NEWPRP=T.  PLFUEL=T.  PRNT=60. . SI ZE=7. 5. 
TBOILC=390«S 
000001 

840* 
TEST  PROGRAM  FOR  WOOD  CRIB  FIRE.  REGRES  SPECIFIED 
SVARS  CFLPC=44«4.CVGROS=i8«8E6.FLOAD=:S0«0«HFLPC=5«4.OFLPC=38.2. 
REGRES-1«5E-5.SHAPE=2«0.SIZE=0.05.WFLPC=12.0S 
TEST  PROGRAM  1  FOR  WOOD  CRIB  FIRE.  NILSSON'S  FORMULAS 

SVARS  REGRES=0*0.SH=0«10S 
TEST  PROGRAM  2  FOR  WOOD  CRIB  FIRE.  NILSSON'S  FORMULAS 

SVARS  FLOAD=20«.SH=0«20S 
TEST  PROGRAM  FOR  PVTFIX  ROUTINE.  VARIABLE  WALL  PROPERTIES 

SVARS  NEWPRP=T.VTSPEC=TS 
004006 

£73*       0*2i       372*       0«2i       373«        0«&6 
273.       1090.      372»       1090.      373.       47300. 
3e4«        5000*       413.        5000*       4&4«        840* 
TEST  PROGRAM  FOR  PFLFIX  ROUTINE.  POOL  OPTION 
SVARS  AWDOW=10«.CFLPC=e5e7.CVGROS=46«5E6.FLSPEC=T.HFLPC=14.3.NEWPRP=T. 
OFLPC=0«0.PLFUEL=T.SIZE=5.0.WFLPC=0.0$ 
OOiOOl 

0*17 
640* 
TEST  PROGRAM  FOR  RPFIX  ROUTINE 

SV^.RS  MTIME=1903«  .  N£WPRP=T  .  RPSPEC=TS 
000000000003 
0*0        0«i2       S20«       0*12       i2ie        0.25 


a  073. 

0«26 

363* 

47300 

i073« 

840* 
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TEST  PROGRAM  MITH  DECISERATE  ERROR  TO  CHECK  KTRACE  OPERATION 
GEOMETRY  AND  VENTILATION 


C0MPF2  VERSION  1.1  -  RUN  NO.     3 


WALL  SURFACE  AREA  =       80.0  M2 
FLOOR  AREA  =      20.00  M2 
WINDOW  HEIGHT  =       1 . 50  M 

AREA  =        4.00  M2 
OPENING  FACTOR  =     4.899  M2 . S 
DISCHARGE  COEFF.=   .68 


-FUEL  LOAD  PROPERTIES 


FIRE  LOAD  PER  FLOOR  AREA  =   20.0  KG/M2 

TOTAL  ENTHALPY  OF  PYROLYSIS=      2.40*06  J/KG 

BOILING  TEMPERATURE=2000.  OEG  C 


PERCENT  8Y  WEIGHT 

PERCENT 

PERCENT 

PERCENT 

PERCENT 


FUEL  COMPOSITION 

CARBON    =  85. 

HYDROGEN  =  14. 

OXYGEN     = 

NITROGEN  =     . 

WATER      = 

R  =  14.78 

RO-   3.43 

HEAT  OF  COMBUSTION  OF  DRY  FUEL  = 

LOWER  ACTUAL  HEAT  OF  COMBUSTION  = 

MOLECULAR  WEIGHT  OF  UNBURNT  PYROLYSATES  =    28.97 

CP  OF  PYROLYSIS  GAS  =  (   .1I27»TGAS  ♦    1010.)  J/KG- 

MAXIMUM  FRACTION  OF  PYROLYSATES  BURNED  =     .70 

GREY-GAS  FLAME  EMISSIVITY  =    .900 

FUEL  AREA=      174.07  M2 


46.50+06  J/KG 
43.36+06  J/KG 


WALL  THERMAL  PROPERTIES- 


THICKNESS 
DENSITY  = 


.038    M 
790.     KG/M3 


THERMAL    CONDUCTIVITY 
EMISSIVITY    =        .50 


.170  W/M-K 


THERMOPLASTIC  POOL  FIRE 


PAGE  NO.    1      RUN  NO. 


TEMP 
GAS.C 


WALL  TEMPS 


RC 
KG/S 


EXC .PYR. 
KG/S 


FUEL      AIR  IN 
PCT       KG/S 


N.P.  VELOCITY  MOL.WT  FUEL 

M/S  CNTRL 


ITERATION    FAILURE  — 

TGAS=  1890.00 


TGAS.LT.TBOIL  TGAS=       1890.0  GO    TO    NEXT    CASE 


RUN    NO.  3 

TGASl  TGAS2 

1800.00  .00 

1810.00  1800.00 

1900.00  1810.00 


.43*08 
.40+08 
.10+08 


2.43+08 
2.40+08 


DERI  VI  K  KO  KH 

1.35+05  1  3  0 

-2.93+05  2  3  0 

-3. 33+05  3  3  0 


1730.86  1791.43 
1740.49  1801.51 
1827.13       1892.16 


4.028+05  -1. 
2.891+04  -1. 
2.405+03    -9, 


009+08 
008+08 
183+07 


4.052+05-59.938 
4.083+0  5-59.068 
4.363+05-50.558 


.009 
.001 
.000 
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TEST  PROGRAM  FOR  POOL  FIRE,  TRANSIENT  CASE.  SIZE=7.5  M2 

GEOMETRY  AND  VENTILATION 

WALL  SURFACE  AREA  =       80.0  M2 
FLOOR  AREA  =     20.00  M2 
WINDOW  HEIGHT  =       1.50  M 

AREA  =       4.00  M2 
OPENING  FACTOR  =    4.899  M2.5 
DISCHARGE  COEFF.=   .68 


— FUEL  LOAD  PROPERTIES 

FIRE  LOAD  PER  FLOOR  AREA  =   20.0  KG/M2 
TOTAL  ENTHALPY  OF  PYROLYSIS=     2.40+06  J/KG 
BOILING  TEMPERATURE=  390.  DEG  C 

FUEL  COMPOSITION 

CARBON    =  85.7  PERCENT  BY  WEIGHT 

HYDROGEN  =  14.3  PERCENT 

OXYGEN    =    .0  PERCENT 

NITROGEN  =    .0  PERCENT 

WATER     =    .0  PERCENT 

R  =  14.78 

R0=   3.43 

HEAT  OF  COMBUSTION  OF  DRY  FUEL  =    46.50+06  J/KG 

LOWER  ACTUAL  HEAT  OF  COMBUSTION  =    43.36+06  J/KG 

MOLECULAR  WEIGHT  OF  UNBURNT  PYROLYSATES  =   28.97 

CP  OF  PYROLYSIS  GAS  =  (  .1127*TGAS  +   1010.)  J/KG-K 

MAXIMUM  FRACTION  OF  PYROLYSATES  BURNED  =     .70 

GREY-GAS  FLAME  EMISSIVITY  =   .900 

FUEL  AREA=        7.50  M2 

—-WALL     THERMAL    PROPERTIES r- 

THICKNESS    =       .038    M 
DENSITY    =  790.    KG/M3 

THERMAL    CONDUCTIVITY    =  .170     W/M-K 

HEAT    CAPACITY    =  840.     J/KG-K, 

EMISSIVITY    =       .50 


I 
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THERMOPLASTIC  POOL  FIRE 


33 
34 


TIME 

S 


60. 

120. 

180. 

240. 

300. 

360. 

420. 

480. 

540. 

600. 

660. 

720. 

780. 

84  0. 

900. 

960. 
1020. 
1080^ 
1  140. 
1200. 
1260. 
1  320. 
1  380. 
1440. 
1500. 
1560. 
1620. 
1680. 
1740. 
1  800. 
1  86  0. 
1920. 
1980. 
2040. 
2100. 


TEMP 
GASiC 

739. 

826. 
852. 
870. 
882. 
891. 
898. 
903. 
907. 
91  1. 
914. 
917. 
919. 
921. 
923. 
925. 
926. 
358. 
255. 
211. 
184. 
164. 
149. 
138. 
129. 
121. 
1  14. 
108. 
103. 

99. 

95. 


82. 
80. 


WALL    TEMPS 

RP 

RC 

EXC .PVR. 

FUEL 

AIR   n 

C 

KG/S 

KG/S 

KG/S 

PCT 

KG/S 

513. 

26. 

25. 

•  306 

.108 

.  198 

95 

.4 

2.29 

712. 

34. 

25. 

.36  2 

.  106 

.  256 

90 

.0 

2.23 

763. 

53. 

25. 

.381 

.105 

.277 

84 

.3 

2.21 

798. 

83. 

26. 

.396 

.  104 

.292 

78 

.3 

2.20 

820. 

117. 

27. 

.406 

.104 

.302 

72 

.2 

2.19 

836. 

151. 

30. 

.413 

.103 

.310 

66 

.0 

2.18 

847. 

183. 

34. 

.419 

.103 

.316 

59 

.7 

2.18 

857. 

213. 

39. 

.424 

.103 

.321 

53 

•  4 

2.17 

86  4. 

240. 

46. 

.428 

.  103 

.325 

47 

.0 

2.17 

871  . 

265. 

53. 

.431 

.103 

.328 

40 

.5 

2.17 

876. 

288. 

61. 

.434 

.102 

.331 

34 

.0 

2.16 

881. 

309. 

69. 

.436 

.102 

.334 

27 

.5 

2.16 

885. 

329. 

77. 

.4  39 

.102 

.336 

20 

.9 

2.16 

888. 

347. 

85. 

•  440 

.102 

.338 

14 

.3 

2.16 

892. 

364. 

93. 

.442 

.102 

•  340 

7 

.6 

2.16 

394. 

380. 

101. 

.444 

.102 

.342 

1 

.0 

2.15 

897. 

394. 

108. 

.445 

.102 

.34  3 

.0 

2.15 

534. 

407. 

115. 

.000 

•  000 

.000 

.0 

2.47 

409. 

414. 

122. 

.000 

•  000 

.000 

.0 

2.39 

350. 

412. 

128. 

.000 

•  000 

.000 

.0 

2.32 

310. 

402. 

133. 

.000 

•  000 

.000 

.0 

2.26 

280. 

387. 

137. 

.000 

•  000 

.000 

.0 

2.21 

257. 

371. 

141. 

.000 

•  000 

.000 

.0 

2.16 

238. 

355. 

143. 

.000 

.000 

.000 

.  0 

2.12 

222. 

339. 

145. 

•  000 

.000 

.000 

.0 

2.09 

209. 

324. 

145. 

.003 

.000 

.000 

•  0 

2.05 

196. 

310. 

145. 

.000 

.000 

.000 

.0 

2.02 

188. 

296. 

144. 

.000 

.000 

.000 

.0 

1  .99 

179. 

284. 

142. 

.000 

.000 

.000 

.0 

1.97 

171. 

272. 

140. 

.000 

.000 

.000 

.0 

1.95 

164. 

261. 

137. 

.000 

.000 

.000 

.0 

1.92 

157. 

251. 

135. 

.000 

.000 

.000 

.0 

1  .90 

151. 

24  1. 

132. 

.000 

.000 

.000 

.0 

1.88 

146. 

232. 

129. 

.000 

.000 

.000 

.0 

1.86 

141. 

223. 

126. 

.000 

.000 

.000 

.0 

1.85 

136. 

215. 

123. 

.000 

,000 

.000 

.0 

1.83 

60. 

120. 

180. 

240. 

300. 

360. 

420. 

480. 

540. 

600. 

660. 

720. 

780. 

840. 

900. 

960. 
1020. 
1080. 
1140. 
1200. 
1260. 
1320. 
1380^ 
1440^ 
ISOO. 
1560. 
1620. 
1680. 
1740. 
ISOO. 
1860. 
1920. 
1980. 
2040. 
2100. 


THERMOPLASTIC  POOL  FIRE 


HEAT  BALANCE 


GAS    FLO* 

WND    RAD 

WALL    C> 

PCT 

PCT 

PCT 

44.886 

5.322 

1  1.092 

51.859 

7.174 

4.559 

54.089 

7.950 

3.297 

55.726 

8.546 

2.523 

56.787 

8.943 

2.091 

57.S86 

9.250 

1  .787 

58.196 

9.486 

1.574 

58.693 

9.681 

1.41  1 

59.095 

9.840 

1  .283 

59.440 

9.978 

1  .180 

59.737 

10.097 

1.094 

59.996 

10.201 

1.022 

6  0.2  24 

10.294 

.96  0 

60.427 

10.376 

.907 

60.610 

10.451 

.860 

60.7  76 

10.518 

.818 

60.927 

10.580 

.781 

96.257 

3.743 

-40.770 

97.357 

2.643 

-52.654 

97.717 

2.283 

-58.391 

97.914 

2.086 

-62.064 

98.036 

1  .964 

-64.780 

98.120 

1.880 

-67.042 

98.180 

1.820 

-68.677 

98.223 

1.777 

-69.805 

98.258 

1  .742 

-71 .093 

98.285 

1.715 

-71 .919 

98.306 

1.694 

-72.779 

98.324 

1.676 

-73.334 

98.338 

1.662 

-73.959 

98.349 

1.651 

-74.362 

98.359 

1.641 

-74.862 

98.367 

1.633 

-75.260 

98.374 

1.626 

-75.620 

98.379 

1.621 

-75.777 

98.384 

1.616 

-76.079 

WALL    RAO 
PCT 

28.91 0 

22.952 

20.091 

17.707 

16.097 

14.79  4 

13.785 

12.954 

12.252 

1  I .659 

1 1.144 

10.691 

10.29  0 

9.929 

9.604 

9.308 

9.037 

-59.079 

-47.486 

-41.795 

-37.859 

-35.129 

-33.129 

-31  .50  5 

-30.128 

-29.103 

-28.159 

-27.407 

-26.712 

-26.146 

-25.603 

-25.162 

-24.756 

-24.394 

-24.028 

-23.735 


4.694*06 

4.582+06 

4,543+06 

4.514+06 

4.496+06 

4. 481+06 

4.470+06 

4.461+06 

4. 453+06 

4. 447+06 

4. 442+06 

4,437+06 

4.433+06 

4.429+06 

4.426+06 

4.423+06 

4.4  20+06 

0.000 

0.000 

0.  000 

0.  000 

0.000 

0.  000 

0.000 

0.  0  00 

0.000 

0,  000 

0,000 

0.000 

0.  000 

0.  000 

0.000 

0.  000 

3.000 

0,  000 

0.000 


a-WALL 
SUM 


1  .127+08 
1  .883  +  08 
2.520+08 
3.068+08 
3,559+08 
4.005+08 
4.417+08 
4.801+08 
5.163+08 
5.505+08 
5.832+08 
6.143+08 
6.443+08 
6.731+08 
7.008+08 
7.277+08 
7.537+08 
7.537+08 
7.537+08 
7.537+08 
7.537+08 
7.537+08 
7.537+08 
7.537+08 
7.537+08 
7.537+08 
7.637+08 
7.537+08 
7.537+08 
7.537+08 
7.537+08 
7.537+08 
7.537+08 
7.537+08 
7.537+08 
7.537+08 


N.P. 

VELOCITY 

MOL. WT 

FUEL 

M/S 

CNTRL 

.38 

1.86 

29.87 

F 

.37 

1.87 

29.85 

F 

.37 

1.87 

29.85 

F 

.37 

1.87 

29.84 

F 

.36 

1.87 

29.84 

F 

.36 

1.87 

29.83 

F 

.36 

1.87 

29.83 

F 

.36 

1.87 

29.83 

F 

.36 

1.87 

29.83 

F 

•  36 

1.87 

29.83 

F 

•  36 

1.87 

29.82 

F 

•  36 

1.87 

29.82 

F 

.36 

1.87 

29.82 

F 

.36 

1.87 

29.82 

F 

.36 

1.87 

29,82 

F 

.36 

1.87 

29.82 

F 

.36 

1.87 

29.82 

F 

.43 

1.78 

26.90 

T 

.45 

1.66 

26.95 

T 

.45 

1.59 

26.94 

T 

.46 

1.54 

26.90 

T 

.46 

1.49 

26.87 

T 

.46 

1.45 

26.84 

T 

.47 

1  .42 

26.80 

T 

.47 

1.39 

26.78 

T 

.47 

1,36 

26.75 

T 

.47 

1,34 

26.73 

T 

.47 

1  .31 

26.70 

T 

.47 

1.30 

26.67 

T 

.47 

1.28 

26.64 

T 

.48 

1.26 

26,62 

T 

.48 

1.24 

26,6  0 

T 

.48 

1.23 

26,58 

T 

.48 

1.22 

26,56 

T 

.48 

1  .20 

26.54 

T 

.48 

1.19 

26.52 

T 

PAGE     NO. 

1               RUN    NO. 

y02 

rN2 

YC02 

YH20 

YPYR 

PCT 

PCT 

PCT 

PCT 

PCT 

MASS 

MASS 

MASS 

MASS 

MASS 

.063 

.679 

.  131 

.054 

.076 

.058 

.663 

.  128 

.052 

.099 

.358 

.657 

.  127 

.052 

.  107 

.357 

.652 

.  126 

.052 

.113 

.057 

.650 

.126 

.051 

.116 

.057 

.647 

.  125 

,051 

.  1  19 

.057 

.646 

.  125 

,051 

.  122 

.057 

.644 

.  125 

.051 

.124 

.0  56 

.643 

.124 

.051 

.125 

.0  56 

.642 

.  124 

.051 

.126 

.0  56 

•  64  1 

.124 

.051 

.128 

.056 

.64  1 

.  124 

.051 

.129 

.056 

•  640 

.  124 

.051 

.130 

.056 

.639 

,124 

.051 

.130 

.0  56 

.639 

•  123 

.051 

.131 

.0  56 

.638 

.  123 

.051 

.  132 

.0  56 

.638 

.  123 

,050 

.  132 

.230 

.770 

.000 

.000 

.330 

.230 

.770 

.000 

.000 

.000 

.230 

.770 

.000 

.000 

.300 

.230 

.770 

.033 

.300 

.030 

.230 

.770 

.000 

.000 

.000 

.230 

.770 

.000 

.003 

.000 

.230 

.770 

.030 

.330 

.000 

.230 

.770 

.000 

.000 

.000 

.230 

.770 

,000 

.000 

.000 

.230 

.770 

.000 

.003 

.000 

.230 

.773 

•  300 

.000 

.000 

.230 

.770 

.300 

.000 

.300 

.230 

•  770 

.000 

.000 

.000 

.230 

.770 

.000 

.000 

.000 

.230 

.770 

.000 

.000 

.000 

.230 

.770 

.030 

.000 

.000 

.230 

.770 

,000 

.000 

.000 

.230 

.770 

.000 

.000 

.000 

.230 

.770 

.000 

.000 

.000 
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TEST  PROGRAM  FJDR  WOOD  CRIB  FIRE,  REGRES  SPECIFIED 

GEOMETRY  AND  VENTILATION 

WALL  SURFACE  AREA  =       80.0  M2 
FLOOR  AREA  =     20,00  M2 
WINDOW  HEIGHT  =       1.50  M 

AREA  =       4.00  M2 
OPENING  FACTOR  =    4.399  M2.5 
DISCHARGE  COEFF.=   .68 

FUEL  LOAD  PROPERTIES 

FIRE  LOAD  PER  FLOOR  AREA  =   10.0  KG/M2 

FUEL  COMPOSITION 

CARBON    =  44.4  PERCENT  BY  WEIGHT 

HYDROGEN  =   5.4  PERCENT 

OXYGEN    =  38.2  PERCENT 

NITROGEN  =    .0  PERCENT 

WATER     =  12.0  PERCENT 

R  =   5.32 

R0=   1.23 

HEAT  OF  COMBUSTION  OF  DRY  FUEL  =     18.80+06  J/KG 

LOWER  ACTUAL  HEAT  OF  COMBUSTION  =     15.07+06  J/KG 

MOLECULAR  WEIGHT  OF  UNBURNT  PYROLYSATES  =   28.97 

CP  OF  PYROLYSIS  GAS  =  (  .1127«TGAS  +   1010.)  J/KG-K 

MAXIMUM  FRACTION  OF  PYROLYSATES  BURNED  =     .70 

GREY-GAS  FLAME  EMISSIVITY  =   .900 

RATE  OF  REGRESSION  =   15.00-06  M/S 

FUEL  DIMENSION  =   .050  M 

SHAPE  FACTOR  =  2.00 

WALL  THERMAL  PROPERTIES 

THICKNESS  =   .038  M 
DENSITY  =    790.  KG/M3 

THERMAL  CONDUCTIVITY  =      .170  W/M-K 

HEAT  CAPACITY  =      840.  J/KG-K 

EMISSIVITY  =   .50 
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CRIB    FIRE 

TIME 

TEMP 

WALL    TEMPS 

RP 

RC 

EXC.PYR. 

FUEL 

AIR     IN 

S 

GAS.C 

C 

KG/S 

KG/S 

KG/S 

PCT 

KG/S 

i 

0. 

510. 

296. 

26. 

25. 

.240 

.168 

.072 

92 

8 

2,35 

2 

60. 

560. 

422. 

30. 

25. 

.231 

.162 

.069 

85 

9 

2,35 

3 

120. 

572. 

461. 

41. 

25. 

.222 

.156 

•  067 

79.2 

2.36 

4 

180. 

574. 

482. 

58. 

25. 

.214 

.  150 

.064 

72 

8 

2.36 

5 

240. 

571. 

490. 

78. 

26. 

.205 

.  143 

.061 

66 

6 

2.37 

6 

30  0. 

563. 

492. 

98. 

28. 

.1  96 

.137 

.059 

60.8 

2.37 

7 

360. 

553. 

489. 

117. 

30. 

.187 

.131 

.056 

55 

2 

2.38 

a 

420. 

54  0. 

483. 

134. 

33. 

.178 

.125 

.053 

49.3 

2.38 

9 

480. 

526. 

474. 

149. 

37. 

.169 

.119 

.051 

44 

7 

2.39 

10 

540. 

510. 

463. 

162. 

42. 

.161 

.1  12 

.048 

39.9 

2.39 

1  1 

600. 

4  94. 

451. 

174. 

46. 

.152 

.106 

.045 

35 

4 

2.40 

12 

660. 

476. 

437. 

184. 

51. 

.143 

.  100 

.043 

31 

1 

2.40 

13 

720. 

457. 

423. 

192. 

56. 

.134 

.0  94 

.040 

27.1 

2.40 

14 

780. 

438. 

408. 

199. 

61. 

.125 

.0  87 

.037 

23.3 

2,40 

15 

840. 

417. 

392. 

204. 

66. 

.116 

.081 

.035 

19.3 

2.40 

16 

900. 

396. 

375. 

209. 

70. 

.107 

.0  75 

.032 

16 

6 

2.40 

17 

960. 

374. 

358. 

212. 

74. 

.098 

.069 

.029 

13 

7 

2.39 

18 

1020. 

352. 

34  1. 

214. 

78. 

.0  89 

.062 

.027 

11 

0 

2.38 

19 

1080. 

328. 

324. 

216. 

81. 

.080 

.056 

.024 

8.6 

2.37 

20 

1  140. 

304. 

307. 

216. 

84. 

.071 

.049 

.021 

6.5 

2,35 

21 

1200. 

279. 

289. 

216. 

87. 

.061 

.043 

.018 

4 

7 

2.33 

22 

1260. 

253. 

27  1, 

214. 

89. 

.052 

.036 

.016 

3.1 

2.29 

23 

1320. 

226. 

252. 

213. 

91. 

.042 

.030 

.013 

1.9 

2.25 

24 

1380. 

199. 

233. 

210. 

93. 

.033 

.023 

.010 

.9 

2.19 

25 

1440. 

169. 

213. 

207. 

94. 

.022 

.016 

.007 

2 

2,10 

26 

1500. 

134. 

190. 

203. 

95. 

.011 

.008 

.003 

0 

1  .96 

27 

1560. 

100. 

166. 

198. 

96. 

.000 

.000 

.000 

.0 

1,74 

28 

1620. 

92. 

152. 

193. 

96. 

.000 

.000 

.  000 

0 

1,68 

29 

1680. 

87. 

143. 

137. 

96. 

.000 

.000 

.000 

0 

1.63 

30 

1740. 

83. 

135. 

181. 

96. 

.000 

.000 

.000 

0 

1.59 

31 

1800. 

79. 

129. 

175. 

95. 

.000 

.000 

.000 

0 

1.55 

CRIB    FIRE 

TIME 

HEAT     BALANCE 

Q- 

-WALL 

GAS 

FLOW 

WNO     RAD 

WALL     CNV 

WALL     RAO 

Q-FIRE 

SUM 

PCT 

PCT 

PCT 

PCT 

w 

J 

0. 

55.544 

3.303 

19 

249 

21  . 

906 

2.531*06 

6.250+07 

60. 

63 

693 

4.410 

1  1  .122 

20. 

80  0 

2.4  38'f06 

1 

,092+08 

120. 

6  7.4  96 

4.843 

a 

.615 

19. 

040 

2.3454-06 

1.481+08 

ISO. 

70.479 

5.1  13 

7 

081 

17. 

334 

2.252+06 

1 

811+08 

240. 

72.815 

5.238 

6 

.119 

15. 

82  8 

2. 159+06 

2.095+08 

300. 

74 

815 

5.278 

5 

.434 

14. 

47  3 

2. 066+06 

2.342+08 

360. 

76 

584 

5.253 

4.919 

13. 

244 

1.973+06 

2 

557+08 

420. 

78.203 

5.181 

4 

507 

12. 

109 

1. 880+06 

2.745+08 

480. 

79.721 

S.074 

4 

.160 

1  1  . 

046 

1. 786+06 

2 

.908+08 

540. 

81.173 

4.941 

3 

.852 

10. 

036 

1.693+06 

3 

049+08 

600. 

82 

590 

4.787 

3 

.564 

9. 

062 

1. 599+06 

3 

.170+08 

660. 

83 

.993 

4.61  9 

3 

.280 

8. 

112 

1.505+06 

3 

.273+08 

720. 

85 

,405 

4.440 

2 

.987 

7. 

173 

1  .4  11  +  06 

3 

,359+08 

780. 

86.848 

4.254 

2.673 

6. 

232 

1.317+06 

3.429+08 

840. 

88 

347 

4.064 

2 

.324 

5. 

281 

1.222+06 

3 

485  +  08 

900. 

89 

926 

3.872 

1 

.924 

4. 

299 

1 . 127+06 

3 

,527+08 

960. 

91 

618 

3.681 

1.461 

3. 

268 

1. 032+06 

3 

,556+08 

1020. 

93 

446 

3.492 

.919 

2. 

150 

9.368+05 

3 

,573+08 

1080. 

95.478 

3.308 

.318 

. 

902 

8.407+05 

3 

580+08 

1  140. 

96 

899 

3.101 

- 

.220 

-. 

632 

7.440+05 

3.580+08 

1200. 

97.133 

2.867 

-1 

.288 

-2. 

188 

6. 465+05 

3 

580+08 

1260. 

97.353 

2.647 

-3 

.035 

-3. 

814 

5.479+05 

3.580+08 

1320. 

97 

.559 

2.441 

-5 

.736 

-5. 

614 

4.476+05 

3 

580+08 

1380. 

97 

.750 

2.250 

-10 

.04  1 

-7. 

776 

3.446+05 

3.580+08 

1440. 

97.922 

2.078 

-17 

571 

-10. 

734 

2.366+05 

3 

580+08 

1500. 

98.072 

1.928 

-34 

.79  1 

-16. 

210 

1. 134+05 

3.580+08 

1560. 

98 

.139 

1.861 

-73 

.76  0 

-26. 

340 

0.000 

3.580+08 

1620. 

98 

.135 

1.865 

-74 

.639 

-25. 

412 

0.000 

3 

580+08 

1680. 

98 

.127 

1.873 

-75 

.239 

-24. 

810 

0.000 

3. 560+08 

1740. 

98 

.1  16 

1.884 

-75 

.698 

-24. 

354 

0.000 

3.580+08 

1800. 

98 

.104 

1.896 

-76 

.063 

-23. 

994 

0.0  00 

3 

580+08 

PAGE     NO. 

1                RUN 

NO. 

N.P. 

VELOCITY 

MOL.WT 

FUEL 

M/S 

CNTRL 

.41 

1,80 

29,90 

T 

.40 

1,82 

29,86 

T 

.40 

1,83 

29,83 

T 

.40 

1,83 

29,79 

T 

.40 

1.83 

29.76 

T 

.40 

1.83 

29.72 

T 

.41 

1.83 

29.69 

T 

.41 

1,82 

29,65 

T 

.41 

1.82 

29.62 

T 

.41 

1,81 

29.58 

T 

.41 

1,81 

29.55 

T 

.42 

1.80 

29.51 

T 

.42 

1,79 

29.48 

T 

.42 

1.78 

29.44 

T 

.42 

1.76 

29.41 

T 

.43 

1.75 

29.37 

T 

,43 

1.73 

29.33 

T 

,43 

1.71 

29,30 

T 

,44 

1.69 

29,26 

T 

,44 

1.66 

29.23 

T 

,45 

1.63 

29.19 

T 

,45 

1,59 

29.15 

T 

,45 

1,54 

29.12 

T 

,46 

1.48 

29.08 

T 

,47 

1,40 

29.03 

T 

,47 

1.29 

28.98 

T 

.48 

1.13 

23.92 

T 

.48 

1.08 

28.92 

T 

.48 

1.05 

28.92 

T 

.49 

1.02 

28.92 

T 

.49 

.99 

28.92 

T 

PAGE     NO. 

1               RUN 

NO. 

yo2 

»N2 

YC02 

yH20 

VPYR 

PCT 

PCT 

PCT 

PCT 

PCT 

MASS 

MASS 

MASS 

MASS 

MASS 

.129 

•  699 

,  106 

.043 

.028 

.132 

,701 

.  102 

.041 

.027 

.136 

,704 

.09  8 

.040 

.026 

.139 

,706 

.09  4 

.038 

.025 

.143 

.709 

.091 

.037 

.024 

.147 

,71  1 

.087 

.035 

.023 

.150 

,714 

,083 

.034 

.022 

.154 

,716 

.  079 

,032 

.021 

.158 

.719 

.075 

,030 

,020 

.161 

.722 

,072 

.029 

.019 

.165 

.724 

,068 

,027 

,018 

,169 

.727 

,064 

.026 

,017 

,  172 

.729 

,06  0 

.024 

.016 

.176 

,732 

,056 

.023 

.015 

.180 

,735 

,052 

.021 

.014 

.  183 

,737 

,049 

.020 

.013 

.187 

.740 

,045 

.018 

.012 

.191 

.74  2 

.04  1 

.017 

.01  1 

,194 

.74  5 

,037 

.015 

.010 

.  198 

.748 

.033 

.013 

.009 

.202 

.75  0 

.029 

.012 

.008 

.206 

.753 

,025 

.010 

.007 

.210 

.756 

,021 

,009 

.006 

.214 

,759 

,017 

,007 

.004 

.218 

.762 

.012 

.005 

.003 

.224 

.766 

.006 

.003 

.002 

.230 

.770 

,000 

,000 

.000 

.230 

.770 

,000 

.000 

.000 

.230 

,770 

.000 

.000 

.000 

.230 

,770 

.  000 

.000 

.000 

.230 

,770 

.000 

.000 

.000 
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TEST  PROGRAM  1  FOR  WOOD  CRIB  FIRE.  NILSSON'S  FORMULAS 

GEOMETRY  AND  VENTILATION 

WALL  SURFACE  AREA  =       80.0  M2 
FLOOR  AREA  =      20.00  M2 
WINDOW  HEIGHT  =       1.50  M 

AREA  =       4.00  M2 
OPENING  FACTOR  =    4.899  M2.5 
DISCHARGE  COEFF.=   .68 

FUEL  LOAD  PROPERTIES 

FIRE  LOAD  PER  FLOOR  AREA  =   10.0  KG/M2 

FUEL  COMPOSITION 

CARBON    =  44.4  PERCENT  BY  WEIGHT 

HYDROGEN  =   5.4  PERCENT 

OXYGEN    =  38.2  PERCENT 

NITROGEN  =    .0  PERCENT 

WATER      =  12.0  PERCENT 

R  =   5.32 

R0=   1.23 

HEAT  OF  COMBUSTION  OF  DRY  FUEL  =     18.80-I-06  J/KG 

LOWER  ACTUAL  HEAT  OF  COMBUSTION  =    15.07+06  J/KG 

MOLECULAR  WEIGHT  OF  UNBURNT  PYROLYSATES  =   28.97 

CP  OF  PYROLYSIS  GAS  =  (   .1127*TGAS  +   1010.)  J/KG-K 

MAXIMUM  FRACTION  OF  PYROLYSATES  BURNED  =     .70 

GREY-GAS  FLAME  EMISSIVITY  =   .900 

RATE  OF  REGRESSION  =   00.00      M/S 

FUEL  DIMENSION  =   .050  M 

SHAPE  FACTOR  =  2.00 

CRIB  SPACING/HEIGHT  RATIO=   .100 

WALL  THERMAL  PROPERTIES 

THICKNESS  =   .038  M 
DENSITY  =    790.  KG/M3 

THERMAL  CONDUCTIVITY  =      .170  W/M-K 

HEAT  CAPACITY  =      840.  J/KG-K 

EMISSIVITY  =    .50 
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CRIB    FIRE 

PAGE    NO. 

1              RUN 

NO. 

TIME 

TEMP 

KACL    TEMPS 

RP 

RC 

EXC .PYR. 

FUEL 

AIR     IN 

N.P. 

VELOCITY 

MOL. WT 

FUEL 

S 

GAS.C 

C 

KG/S 

K.G/S 

KG/S 

PCT 

KG/S 

M/S 

CNTRL 

1 

0. 

408. 

217. 

25. 

25. 

.176 

.  123 

.053 

94 

.7 

2.36 

.42 

1.74 

29.65 

T 

2 

60. 

452. 

315. 

28. 

25. 

.176 

.  123 

.053 

89 

.4 

2.38 

.42 

1.78 

29.65 

T 

3 

120. 

472. 

357. 

36. 

25. 

.  176 

.123 

.053 

84 

.2 

2.38 

.41 

1.79 

29.65 

T 

4 

180. 

487. 

386. 

49. 

25. 

.176 

.123 

.053 

78 

.9 

2.38 

.41 

1.80 

29.65 

T 

S 

240. 

498. 

408. 

64. 

26. 

.176 

.123 

.053 

73 

.6 

2.38 

.41 

1.80 

29.65 

T 

6 

300. 

506. 

424. 

80. 

27. 

.176 

.123 

.053 

68 

3 

2.38 

.41 

1  .81 

29.65 

T 

7 

360. 

513. 

438. 

96. 

29. 

.176 

.  123 

.053 

63 

0 

2.38 

.41 

1.81 

29.65 

T 

8 

420. 

519. 

449. 

111. 

31  . 

.176 

.  123 

.053 

57 

.8 

2.39 

.41 

1.81 

29.65 

T 

9 

480. 

524. 

458. 

126. 

34. 

.176 

.123 

.053 

52 

.5 

2.39 

.41 

1.32 

29.65 

T 

10 

540. 

524. 

463. 

139. 

38. 

.173 

.  121 

.052 

47 

.3 

2.39 

.41 

1.82 

29.64 

T 

1  1 

600. 

512. 

458. 

152. 

42. 

.165 

.115 

.049 

42 

3 

2.39 

.41 

1.81 

29.60 

T 

12 

660. 

497. 

449. 

164. 

46. 

.156 

.109 

.047 

37 

7 

2.40 

.41 

1.81 

29.56 

T 

13 

720. 

481  . 

438. 

174. 

51  . 

.147 

.103 

.044 

33 

3 

2.40 

.42 

1.80 

29.53 

T 

!♦ 

780. 

464. 

425. 

183. 

55. 

.138 

.0  97 

.041 

29.  1 

2.40 

.42 

1.79 

29.49 

T 

15 

840. 

445. 

411. 

191. 

60. 

.129 

.090 

.039 

25 

2 

2.40 

.42 

1.78 

29.46 

T 

16 

900. 

425. 

396. 

198. 

64. 

.120 

.084 

.036 

21 

6 

2.40 

.42 

1.77 

29.42 

T 

17 

960. 

405. 

381  . 

203. 

68. 

.111 

.078 

.033 

18 

3 

2.40 

.43 

1.76 

29.39 

T 

18 

1020. 

384. 

364. 

207. 

72. 

.102 

.072 

.031 

15 

2 

2.39 

.43 

1  .74 

29.35 

T 

19 

1080. 

362. 

348. 

210. 

76. 

.093 

.065 

.023 

12.4 

2.39 

.43 

1.72 

29.32 

T 

20 

1  140. 

339. 

33  0. 

212. 

79. 

.084 

.059 

.025 

9.9 

2.38 

.44 

1.70 

29.28 

T 

21 

1200. 

315. 

313. 

213. 

82. 

.075 

.053 

.023 

7.6 

2.36 

.44 

1.67 

29.25 

T 

22 

1260. 

291. 

296. 

213. 

85. 

.066 

.046 

.020 

5.6 

2.34 

.44 

1.64 

29.21 

T 

23 

1320. 

266. 

279. 

213. 

87. 

.057 

.040 

.017 

3.9 

2.31 

.45 

1.61 

29.17 

T 

24 

1380. 

240. 

26  1  . 

211. 

90. 

.047 

.033 

.014 

2.5 

2.27 

.45 

1.57 

29.  14 

T 

25 

1440. 

213. 

242. 

209. 

91. 

.038 

.027 

.011 

1.4 

2.22 

.46 

1.51 

29.10 

T 

26 

1500. 

185. 

222. 

206. 

93. 

.028 

.020 

.008 

5 

2.15 

.46 

1.45 

29.06 

T 

27 

1560. 

154. 

20  1. 

203. 

94. 

.017 

.012 

.005 

0 

2.04 

.47 

1.36 

29.01 

T 

28 

1620. 

109. 

175. 

199. 

94. 

.002 

.001 

.001 

0 

1  .80 

.48 

1.17 

28.93 

T 

29 

1680. 

95. 

157. 

194. 

95. 

.000 

.000 

.000 

0 

1.70 

.48 

1.10 

28.92 

T 

30 

1740. 

89. 

146. 

189. 

95. 

.000 

.000 

.000 

.0 

1.65 

.43 

1.06 

28.92 

T 

31 

1800. 

64. 

138. 

183. 

95. 

.000 

.000 

.000 

0 

1.60 

.48 

1.03 

23.92 

T 

32 

1860. 

81  . 

131. 

177. 

95. 

.000 

.000 

.000 

.0 

1  .56 

.49 

1.00 

28.92 

T 

33 

1920. 

77. 

125. 

171. 

94. 

.000 

.000 

.000 

.0 

1  .53 

.49 

.98 

23.92 

T 

CRIB    FIRE 

PAGE    NO. 

1             RUN 

NO. 

TIME 

HEAT    BALANCE 

0- 

-WALL 

V02 

YN2 

YC02 

YH20 

Y?YR 

GAS 

FLOW 

UNO    RAO 

WALU     CNV 

WALL    RAD 

0-FIRE 

SUM 

PCT 

PCT 

PCT 

PCT 

PCT 

PCT 

PCT 

PCT 

PCT 

W 

J 

MASS 

MASS 

MASS 

MASS 

MASS 

a. 

57.683 

2.530 

22 

.491 

17.294 

1.856«^06 

4 

.431+07 

.154 

.717 

.079 

.032 

.021 

60. 

64.928 

3.282 

14.508 

17.280 

l.856<-06 

7.971+07 

.155 

.717 

.079 

.032 

.021 

120. 

68.196 

3.673 

11 

.511 

16.622 

1.856+06 

.1 10+08 

.155 

.717 

.078 

.032 

.021 

180. 

70.573 

3.978 

9 

.568 

15.88  1 

1  .  8564-06 

1 

394+08 

.  1  55 

.717 

.078 

.032 

.021 

240. 

72.319 

4.215 

8 

266 

15.200 

1.856+06 

.655  +  08 

.155 

.717 

.078 

.032 

.021 

300. 

73.7  06 

4.41  1 

7 

306 

14.578 

1 . 856*06 

1 

899+08 

.155 

.717 

.078 

.032 

.021 

360. 

74.838 

4.576 

6 

.571 

14.020 

1.8S6+06 

2 

128+08 

.  155 

.717 

.078 

.032 

.021 

420. 

75.787 

4,71  7 

5.985 

13.515 

1 .856+06 

2 

.345+08 

.155 

.717 

.078 

.032 

.021 

430. 

76.597 

4.341 

5 

.508 

13.059 

1.856+06 

2 

.552+08 

.155 

.717 

.073 

.032 

.021 

540. 

77.663 

4.912 

5.042 

12.382 

1.829+06 

2.743+08 

.  156 

.718 

.077 

.031 

.020 

600. 

79 

542 

4.86  2 

4 

.485 

11.111 

1. 736+06 

2 

906+08 

.  160 

.720 

.073 

.030 

.019 

660. 

81 

169 

4.755 

4 

.077 

10.00  1 

1.643+06 

3 

045+08 

.163 

.723 

.070 

.023 

.018 

720. 

82.701 

4.618 

3.722 

8.962 

1.550  +  06 

3 

.162+08 

.167 

.726 

.066 

.027 

.017 

780. 

84.190 

4.461 

3 

.388 

7.965 

1.456+06 

3 

.262+08 

.171 

.728 

.06  2 

.025 

.016 

840. 

85 

673 

4.291 

3 

.052 

6.989 

1.362+06 

3 

344+08 

.174 

.731 

.058 

.023 

.015 

900. 

87 

180 

4.113 

2 

.697 

6.017 

1.268+06 

3 

410+08 

.178 

.733 

.054 

.022 

.014 

960. 

88.743 

3.930 

2.307 

5.036 

1. 174+06 

3 

462+08 

.181 

.736 

.051 

.020 

.013 

1020. 

90.393 

3.745 

1 

.864 

4.022 

1.080+06 

3 

500+08 

.185 

.738 

.047 

.019 

.012 

1080. 

92.159 

3.560 

1 

.351 

2.951 

9.850+05 

3.525+08 

.189 

.74  1 

.043 

.017 

.01  1 

1  140. 

94.085 

3.378 

.759 

1.783 

8.897+05 

3.539+08 

.193 

.744 

.039 

.016 

.010 

1200. 

96.227 

3.200 

.135 

.453 

7.938+05 

3 

542+08 

.196 

.746 

.035 

.014 

.009 

1260. 

97 

,021 

2.979 

- 

.524 

-1.158 

6.973+05 

3 

542+08 

.200 

.749 

.03  1 

.013 

.008 

1320. 

97 

247 

2.753 

-1 

.857 

-2.747 

5.998+05 

3.542+08 

.204 

.752 

.027 

.011 

.007 

1380. 

97 

458 

2.542 

-3 

97  0 

-4.448 

5. 009+05 

3 

542+08 

.208 

.754 

.023 

.009 

.006 

1440. 

97.655 

2.345 

-7.264 

-6.383 

4.001+05 

3 

542+08 

.212 

.75  7 

.019 

.008 

•  005 

1500. 

97 

836 

2.164 

-12.695 

-8.823 

2.958+05 

3 

542+08 

.216 

.760 

.015 

.006 

.004 

1560. 

97.997 

2.003 

-23 

.094 

-12.532 

1.840+05 

3 

542+08 

.221 

.763 

.010 

.004 

.003 

1620. 

98.133 

1.867 

-62 

659 

-23.888 

2.191+04 

3.542+08 

.229 

.769 

.00  1 

.001 

.000 

1680. 

98 

.137 

1.863 

-74 

.329 

-25.735 

0.000 

3.542+08 

.230 

.770 

.000 

.000 

.000 

1740. 

98.130 

1.870 

-75 

024 

-25.025 

0.0  00 

3 

542+08 

.230 

.770 

.000 

.000 

.000 

1800. 

98 

.120 

1.880 

-75.537 

-24.512 

0.000 

3. 

542+08 

.230 

.770 

.000 

.000 

.000 

1860. 

98 

.108 

1.892 

-75 

.938 

-24.116 

0.000 

3 

542+08 

.230 

.770 

.000 

.000 

.000 

1920. 

98.095 

1.905 

-76 

262 

-23.798 

0.000 

3.542+08 

.230 

.770 

.000 

.000 

.000 
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TEST    PROGRAM    2    FOR    WOOD    CRIB    FIRE,     NIUSSON'S    FORMULAS 

GEOMETRY  ANO  VENTILATION 

WALL  SURFACE  AREA  =:       80.0  M2 
FLOOR  AREA  =      20.00  M2 
WINDOW  HEIGHT  =       1.50  M 

AREA    =^  4.00    M2 

OPENING    FACTOR    =  4.899    M2 . 5 

DISCHARGE    COEFF.=        .68 

FUEL  LOAD  PROPERTIES 

FIRE  LOAD  PER  FLOOR  AREA  =   20.0  KG/M2 

FUEL  COMPOSITION 

CARBON    =  44.4  PERCENT  BY  WEIGHT 

HYDROGEN  =   5.4  PERCENT 

OXYGEN    =  38.2  PERCENT 

NITROGEN  =    .0  PERCENT 

WATER     =  12.0  PERCENT 

R  =   5.32 

R0=   1.23 

HEAT  OF  COMBUSTION  OF  DRY  FUEL  =     18.80+06  J/KG 

LOWER  ACTUAL  HEAT  OF  COMBUSTION  =     15.07+06  J/KG 

MOLECULAR  WEIGHT  OF  UNBURNT  PYROLYSATES  =   28.97 

CP  OF  PYROLYSIS  GAS  =  (  .1127*TGAS  ♦   1010.)  J/KG-K 

MAXIMUM  FRACTION  OF  PYROLYSATES  BURNED  =     .70 

GREY-GAS    FLAME    EMISSIVITY    =        .900 

RATE  OF  REGRESSION  =   00.00      M/S 

FUEL  DIMENSION  =    .050  M 

SHAPE  FACTOR  =  2.00 

CRIB  SPACING/HEIGHT  RATIO=   .200 

WALL  THERMAL  PROPERTIES — ■ 

THICKNESS  =   .038  M 
DENSITY  =    790.  KG/M3 

THERMAL     CONDUCTIVITY     =  .170     W/M-K 

HEAT    CAPACITY    =  840.     J/KG-K 

EMISSIVITY    =       .50 
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CRIB  FIRE 

PAGE  NO. 

1      RUN 

NO^ 

TIME 

TEMP 

WALL  TEMPS 

RP 

RC 

EXC.PYR. 

FUEL 

AIR  IN 

N.P. 

VELOCITY 

MOL.WT 

FUEu 

S 

GAS.C 

C 

K6/S 

KG/S 

KG/S 

PCT 

KG/S 

M/S 

CNTRL 

J 

0. 

739. 

514. 

26. 

25. 

.479 

.2  88 

.  191 

92 

.8 

2.19 

.37 

.83 

30.60 

F 

2 

60. 

858. 

745. 

34. 

25. 

.461 

.2  86 

.176 

35 

.9 

2.17 

.36 

.85 

30.60 

F 

3 

120. 

902. 

815. 

54. 

2S. 

.444 

.236 

.  158 

79 

.2 

2.17 

.36 

.86 

30.60 

F 

4 

160. 

936. 

366. 

35. 

26. 

.426 

.286 

.141 

72 

.8 

2.17 

.36 

.87 

30.59 

F 

5 

240. 

960. 

900. 

122. 

27. 

.409 

.286 

.123 

66 

.7 

2.17 

.36 

.87 

30.59 

T 

6 

300. 

955. 

90S. 

159. 

30. 

.391 

.274 

.117 

60 

.8 

2.18 

.36 

.88 

30.52 

T 

7 

360. 

942. 

397. 

194. 

34. 

.374 

.261 

.112 

55 

.2 

2.20 

.36 

.88 

30.45 

T 

8 

420. 

924. 

884. 

227. 

40. 

.355 

.249 

.107 

49 

.9 

2.21 

.37 

.88 

30.38 

T 

9 

480. 

903. 

366. 

255. 

47. 

.338 

.237 

.  101 

44 

.8 

2.23 

.37 

.83 

30.31 

T 

to 

540. 

879. 

846. 

281  . 

55. 

.321 

.224 

.096 

40 

.0 

2.25 

.37 

.38 

30.24 

T 

11 

600. 

853. 

82  2. 

302. 

63. 

.303 

.212 

.091 

35 

.5 

2.26 

.37 

.88 

30.17 

T 

12 

660. 

824. 

797. 

32  1. 

72. 

.285 

.200 

.086 

31 

.2 

2.23 

.38 

.83 

30.09 

T 

13 

720. 

794. 

769. 

336. 

80. 

.267 

.187 

.080 

27 

.2 

2.30 

.38 

.33 

30.02 

T 

1* 

780. 

762. 

740. 

349. 

89. 

.250 

.175 

.075 

23 

.4 

2.31 

.33 

.38 

29.95 

T 

IS 

640. 

728. 

709. 

360. 

97. 

.232 

.  162 

.070 

20 

.0 

2.33 

.39 

.87 

29.87 

T 

16 

900. 

693. 

678. 

368. 

104. 

.214 

.  150 

.064 

16 

.8 

2.35 

.39 

.87 

29.80 

T 

17 

960. 

656. 

644. 

374. 

111. 

.196 

.137 

.059 

13 

.8 

2.37 

.40 

.86 

29.73 

T 

IS 

1020. 

617. 

610. 

378. 

117. 

.178 

.125 

.053 

11 

.1 

2.38 

.40 

.85 

29.65 

T 

19 

1080. 

576. 

S75. 

380. 

122. 

.160 

.1  12 

.048 

8 

.7 

2.39 

.40 

.84 

29.58 

T 

20 

1  140. 

533. 

540. 

380. 

127. 

.142 

.0  99 

.042 

6 

.6 

2.41 

.41 

.83 

29.51 

T 

21 

1200. 

489. 

504. 

379. 

131  . 

.123 

.086 

.037 

4 

.8 

2.41 

.42 

.81 

29.43 

T 

22 

1260. 

442. 

467. 

376. 

134. 

.105 

.073 

.031 

3 

.2 

2.42 

.42 

.78 

29.36 

T 

23 

1320. 

393. 

429. 

372. 

137. 

.036 

.060 

.026 

1 

.9 

2.41 

.43 

.75 

29.28 

T 

2« 

1380. 

342. 

390. 

366. 

139. 

.066 

.045 

.020 

.9 

2.39 

.44 

.71 

29.21 

T 

25 

144  0. 

286. 

350. 

359. 

141  . 

.0  46 

.032 

.014 

.2 

2.35 

.44 

•  64 

29.12 

T 

26 

1500. 

223. 

307. 

350. 

142. 

.023 

.016 

.0  07 

.0 

2.25 

.46 

.54 

29.03 

T 

27 

1560. 

155. 

261. 

341. 

143. 

.000 

.000 

.000 

.0 

2.06 

.47 

.37 

28.92 

T 

28 

1620. 

140. 

236. 

330. 

143. 

.000 

.000 

.000 

.0 

1  .99 

.47 

.31 

28.92 

T 

29 

1680. 

129. 

219. 

318. 

143. 

.000 

.000 

.000 

.0 

1.94 

.47 

.27 

28.92 

T 

30 

1740. 

121. 

20S. 

306. 

142. 

.000 

.000 

.000 

.0 

1  .69 

.48 

.24 

28.92 

T 

31 

1800. 

lis. 

194. 

294. 

141  . 

.000 

.000 

.000 

.0 

1.85 

.48 

.21 

28.92 

T 

32 

186  0. 

109. 

184. 

282. 

139. 

.000 

.000 

.000 

.0 

1.31 

.48 

.18 

28.92 

T 

33 

1920. 

105. 

175. 

270. 

137. 

.000 

.000 

.000 

.0 

1.78 

.43 

.15 

28.92 

T 

34 

1980. 

101  . 

168. 

260. 

135. 

.000 

.000 

.000 

.0 

1.75 

.48 

.13 

28.92 

T 

35 

2040. 

97. 

161. 

250. 

133. 

.000 

.000 

.000 

.0 

1.72 

.48 

.11 

23.92 

T 

36 

2100. 

94. 

155. 

240. 

130. 

.000 

.000 

.000 

.0 

1.69 

.48 

.09 

28.92 

T 

37 

2160. 

91. 

149. 

231. 

128. 

.000 

.000 

.000 

.0 

1.66 

.48 

.07 

28.92 

T 

38 

2220. 

88. 

144. 

223. 

125. 

.000 

.000 

.000 

.0 

1.64 

.43 

.05 

28.92 

T 

39 

2280. 

85. 

139. 

214. 

122. 

.000 

.000 

.000 

.0 

1.61 

.48 

•  04 

28.92 

T 

40 

2340. 

83. 

135. 

207. 

119. 

.000 

.000 

.000 

.0 

1.59 

.49 

•  02 

28.92 

T 

41 

2400. 

81  . 

131. 

199. 

117. 

.000 

.000 

.000 

.0 

1  .56 

.49 

.00 

28.92 

T 

42 

2460. 

79. 

127. 

192. 

114. 

.000 

.000 

.000 

.0 

1.54 

.49 

.99 

28.92 

T 

CRIB  FIRE 

PAGE  NO. 

1      RUN 

NO. 

TIME 

HEAT  BALANCE 

Q 

-WALL 

Y02 

YN2 

YC02 

YH20 

YPYR 

GAS 

FLOII 

UNO  RAD 

*ALL  CNV 

WALL  RAO 

Q-FIRE 

SUM 

PCT 

PCT 

PCT 

PCT 

PCT 

PCT 

PCT 

PCT 

PCT 

w 

J 

MASS 

MASS 

MASS 

MASS 

MASS 

0. 

SI 

239 

5 

.446 

1  1 

.986 

31  .325 

4.342+06 

1 

.128+08 

.055 

.632 

.  176 

.0  74 

.071 

60. 

59.890 

a 

.578 

4 

.838 

26.713 

4.305+06 

1 

.943+08 

.056 

.635 

.  177 

•  074 

.067 

120. 

62.785 

10 

.011 

3 

.335 

23.867 

4.301  +  06 

2 

.647+03 

.056 

.639 

.  178 

.0  73 

.061 

lao. 

64 

837 

11 

.207 

2 

.533 

21 .406 

4.303+06 

3 

.255+08 

.057 

.644 

.  179 

.073 

.054 

240. 

66 

133 

12 

.113 

2 

.058 

19.70  6 

4. 310+06 

3 

.828+08 

.057 

.648 

.  180 

.073 

.047 

300. 

68.439 

12 

.473 

1 

.715 

1  7  .  38  3 

4. 125+06 

4 

.300+08 

.064 

.653 

.173 

.070 

.046 

360. 

70.315 

12 

.500 

1 

.524 

15.674 

3.939+06 

4 

,707+08 

.071 

.658 

.  166 

.067 

.044 

420. 

72.040 

12 

.361 

1 

.389 

14.222 

3.754+06 

5 

,058+08 

.078 

.663 

.  158 

.064 

.042 

480. 

73.682 

12 

.102 

1 

.287 

12.939 

3.568+06 

5 

,363+08 

.036 

.669 

.  ISO 

.061 

.040 

540. 

75 

291 

1  1 

.7S4 

1 

.203 

1  1.759 

3.331+06 

5 

,626+03 

.093 

.674 

.  142 

.058 

.037 

600. 

76.901 

11 

.339 

1 

.123 

10.638 

3. 195+06 

5 

,351+08 

.101 

.679 

.  135 

.054 

.035 

660. 

78 

537 

10 

.370 

1 

.056 

9.543 

3.0  03+06 

6 

043+08 

.  103 

.684 

.  127 

.051 

.033 

720. 

80. 

222 

10 

.360 

.978 

8.444 

2.821+05 

6 

202+08 

.1  16 

.690 

.  119 

.048 

.031 

780. 

81 

9  79 

9 

.318 

.888 

7.317 

2.633+06 

6.332+08 

.123 

.695 

.111 

.045 

.029 

840. 

83 

833 

9 

.251 

.778 

6.136 

2.445+06 

6.433+08 

.131 

.700 

•  103 

.042 

.027 

900. 

85 

814 

8 

.668 

.64  0 

4.374 

2. 256+06 

6 

508+03 

.1  39 

.706 

.09  5 

.033 

.025 

960. 

87. 

96  0 

a 

.074 

.464 

3.49  6 

2.067+06 

6 

557+08 

.  146 

.711 

.037 

.035 

.023 

1020. 

90. 

318 

7 

.476 

.243 

1.956 

1.877+06 

6 

582+08 

.  154 

.716 

.  079 

.032 

.021 

1080. 

92. 

946 

6 

.878 

.010 

.165 

1.686+06 

6 

583+08 

.162 

.722 

.071 

.029 

.019 

1140. 

93. 

846 

6 

.154 

- 

.314 

-1.899 

1.494+06 

6 

583+08 

.169 

.727 

.063 

.026 

.017 

1200. 

94, 

558 

5 

.442 

-1 

.020 

-4.049 

1. 300+06 

6 

583+08 

.177 

.733 

.055 

.022 

.015 

1260. 

95. 

238 

4 

.762 

-2 

.239 

-6.376 

1. 103+06 

6, 

533+08 

.185 

.738 

.047 

.019 

.012 

1320. 

95. 

883 

4 

.117 

-4 

.264 

-8.925 

9.040+05 

6 

583+08 

.192 

.744 

.039 

.016 

.010 

1380. 

96. 

485 

3 

.515 

-7 

.700 

-1  1.364 

6.995+05 

6.583+08 

.200 

.749 

.031 

.012 

.008 

1440. 

97. 

049 

2 

.951 

-14 

.093 

-15.683 

4.853+05 

6 

583+03 

.209 

•  755 

.022 

.009 

.006 

1500. 

97.574 

2 

.426 

-28 

.573 

-21 .936 

2. 439+05 

6 

583+03 

.219 

.762 

.012 

.005 

.003 

1560. 

97. 

982 

2 

.018 

-66 

.264 

-3  3.90  9 

0.000 

6 

583+08 

.230 

.770 

.000 

.000 

.000 

1620. 

98. 

049 

.951 

-68 

.413 

-31 .701 

0.000 

6 

583+03 

.230 

.770 

.000 

.000 

.000 

1680. 

93 

085 

.915 

-69 

.814 

-30.274 

0.000 

6. 

583+08 

.230 

.770 

.000 

.000 

.000 

1740. 

98. 

107 

.893 

-70 

.393 

-29.183 

0.000 

6.583+08 

.230 

.770 

.000 

.000 

.000 

1800. 

98. 

122 

.878 

-71 

.769 

-28.312 

0.000 

6 

583+08 

.230 

.770 

.000 

.000 

.000 

1860. 

98 

131 

.869 

-72 

.495 

-27.590 

0.000 

6.533+03 

.230 

.770 

.000 

.000 

.000 

1920. 

98. 

136 

.864 

-73 

.116 

-26.981 

0.000 

6 

583+08 

.230 

.770 

.000 

.000 

.000 

1980. 

98. 

138 

.362 

-73 

.665 

-26.462 

0.000 

6 

583+03 

.230 

.770 

.000 

.000 

.000 

2040. 

98. 

138 

.362 

-73 

.953 

-25.954 

0.000 

6 

583+03 

.230 

.770 

.000 

•  000 

.000 

2  100. 

98. 

136 

.864 

-74 

.573 

-25.606 

0.000 

6 

583+08 

.230 

.770 

.000 

.000 

.000 

2160. 

98. 

133 

.867 

-74 

.763 

-25.201 

0.000 

6 

583+08 

.230 

.770 

.000 

.000 

.000 

2220. 

98 

128 

.372 

-75 

.  1  14 

-24.888 

0.000 

6. 

583+08 

.230 

.770 

.000 

.000 

.000 

2280. 

98. 

123 

.877 

-75 

.413 

-24.604 

0.000 

6. 

583+08 

.230 

.770 

.00  0 

.000 

.000 

2340. 

98. 

1  16 

.884 

-75 

.679 

-24.347 

0.  000 

6 

583+08 

.230 

.770 

.000 

.000 

.000 

2400. 

98. 

109 

.891 

-75 

.769 

-24.079 

0.000 

6, 

583+08 

.230 

.770 

.000 

.000 

.000 

24-0. 

98.100 

.900 

-75 

.979 

-23.867 

0.000 

6, 

583+08 

.230 

.770 

.000 

•  000 

.000 
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TEST  PROGRAM  FOR  PVTFIX  ROUTINE.  VARIABLE  WALL  PROPERTIES 

GEOMETRY  AND  VENTILATION 

WALL  SURFACE  AREA  =       80.0  M2 
FLOOR  AREA  =      20.00  M2 
WINDOW  HEIGHT  =       1.50  M 

AREA    =  4.00     M2 

OPENING     FACTOR    =  4.899     M2.5 

DISCHARGE    COEFF.=        .68 

FUEL    LOAD     PROPERTIES 

FIRE    LOAD    PER     FLOOR     AREA    =        20.0    KG/M2 

FUEL    COMPOSITION 

CARBON  =    44.4    PERCENT     BY    WEIGHT 

HYDROGEN    =        5.4    PERCENT 

OXYGEN  =    38.2    PERCENT 

NITROGEN    =  .0    PERCENT 

WATER  =     12.0    PERCENT 

R    =       5.32 

R0=        1.23 

HEAT  OF  COMBUSTION  OF  DRY  FUEL  =     18.80+06  J/KG 

LOWER  ACTUAL  HEAT  OF  COMBUSTION  =     15,07+06  J/KG 

MOLECULAR  WEIGHT  OF  UNBURNT  PYROLYSATES  =   28.97 

CP  OF  PYROLYSIS  GAS  =  <  .1127*TGAS  +   1010.)  J/KG-K 

MAXIMUM  FRACTION  OF  PYROLYSATES  BURNED  =     .70 

GREY-GAS     FLAME    EMISSIVITY    =        .900 

WALL    THERMAL    PROPERTIES 

THICKNESS    =        .038    M 
DENSITY    =  790.     KG/M3 

THERMAL    CONDUCTIVITY    ARRAY     (W/M-K) 
TEMPERATURE  CONDUCTIVITY 

273.0  .210 

372.0  .210 

373.0  .160 

1073.0  .260 


HEAT    CAPACITY  ARRAY     (J/KG-K) 

TEMPERATURE  HEAT    CAPACITY 

273.0  1090. 

372.0  1090. 

373.0  47300. 

383.0  47300. 

384.0  5000. 

413.0  5000. 

414.0  840. 

1073.0  840. 


EMISSIVITY    =        .50 
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VENTILATION  SPECIFIED.  FUEL  PYROLYSIS  ADJUSTED  FOR  WORST  CONDITIONS 


PAGE     NO. 


TIME 

TEMP 

S 

GAS.C 

1 

0. 

720. 

2 

60. 

853. 

3 

120. 

867. 

4 

180. 

913. 

5 

240. 

933. 

6 

300. 

947. 

7 

360. 

959. 

8 

420. 

969. 

9 

480. 

978. 

10 

540. 

986. 

11 

600. 

992. 

12 

660. 

998. 

13 

720. 

1002. 

14 

780. 

1007. 

15 

840. 

1010. 

16 

900. 

1013. 

17 

960. 

1015. 

18 

1020. 

426. 

19 

1080. 

296. 

20 

1  140. 

239. 

21 

1200. 

203. 

22 

1260. 

178. 

23 

1320. 

160. 

24 

1380. 

145. 

25 

1440. 

134. 

26 

1500. 

124. 

27 

1560. 

1  15. 

28 

1620. 

108. 

29 

1680. 

102. 

30 

1740. 

96. 

31 

1800. 

92. 

32 

1860. 

87. 

33 

1920. 

84. 

34 

1980. 

80. 

35 

2040. 

77. 

WALL  TEMPS 

RP 

RC 

EXC.PYR. 

FUEL 

AIR  IN 

C 

KG/S 

KG/S 

KG/S 

PCT 

KG/S 

453. 

26. 

25. 

.419 

.293 

.  126 

93 

.7 

2.23 

729. 

32. 

25. 

.414 

.290 

.124 

87 

.5 

2.20 

787. 

48. 

25. 

.412 

•  289 

.  124 

81 

.3 

2.19 

329. 

75. 

26. 

.411 

.2  38 

.123 

75 

.2 

2.19 

860. 

96. 

27. 

.410 

.287 

.  123 

69 

.  0 

2.18 

881  . 

131. 

29. 

.409 

.287 

.  123 

62 

.9 

2.18 

898. 

146. 

32. 

.409 

.2  86 

.  123 

56 

.7 

2.17 

914. 

214. 

36. 

.408 

.286 

.122 

50 

.6 

2.17 

927. 

252. 

41. 

.408 

.285 

.  122 

44 

.5 

2.17 

938. 

278. 

47. 

.407 

.285 

.  122 

38 

.  4 

2.17 

947. 

312. 

S3. 

.407 

.285 

.122 

32 

.3 

2.  17 

955. 

343. 

57. 

.407 

.285 

.122 

26 

.2 

2.  16 

961. 

363. 

61  . 

.407 

.285 

.122 

20 

.  1 

2.16 

967. 

379. 

64. 

.406 

.284 

.  122 

14 

.0 

2.16 

972. 

392. 

67. 

.406 

.284 

.122 

7 

.9 

2.16 

976. 

403. 

69. 

.406 

.284 

.  122 

1 

.8 

2.16 

980. 

416. 

72. 

.406 

.284 

.  122 

.0 

2.16 

606. 

428. 

73. 

.000 

.000 

.000 

.0 

2.47 

459. 

433. 

74. 

.000 

.000 

.000 

.0 

2.39 

385. 

419. 

76. 

.000 

.000 

.000 

.0 

2.30 

335. 

392. 

77. 

.000 

.000 

.000 

.  0 

2.22 

298. 

364. 

78. 

.000 

.000 

.000 

.0 

2.15 

269. 

337. 

79. 

.000 

.000 

.000 

.0 

2.08 

245. 

312. 

80. 

.000 

.000 

.000 

.0 

2.02 

226. 

290. 

81. 

.000 

.000 

.000 

.0 

1.96 

209. 

271. 

81. 

.000 

.000 

.000 

.0 

1.91 

194. 

255. 

82. 

.000 

.000 

.000 

.0 

1.85 

181. 

240. 

82. 

.000 

.000 

.000 

.0 

1.80 

170. 

227. 

82. 

.000 

.000 

.000 

.0 

1.76 

160. 

216. 

82. 

.000 

.000 

.000 

.0 

1.71 

152. 

205. 

83. 

.000 

.0  00 

.000 

.0 

1.67 

144. 

196. 

83. 

.000 

.000 

.000 

.0 

1  .63 

137. 

187. 

83. 

.000 

.0  00 

.000 

.0 

1.60 

131. 

179. 

83. 

.000 

.000 

.000 

.0 

1.56 

125. 

172. 

83. 

.000 

.000 

.000 

.0 

1.53 

VENTILATION     SPECIFIED,     FUEL    PYROLYSIS    ADJUSTED    FOR    WORST    CONDITIONS 


60. 

120. 

180. 

240. 

300. 

360. 

420. 

480. 

540. 

600. 

660. 

720. 

780. 

840. 

900. 

960. 
1020. 
1080. 
1  140. 
1200. 
1260. 
1320. 
1380. 
1440. 
1  500. 
1560. 
1620. 
1680. 
1740. 
1800. 
1860. 
1920. 
1980. 

040. 


GAS    FLOW 
PCT 

48.339 
58.154 
60.694 
62.608 
64.103 
65.155 
66.052 
66.825 
67.492 
68.080 
68.577 
68.996 
69.354 
69.673 
69.940 
70.160 
70.342 
95.334 
96.914 
9  7.433 
97.701 
97.857 
97.958 
98.025 
98.071 
96.101 
98.121 
98.133 
98.138 
98. 136 
98.135 
98.126 
98.1 19 
98.108 
98.095 


HEAT     BALANCE 
WNO    RAO  WALL    CNV 

PCT  PCT 


4.946 

8.318 

9.415 

10.308 

11.047 

1 1.590 

12.068 

12.493 

12.868 

13.204 

13.494 

13.741 

13.956 

14. 149 

14.313 

14.448 

14.561 

4.666 

3.086 

2.567 

2.299 

2.143 

2.042 

1.975 

1.929 

1  .899 

1.879 

1  .867 

1.862 

1  .862 

1.865 

1.872 

1.881 

1.892 

1.905 


14.745 

5.410 

4.060 

3.244 

2.695 

2.344 

2.081 

1.862 

1.686 

1  .545 

1.429 

1.335 

1.254 

1.187 

1.132 

1.088 

1  .052 

-34.635 

-47.715 

-54.773 

-59.342 

-62.819 

-65.406 

-67.489 

-69.232 

-70.589 

-71.719 

-72.665 

-73.461 

-74.133 

-74.703 

-75. 188 

-75.603 

-75.960 

-76.269 


WALL    RAD 
PCT 

32.022 

28.134 

25.768 

23.808 

22. 137 

2  0.86  0 

19.804 

18.80  7 

17.929 

17.190 

16.52  9 

15.963 

15.443 

15.000 

14.625 

14.315 

14.056 

-65.520 

-52.475 

-45.402 

-40.481 

-37.109 

-34.51 2 

-32.461 

-30.849 

-29.463 

-28.347 

-27,396 

-26.598 

-25.925 

-25.355 

-24,871 

-24,459 

-24.108 

-23.809 


4.419+06 

4.3634-06 

4. 347+06 

4.334+06 

4.324+06 

4.3  17+06 

4.310+06 

4. 305+06 

4. 301+06 

4. 296+06 

4.293+06 

4. 290+06 

4.288+06 

4. 285+06 

4.284+06 

4.282+06 

4.281+06 

0.000 

0.000 

0.000 

0.  000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.  000 

0.000 

0.000 

0.  000 

0.0  00 

0.000 


0-WALL 

SUM 


1.240+08 
2.1 18+08 
2.896+08 
3.599+08 
4.244+08 
4.845+08 
5,41 1+08 
5.945+08 
6.451+0  8 
6.934+08 
7.396+08 
7.841+08 
6.271+08 
8.667+08 
9.092+08 
9.486+08 
9.676+08 
9.876  +  08 
9.876+08 
9.676+08 
9.876  +  08 
9.876+06 
9.876+08 
9.876+08 
9.876+08 
9.876+08 
9.876+08 
9.876+08 
9.676+06 
9.876+08 
9.876+08 
9.876  +  08 
9.876+08 
9.876  +  08 
9.S76+08 


N.P. 

VELOCITY 

MOL.WT 

FUEL 

M/S 

CNTRL 

,38 

1.84 

30.59 

F 

.37 

1.86 

30.59 

F 

.37 

1.87 

30.59 

F 

,36 

1.67 

30.59 

F 

.36 

1.67 

30.59 

F 

.36 

1.67 

30.59 

F 

.36 

1.87 

30.59 

F 

•  36 

1.88 

30.59 

F 

.36 

1.88 

30.59 

F 

.36 

1.88 

30.59 

F 

.36 

1.88 

30.59 

F 

.36 

1.88 

30.59 

F 

.36 

1.88 

30.59 

F 

.36 

1.88 

30.59 

F 

,36 

1.88 

30.59 

F 

.36 

1.88 

30.59 

F 

.36 

1.88 

30.59 

F 

.43 

1.79 

28.92 

T 

.45 

1.67 

26.92 

T 

.45 

1.58 

28.92 

T 

.46 

1.50 

28.92 

T 

.47 

1  .44 

28,92 

T 

.47 

1.38 

28.92 

T 

.47 

1.33 

28.92 

T 

.47 

1.29 

28,92 

T 

.46 

1.25 

28.92 

T 

.48 

1.21 

28.92 

T 

.46 

1.17 

28.92 

T 

.48 

1  .14 

28.92 

T 

.48 

1.11 

28.92 

T 

.48 

1  .08 

28.92 

T 

.48 

1.05 

28.92 

T 

.46 

1.03 

28.92 

T 

.49 

1.00 

28.92 

T 

.49 

.96 

26.92 

T 

PAGE  NO. 

1      RUN 

NO. 

Y02 

YN2 

YC02 

YH20 

YPYR 

PCT 

PCT 

PCT 

PCT 

PCT 

MASS 

MASS 

MASS 

MASS 

MASS 

.057 

.648 

.180 

.073 

.047 

.057 

.646 

.  180 

.073 

.047 

.057 

.64  8 

.  180 

.073 

.047 

,057 

.648 

.  180 

.073 

.047 

,057 

.648 

.  180 

.073 

.047 

.057 

.648 

.  180 

.073 

,047 

.057 

.648 

,  180 

.073 

.047 

.057 

.648 

.  180 

.073 

.047 

.057 

.648 

.180 

.073 

.047 

.057 

.648 

.  180 

.073 

.047 

.057 

.646 

.  160 

.073 

.047 

.057 

.648 

.  180 

.073 

.047 

.057 

.648 

.  180 

.073 

.047 

.057 

.64  8 

.  180 

.073 

.047 

.057 

.648 

,  180 

.073 

.047 

.057 

.64  8 

.  180 

.0  73 

.047 

.057 

,64  8 

.  180 

.073 

,047 

.230 

.770 

.000 

.000 

,000 

.230 

.770 

.000 

.000 

,000 

.230 

.770 

.000 

.000 

.000 

,230 

.770 

.000 

.000 

.000 

,230 

,770 

,000 

.000 

.000 

,230 

.770 

.000 

.000 

,000 

,230 

.770 

.000 

.000 

,900 

.230 

.770 

,000 

,000 

.000 

.230 

.770 

.000 

.000 

.000 

.230 

.770 

,000 

.000 

.000 

.230 

.770 

.000 

.000 

.000 

.230 

.770 

.000 

.000 

.000 

.230 

,770 

.000 

.000 

.000 

.230 

.770 

.000 

.000 

.000 

.230 

.770 

.000 

.000 

,000 

.230 

.770 

.000 

.000 

,000 

.230 

.7  70 

.000 

.000 

.000 

.230 

.770 

.000 

.000 

.000 
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TEST  PROGRAM  FOR  PFLFIX  ROUTINE.  POOL  OPTION 

-GEOMETRY  AND  VENTILATION 

WALL  SURFACE  AREA  =       80.0  M2 
FLOOR  AREA  =      20.00  M2 
WINDOW  HEIGHT  =       1.50  M 

AREA  =      10.00  M2 
OPENING  FACTOR  =    12.247  M2.5 
DISCHARGE  COEFF.=   .68 


FUEL  LOAD  PROPERTIES 

FIRE  LOAD  PER  FLOOR  AREA  =   20.0  KG/M2 

TOTAL  ENTHALPY  OF  PYROLYSIS=      2.40+06  J/KG 

BOILING  TEMPERATURE=  390.  DEG  C 

FUEL  COMPOSITION 

CARBON     =  85.7  PERCENT  BY  WEIGHT 

HYDROGEN  =  14.3  PERCENT 

OXYGEN    =     .0  PERCENT 

NITROGEN  =    .0  PERCENT 

WATER     =     .0  PERCENT 

R  =  14.78 

R0=   3.43 

HEAT  OF  COMBUSTION  OF  DRY  FUEL  =     46.50+06  J/KG 

LOWER  ACTUAL  HEAT  OF  COMBUSTION  =     43.36+06  J/KG 

MOLECULAR  WEIGHT  OF  UNBURNT  PYROLYSATES  =   28.97 

CP  OF  PYROLYSIS  GAS  =  (  •1127*TGAS  +   1010.)  J/KG-K 

MAXIMUM  FRACTION  OF  PYROLYSATES  BURNED  =     .70 

GREY-GAS  FLAME  EMISSIVITY  =    .900 

FUEL  AREA=        5.00  M2 

WALL  THERMAL  PROPERTIES 

THICKNESS  =    .038  M 
DENSITY  =     790.  KG/M3 

THERMAL    CONDUCTIVITY    =  .170     W/M-K 

HEAT    CAPACITY    =  840.     J/KG-K 

EMISSIVITY    =        .50 
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TEST  PROGRAM  FOR  RPFIX  ROUTINE 

GEOMETRy  AND  VENTILATION 

WALL  SURFACE  AREA  -  80.0  M2 

FLOOR  AREA  =      20.00  M2 
WINDOW  HEIGHT  =       1.50  M 

AREA  =      10.00  M2 
OPENING  FACTOR  =   12.247  M2.5 
DISCHARGE  COEFF.=   .68 

FUEL  LOAD  PROPERTIES 

FIRE  LOAD  PER  FLOOR  AREA  =   20.0  KG/M2 

FUEL  COMPOSITION 

CARBON    =  85.7  PERCENT  BY  WEIGHT 

HYDROGEN  =  14.3  PERCENT 

OXYGEN    =    .0  PERCENT 

NITROGEN  =     .0  PERCENT 

WATER      =     .0  PERCENT 

R  =  14.78 

R0=   3.43 

HEAT  OF  COMBUSTION  OF  DRY  FUEL  =    46.50+06  J/KG 

LOWER  ACTUAL  HEAT  OF  COMBUSTION  =    43.36+06  J/KG 

MOLECULAR  WEIGHT  OF  UNBURNT  PYROLYSATES  =   28.97 

CP  OF  PYROLYSIS  GAS  =  (  .1127*TGAS  +   1010.)  J/KG-K 

MAXIMUM  FRACTION  OF  PYRCLYSATES  BURNED  =     .70 

GREY-GAS  FLAME  EMISSIVITY  =   .900 

RATE  OF  PYROLYSIS   (KG/S) 


TIME 

RP 

0. 

.  120 

120. 

.120 

121  . 

.250 

WALL  THERMAL  PROPERTIES 

THICKNESS  =   .036  M 
DENSITY  =    790.  KG/M3 

THERMAL  CONDUCTIVITY  =      .170  W/M-K 

HEAT  CAPACITY  =      840.  J/KG-K 

EMISSIVITY  =   .50 
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INPUTTED  VALUES  OF  RP  ARE  USEO 


TIME 

TEMP 

■ALU     TEMPS 

RP 

RC 

EXC .PVR. 

FUEL 

AIR     IN 

N.P. 

VELOCITY 

MOL. WT 

FUEL 

S 

GAS.C 

C 

KG/S 

KG/S 

KG/S 

PCT 

KG/S 

H/S 

CNTRL 

1 

0. 

432. 

234. 

26. 

25. 

.120 

.084 

.036 

98.2 

6.11 

.43 

.79 

29,23 

T 

2 

60. 

457. 

323. 

29. 

25. 

.120 

.0  84 

.036 

96. 

.4 

6.13 

.42 

.80 

29.23 

T 

3 

120. 

467. 

355. 

37. 

25. 

.120 

.084 

.036 

94.6 

6.13 

.42 

.81 

29.23 

T 

4 

180. 

805. 

666. 

51  . 

25. 

.250 

.175 

,075 

90. 

.9 

5.98 

.39 

,91 

29.56 

T 

5 

240. 

843. 

756. 

71. 

26. 

.250 

.175 

.075 

87. 

.  1 

5.95 

.39 

.92 

29,56 

T 

6 

300. 

856. 

784. 

97. 

27. 

.250 

.175 

.075 

83. 

.4 

5.94 

.39 

.92 

29.56 

T 

7 

360. 

866. 

804. 

128. 

29. 

.250 

.175 

.075 

79. 

,6 

5.94 

.39 

.92 

29.56 

T 

8 

420. 

872. 

817. 

160. 

32. 

.250 

.175 

.075 

75.9 

5.93    ■ 

.38 

.92 

29.57 

T 

9 

480. 

877. 

827. 

190. 

37. 

.250 

.175 

.  075 

72. 

.1 

5.93 

.38 

.92 

29.57 

T 

10 

54  0, 

881  . 

835- 

218. 

42. 

.250 

.175 

.075 

68.4 

5.93 

.38 

.92 

29.57 

T 

11 

600. 

884. 

84  1. 

244. 

49. 

.250 

.175 

.075 

64.6 

S.92 

.38 

.92 

29.57 

T 

12 

660. 

887. 

846. 

267. 

56. 

.250 

.175 

.075 

60.9 

5.92 

.38 

.92 

29.57 

T 

13 

720. 

889. 

851. 

289. 

64. 

.250 

.175 

.075 

57. 

.  1 

5.92 

.38 

.92 

29.57 

T 

14 

780. 

891. 

855. 

309. 

72. 

.250 

.175 

.075 

53. 

,4 

5.92 

.38 

.92 

29.57 

T 

15 

840. 

893. 

858. 

327. 

80. 

.250 

.175 

.075 

49. 

.6 

5.92 

.38 

.92 

29.57 

T 

16 

900. 

894. 

86  1. 

344. 

88. 

.250 

.175 

.075 

45. 

.9 

5.92 

.38 

.92 

29.57 

T 

17 

960. 

896. 

86  4. 

360. 

95. 

.250 

.175 

.075 

42. 

,  1 

5.92 

.38 

.92 

29.57 

T 

18 

1020. 

897. 

366. 

375. 

102. 

.250 

.  175 

.075 

38.4 

5.91 

.38 

.92 

29.57 

T 

19 

1080. 

898. 

868. 

389. 

109. 

.250 

.  175 

.075 

34. 

.6 

5.91 

.38 

.92 

29.57 

T 

20 

1  140. 

899. 

870. 

401. 

116. 

.250 

,175 

.075 

30. 

,9 

5.91 

.38 

.92 

29.57 

T 

21 

1200. 

900. 

872. 

413. 

122. 

.250 

.175 

.075 

27. 

.  1 

5.91 

.38 

.92 

29.57 

T 

22 

1260. 

901  . 

874. 

425. 

127. 

.250 

.175 

.075 

23. 

,4 

5.91 

.38 

.93 

29.57 

T 

23 

1320. 

902. 

876. 

435. 

133. 

.250 

.175 

.075 

19.6 

5.91 

.38 

.93 

29.57 

T 

2* 

1380, 

903. 

877. 

445. 

138. 

.250 

.175 

.075 

15. 

,9 

5.91 

.38 

.93 

29.57 

T 

25 

1440, 

904. 

878. 

454. 

142. 

.250 

.  175 

.075 

12. 

.  1 

5.91 

.38 

.93 

29.57 

T 

26 

1500. 

904. 

380. 

463. 

147. 

.250 

,175 

.  075 

8. 

.4 

5.91 

.38 

.93 

29.57 

T 

27 

1560, 

905. 

88  1. 

471  . 

151. 

.250 

.175 

.075 

4, 

.6 

5,91 

.38 

.93 

29.57 

T 

28 

1620, 

90S. 

882. 

478. 

154. 

.250 

.  175 

.075 

. 

.9 

5.91 

.38 

.93 

29.57 

T 

29 

1680, 

906. 

883. 

485. 

158. 

.250 

.175 

.075 

. 

,0 

5.91 

.38 

.93 

29.57 

T 

30 

1740. 

206. 

476. 

491  . 

161  . 

.000 

.000 

.000 

,0 

5.58 

.46 

.51 

28.92 

T 

31 

1800. 

148. 

357. 

491. 

164. 

.000 

.000 

.000 

, 

,0 

5.08 

.47 

.34 

28.92 

T 

32 

1860, 

126. 

307. 

481  . 

167. 

.000 

,000 

.000 

.0 

4.80 

.48 

.26 

28.92 

T 

33 

1920, 

1  12. 

27  2. 

464. 

169. 

.000 

.000 

.000 

• 

,0 

4.58 

.48 

.19 

28.92 

T 

INPUTTED 

VALUES    OF    RP     ARE    USED 

PAGE     NQ. 

1              RUN    NO. 

TIME 

HEAT     BALANCE 

Q-WALL 

Y02 

YN2 

YC02 

YH20 

YPYR 

GAS 

FLOW 

WND     RAO 

WALL     CNV 

WALL     RAD 

0-FIRE 

SUM 

PCT 

PCT 

PCT 

PCT 

PCT 

PCT 

PCT 

PCT 

PCT 

w 

J 

MASS 

MASS 

MASS 

MASS 

MASS 

0. 

75, 

.919 

3.719 

1  1 

.044 

9. 

317 

3.642»-06 

4.450+07 

.  179 

.755 

,04  2 

.017 

.006 

60. 

80. 

.960 

4.296 

6 

.596 

8. 

147 

3.642+06 

7.672+07 

.  179 

.755 

.042 

.017 

.006 

120. 

82, 

.877 

4.533 

5 

,171 

7. 

419 

3.642*06 

1  .042  +  08 

.  1  80 

.755 

•  04  2 

.017 

.006 

180. 

72 

.523 

10.016 

3 

,312 

14. 

154 

7.588+06 

1.837+08 

.124 

,739 

.088 

.036 

.012 

240. 

76- 

.015 

1 1.532 

1 

.779 

10. 

659 

7.588+06 

2.404+08 

,  124 

,739 

.  089 

.036 

.012 

300. 

77 

.201 

12.086 

1 

.379 

9. 

318 

7.588+06 

2,891+08 

.124 

,739 

.089 

.036 

.012 

360. 

78 

.063 

12.502 

1 

.  122 

8. 

308 

7.588+06 

3.320+08 

.124 

.739 

.089 

.036 

.012 

420. 

78, 

.650 

12.792 

.963 

7. 

606 

7.588+06 

3.710+08 

.  124 

.739 

.089 

.036 

,012 

480. 

79 

.090 

13.013 

.848 

7. 

049 

7.588+06 

4.070+08 

.124 

.739 

.089 

.036 

.012 

540. 

79 

.441 

13.191 

.762 

6. 

606 

7.588+06 

4.405+08 

.124 

.739 

.089 

.036 

,012 

600. 

79 

.729 

13.339 

.694 

6. 

239 

7.588+06 

4.721+08 

,  123 

.739 

.089 

.036 

.012 

660. 

79 

.970 

13.464 

.639 

5. 

927 

7. 588+06 

5.020+08 

.123 

-739 

.089 

.036 

.012 

720. 

80, 

.176 

13,571 

.594 

5. 

659 

7.588+06 

5.305+08 

.123 

.739 

.089 

.037 

.012 

780. 

80 

.355 

13.665 

.556 

5. 

424 

7.588+06 

5.577+08 

.123 

.73  9 

.089 

.037 

.012 

840. 

80 

.513 

13.748 

.523 

5. 

217 

7.588+06 

5.838+08 

.123 

.739 

.089 

.037 

.012 

900. 

80 

.653 

13.822 

.494 

5. 

031 

7.588+06 

6.090+08 

.123 

.739 

.089 

.037 

,012 

960. 

80, 

.778 

13.389 

,469 

4. 

86  4 

7.588+06 

6.332+08 

,123 

.739 

.089 

.037 

.012 

1020. 

80, 

.892 

13,949 

,447 

4, 

71  2 

7.588+06 

6.567+08 

,123 

.739 

.089 

.037 

.012 

1080. 

SO 

.996 

14.005 

.427 

4. 

573 

7.588+06 

6.795+08 

.  123 

.739 

.089 

.037 

.012 

1140. 

81, 

.090 

14.055 

.409 

4. 

445 

7.588+06 

7.016+08 

.123 

.739 

-089 

.037 

.012 

1200. 

81 

.178 

14. 102 

.39  3 

4. 

32  8 

7.588+06 

7.231+08 

.123 

.739 

.089 

.037 

.012 

1260. 

81, 

.258 

14. 146 

.378 

4. 

219 

7.588+06 

7.440+08 

,123 

.739 

.089 

.037 

.012 

1320. 

81 

.333 

14. 186 

.364 

4. 

11  7 

7.588+06 

7.644+08 

.123 

.739 

.089 

.037 

.012 

1380. 

81 

.402 

14.224 

.352 

4. 

023 

7. 588+06 

7.643+08 

,  123 

.739 

.089 

.0J7 

.012 

1440. 

81 

.466 

14.258 

.340 

3. 

935 

7.588+06 

8.038+08 

.1  23 

.739 

.  089 

.037 

.012 

1500. 

81, 

.526 

14.291 

.330 

3. 

853 

7.588+06 

8.228+08 

.123 

.739 

.089 

.037 

.01  2 

1560. 

81 

.582 

14.322 

,320 

3. 

776 

7.588+06 

8.415+08 

.123 

.739 

.089 

.037 

.01  2 

1620. 

81 

.635 

14.350 

.311 

3. 

704 

7. 588+06 

8.598+08 

.123 

.739 

.089 

.037 

.012 

1680. 

81 

.683 

14.377 

.303 

3. 

637 

7.588+06 

8.777+08 

.123 

.739 

.089 

.037 

.012 

1740. 

97 

.678 

2.322 

-55 

.416 

-44. 

732 

0.000 

8.777+08 

,230 

.770 

.000 

.000 

.000 

1800. 

98, 

.014 

1.986 

-64 

.641 

-35. 

282 

0.000 

8.777+08 

,230 

.770 

.000 

.000 

.000 

1860. 

98, 

.094 

1  .906 

-68 

.298 

-31  . 

506 

0.000 

8.777+08 

.230 

.770 

.000 

.000 

.000 

1920. 

98 

.127 

1.873 

-70 

.817 

-29. 

06  6 

0.  000 

8.777+08 

.230 

.770 

.000 

,000 

.000 
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C  COMPF2 MAIN  PROGRAM 

C  CONSTITUTES  REVISION  OF  PROGRAM  COMPF 

C  COMPF2  VERSION  1.0  PROGRAMMED  1  MARCH  1978  BY  V.  BABRAUSKAS 

C  VERSION  1,1    MINOR  REVISIONS   18  AUG  1978 

C 

COMMON    /CNSTS/       AW ALLNt BWOOW, DENS A. G , GASCNT . KTRACE ,MTI ME 
COMMON    /CP/  CPA.CPC0(2) tCPC02(2} tCPH2C2) .CPH20(2) .CPN2(2) t 

1  CP02(2) .CPPYR(2) 

COMMON    /FUEL/  C.CFLPC. CVGROS. CVNET t H, HFLPC. MWPYR . Nt NFLPC. Ot 

1  OFLPC.Rt  ROtREGRES.SH,SHAPE,SIZE,W,WFLPC, WTFUEL 

COMMON    /GP/  AWDOWiBPFtCDtCNV.DTIME,EMS( 2) .HWDOWtIX.IXC» IXLt 

1  J» JM.JPtJPRI NT tK,KD,KH, KIT ER,KNTRL,MWlN,MWOUT«RC»RP .SIGMA 

COMMON    /LOGIC/       FC.FLSPEC . KRI T. NE WPLT . NEWPRP.PLFUEL. PLOT . PNCH» 
1  RPSPEC. VTSPEC 

COMMON    /PLAST/    TBOILC .DHP , STOICH, SI ZE 1 . EITA . EISCAN 
COMMON    /PRBLM/        ADI A. AFLOOR. AW ALL . DENS W. FLO ADf IRUN. OPENF, 
1  PRNTiSTEADY.THICKW 

COMMON    /TEMP/  OENF,  DENU»  T  AMB»TGAS,TI  NPT  «T  1  (  20  ),»T  2(  20  )  .  TSF  ,TSU 

COMMON    /THERML/     CNDA( 2. 1 0 ) . CPW ( 2. 1 0) » DX. EF, EMSA ( 2 . 10 ) . 
1  NCNOtNCPWtNEMStNQGEN«NRP.QGEN( 2t 10) «RPX(  2 f 50 ) 

COMMON    /TITLE/       TITLE(14) 

COMMON    /CPLOT/       BUFX (500) . BUFY( 500 ), SCALX .SCALY . SPECS ( 30) 
INTEGER    TITLE 

LOGICAL    AOIA.EISCAN.FCtFLSPECiKRIT.NEWPLT.NEWPRP.     PLFUEL. 
1  PL OT.PNCH, RPSPEC. STEADY, STOICH, VTSPEC 

REAL    MWIN.MWOUT.MWPYR.MTIME.N.NFLPC 
C 

DATA    ADIA.CD.CFLPC»CNV»CPA»DENSA.EF.G 

1  /.FALSE. (0.6 8. 44. 4 .5.0. 1005..1.18.0.9.9.8/* 

2  HFLPC.IRUN.IX.MTIME.MWPYR.NEWPRP 

3  /5. 4.0. 10.360. .28.97, .TRUE./. 

4  NFLPCOFLPC.  PLOT.  PNCH.REGRES.SH.  SHAPE 

5  /O.O,     38.2,     .FALSE., .FALSE. ,0.0, 0.0. 2./. 

6  SIGMA, SIZEI, STEADY. TINPT.WFLPC 

7  /5.6697E-8.-10.. .FALSE.. 0. 0. 12.0/ 

C  HEAT    CAPACITIES    ARE     GIVEN    IN    THE    FORM       CP( 1) ♦TEMP+CPC 2) 

DATA       CPCO  /0.1185,  1018./. 

1  CPC02  /0.2114.  931./. 

2  CPH20  /0.3E49.  1814./, 

3  CPH2  /0.68e2,  13966./, 

4  CP02  /0.3704.  931./. 

5  CPN2  /0.1127,  1010./. 

6  CPPYR  /0.1127.  1010./ 

C  PROPERTIES  OF  PYROLYSIS  GASES  ARE  ASSUMED  SAME  AS  FOR  NITROGEN 

C 

NA MEL I  ST  /VARS/  ADI A. AFLOOR, AWALL . AWDQW , BPF .CD, CFLPC , CPPYR , C VGROS . 

1  OENSW.DHP.DTIME.EF.EISCAN.  EITA.FLOAD.FLSPEC.HFLPCHWOOW. 

2  I  RUN.  I  X,KTRACE,MT  I  ME,  MWPYR.NEWPLT.NEWPRP.  NFLPCOFLPC,  PLFUEL. 

3  PLOT, PNCH,PRNT,REGRES,RPSP EC, SH, SHAPE, SIZE, STEADY, STOICH, 

4  TINPT.THICKW, TBOILC, VTSPEC, WFLPC,        SCALX  ,  SCALY 
10    READ     ( 1,900, END=150)    TITLE 

900    FORMAT     (13A6,A2) 

WRITE(2,910)  TITLE 
910  FORMAT( IHl, 13A6,A2) 

EISCAN=. FALSE. 

STOICH=. FALSE. 
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STEADY=. FALSE. 

KITER=  0 

TINPT=  0. 

KNTRL=  1 

KTRACE=  0 

NEWPLT=  .FALSE. 

RPSPEC=  .FALSE. 

FLSPEC=  .FALSE. 

VTSPEC=  .FALSE. 

PLFUEL=  .FALSE. 

IRUN=    IRUN+1  lljj 

20     READ     (l.VARS) 

WRITE     (SsVARS) 

IF(ADIA.OR.EISCAN.OR.STOICH)     STEAOY= .TRUE. 

IF     (PNCH)     PUNCH    900.    TITLE 

IF     (NEWPLT. AND. PLOT)     CALL    PLTRST 

IF     (NEWPRP)     CALL    INC 

NEWPRP=  .FALSE. 
30  CALL  ICONOS 

IF  (KNTRL.EQ.2)  GOTO  10 

IF     ( KT RACE. NE. 1 ) GOTO  50 
40  WRITE  CAt901}  IRUN 
901  FORMAT  {•!  RUN  N0.«.I4// 

1  •   TGASl      TGAS2       Fl         F2        DERIVl     K  KO  • 

2  'KH     J    T2(l)      TSF        QFIRE       QFLOW       QRAOW     • 

3  •     RP       RC   •/) 
50  IF  (KITER.EQ.I)  GO  TO  60 

CALL    ECHOIO 
60    IF     (RPSPECJ     GOTO    70 
IF     (FLSPEC)     GOTO    80 
IF     (VTSPECJ     GOTO    90 

IF     CPLFUEL. AND.. NOT. FLSPEC)     GOTO    100 
CALL    CRIB 
GOTO    110 
70    CALL    RPFIX 

GOTO     110 
80    CALL    PFLFIX 

GOT  0    110 
90    CALL    PVTFIX  \ 

GOT  O     110 
100    CALL    POOL 

110    GO    TO    ( 120* 10.130«120)«KNTRL 
C  KNTRL=     1     INITIAL    VALUE 

C  KNTRL=    2    INPUT    ERROR    DETECTED,     PROCEED    TO    NEW    RUN 

C  KNTRL=    3     ITERATION    FAILURE—PRINT    OUT    STEPS    EVEN    IF    KTRACE=    O. 

C  KNTRL=    4    SIMULATION    TIME    LIMIT     EXCEEDED 

120    CONTINUE 
C  INSERT     HERE    ANY    REWIND    COMMAND    TO    BE    DONE    IF    NO    ERROR 

IF     (PLOT-OR.PNCH)     CALL    DOUT 
GO   TO     10 
130     IF     (KTRACE.EQ.l )     GOTO    10 
140    KTRACE=    1 
KITER=     1 

WRITE     (2<.d03)    TGAS 
903    FORMAT     ( 'O— ITERATI  ON    FAILURE— •/•     TGAS=«  .F16  .2  ) 

GO    TO    30 
150    CONTINUE 
C  ENDFILE    5 

STOP 
END 
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SUBROUTINE    CRIB 

CRIB    FIRE    ROUTINE 

EQUATIONS    FOLLOW    NILSSON'S    DATA    FOR    WOOD    CRIBS. 

OTHER    FUEL    CRIBS    CAN    BE    TREATED    IF    PYROLYSIS    CONSTANTS 

ARE    KNOWN. 

COMMON    /CNSTS/       AWALLNt BW DOW, DENSA, G . GASCNT . KTRACE ,MTI ME 

COMMON    /CP/  CPA,CPC0(2),CPC02(2) tCPH2(2) tCPH20(2) fCPN2(2) t 

1  CP02(2).CPPYR(2J 

COMMON    /FUEL/  C t CFLPC tCVGROS.CVNET t H, HFLPC . MWPYR t N. NFLPC, O, 

1  OFLPC.R.R0tREGRES,SH»SHAPE,SIZE, W.WFLPC. WTFUEL 

COMMON    /GP/  AWDOWtBPF,CD,CNV,OTlME,EMS( 2) .HWDOWtlX.lXC. iXLt 

1  Jt JMtJPtJPRI NT* KtKDtKH.K ITER, KNTRL,MWlN«MWOUT*RC«RP« SIGMA 

COMMON    /LOGIC/       FC,FLSPEC,KRIT,NEWPLT,NEWPRP,PLFUELtPLOT.PNCH, 
1  RPSPECVTSPEC 

COMMON    /PLAST/    TBOILCDHP  ,  STOI  GH,  SI  ZE  1  ,EITA  ,EISCAN 

COMMON    /PRBLM/        ADI A, AFLOOR, AWALL , DENSW.FLOAD.I RUN .OPENF , 
1  PRNT, STEADY, THICKW 

COMMON    /OS/  QC0NWtQFlRE«QFLOW*QRADO,QRADW,QWLSUM 

COMMON    /TEMP/  DENF, DENU»T AMB, TGAS.TINPT ,T 1 ( 20 > , T2( 20 ) »TSF ,T SU 

COMMON    /THERML/    CNDAC 2, I  0 ) , CPW( 2, 1 0)  ,DX . EF, EMSA( 2 . 1 0) . 
1  NCND,NCPW,NEMS,NQGEN,NRP.QGEN(2, 10) tRPX(2.50) 

COMMON    /WOUT/  BWORST,FLREM, HRATIO,RMA,RMF,TTIME»VAVGlN. 

1  WA, WB,YC02. YH20, YN2, Y02,YPYR 

LOGICAL    ADIA,EISCAN,FC,FLSPEC,KRIT,NEWPRP,PLFUEL, 
1  PLOT. PNCH,RPSPEC, SCAN, STEADY, STOI CH.VTSP EC 

REAL    MWIN,MW0UT,MWPYR,MTIME,N,NFLPC 

IF     (STEADY)     GOTO     190 

FC=     .FALSE. 

SCAN=     .FALSE. 

QRAOW=0. 

QCONW=    0. 

F2  =  0. 

F1=0. 

DTGAS=10. 

CALL    HEAONG 

START    TIME    LOOP 

DO    170    J=1,JM 

KH=     0 

0ERIV1=     1. 

TGAS2=    0. 

TGAS1=    0. 

TGASP=    2000. 

TGASN=    TAMB 
20    CONTINUE 

K=    0 
30    CONTINUE 

IF     (FLREM.GT.O.)     GOTO    40 

RC=    0. 

RP=    0. 

GO    TO    50 
40    IF     (REGRES.LE.0.0)     GOTO    45 

USE    THIS    FORMULA    IF    INPUT    REGRES    IS    SPECIFIED 

RP=    REGRES*2.*SHAPE/SIZE*FLREM*«( I .- 1 ./SHAPE )*WTFUEL**( 1. /SHAPE) 

GO    TO    50 
45    CONTINUE 
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<;  FUEL  SURFACE  CONTKOL 

C  ASSUME  CRIB  STICK  DENSITY  RHOCR=  500  KG/M**3 

RHOCR=  500. 

REGREN=  l.24E-3/RHOCR*SIZE**-0.6 

RP1=  RE6REN*2.*SHAPE/SIZE*FLREM**( l.-l ./SHAPE )*WTFUEL** ( 1 ./SHAPE) 
C;  CRIB  POROSITY  CONTROL 

RP2=  0.22«WTFUEL/(RHOCR*SIZE)*SH 
C  ROOM  VENTILATION  CONTROL 

RP3=  0.120*AWOOW*SQRT(HWDOW) 

RP=  AMINI  (RPl,RP2.RP3) 
50  RMF=  RMA-I-RP 

YC02=    3.66667*CFLPC*RC/100./RMF 

YH20=     ( WFL  PC*RP+9 • 0*  HFL  PC*RC ) / 1 0  0 . /RMF 

Y02=     (0.23*RMA-R0«RC)/RMF 

YN2=    0.77*RMA/RMF    +NFLPC*RP/1 00 ./RMF 

YPYR=     (RP-RC)/RMF 

IF(YPYR.LT..0>     YPYR=    0. 

MWOUT=    44.«YC02+1 8.*YH20+28.*YN2+32 .*Y02+MWPYR*YPYR 

HRATIO=     l./( l.+( (TGAS/TAMB>*(MWlN/MWOUT)*( I . +RP/RMA)**2} 
1  **0. 3333333333) 

C  NOTE    HIN    IS    TAKEN    AS    POSITIVE 

HIN=    HWOOW*     HRATIO 

ZW= I .-MWOUT*TAMB/MWI N/TGAS 

IF(ZW)195i55.55 
55    VAVGIN=    0.666667*SQRT(2.*G*HIN*ZW) 

RMA=    C0*VAVG1N*H1N*BW00W*DENSA 

RMF=    RMA+RP 

IF     (RMA/R-RP)     60.60.65 
60    RC=    BPF*RMA/R 

GO    TO    70 
65    RC=    BPF*RP 

FC=    .TRUE. 
70    CONTINUE 

QFLOW=    RMF*(YCO2*(TGAS*(0.5*CPCO2(  i  )  *TGAS+CPC02{  2  )  )-TAMB*(0.5* 

1  CPC02( 1 )*TAMB+CPC02(2)) )    +YH20*( TGAS*( 0 . 5*CPH20( 1 ) *T6AS+ 

2  CPH20(2)  )-'TAMB*(0.5*CPH20{  1  )*TAMB+CPH20(2)  ))     +Y02*(TGAS*( 

3  0.5*CP02< 1 )*TGAS+CP02( 2) )-T AMB* ( 0 .5*CP02{ 1 )*T AMB+CP02 ( 2) ) ) 

4  +YN2*(TGAS*(0.5*CPN2(1 )*TG AS+CPN2( 2 ) )-TAMB*( 0. 5*CPN2( 1 )* 

5  TAMB+CPN2(2) ) )     +YPYR* ( T6AS*( 0 .5*CPPYR( 1 ) *TGAS+CPPYR( 2) ) 

6  -TAMB*(0.5*CPPYR(1 )*TAMB+CPPYR( 2 ) ))  ) 
QFIRE=    RC*CVNET 

IF     (ADIA)     GOTO    90 
CALL    DESOLV 

QRADW=  AWALLN*EMS( I )*SI GMA*(T6AS**4.-TSF**4 . ) 

QCONW=    AWALLN*(TGAS-=TSF)*CNV*(  (TGAS-TSF)  *(T  GAS-TSF  ))  **0.  1  666666  7 
90    CONTINUE 

QRADO=    AWDOW*SIGMA*(TGAS**4.-TAMB**4.) 

K=    K  +  l 

F3=F2 

F2  =  Fl 

Fl=     QFIRE-=QFL0W-QRADO-QRADW-QC0NW 

TGAS3=TGAS2 

TGAS2=T6AS1 

TGASi=TGAS 

IF     (F1.LT.0..ANO.TGAS.LT.TGASP)    TGASP=TGAS 

IF     (Fl.GT.O. .AND.TGAS.GT.TGASN)    TGASN=TGAS 

DERIV2=  DERIVI 

IF  (TGAS1.EQ.TGAS2)  GOTO  130 
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DER1VI=(FI-F2)/{TGAS1-TGAS2) 

IF     (KTRACE.GT.O)     WRITE    (4*99)     TGAS 1 »TGAS2.F 1  ,F2 .DERI  VI ,K .KD . 
1  KH. J.T2( 1 ).TSF»QFIRE,QFLOW,QRAOW,RP,RC 

99    F0RMAT(2F9.2.3( I PE9. 2) ,313. 15 , 2( 0PF9.2) . 3 ( IPEIO .3) ,2{0PF7,3)) 

IF     {.NOT. SCAN)     GOTO    95 

IF     (F1/F2.GE.0.0)     GOTO    93 

SCAN=     .FALSE. 

GOTO    100 
93    TGAS=    TGAS-DTGAS 

IF     (TGAS.LT.TAMB)     GOTO    200 

GOTO  120 
95  IF  (DERIVl .LT..0.ANO.ABS(F2) .GT..0001)GOTO  100 

IF(OERIV2.LT-.0.AND.J.GT.2)  GOTO  100 

TGAS=    TGASl+DTGAS/5. 

GOTO     120 
100    DIF=    ABS(F1/QFL0W) 

IF     (OIF.LT.0.002.AND.ABS(TGAS2-TGAS1) .LT.2. )     GOTO     130 

TGAS  =  (F1*TGAS2-F2*TGAS1  )/(Fl-=F2) 

IF     (K.GT.10.ANO.FI.LT.0..AND.TGAS.GT.TGASP)     GOTO    105 

IF     (K.GT.lO.AND.Fl .GT.O.. AND.TGAS.lt  .TGASN)     GOTO     105 

IF     (K.EQ.l.AND.KH.EQ.O)    TGAS=    TGASl    +10. 

IF     (TGAS.GT.2000.)     GOTO     110 

IF     (TGAS.LT.(TAMB+30.))     GOTO    110 

GOTO     120 
105    TGAS=     (TGASN+TGASP)/2. 

GOTO  120 
110  SCAN=  .TRUE. 

TGAS=  1900. 
120  CONTINUE 

IF  (K-200)  30t30»200 
130  CONTINUE 

CALL  RSTA 

FLREM=    FLREM-RP*DTIME 

IFCFLREM.LT.O.)     FLREM=0. 

IF     (QCONW.GT.O. )     QWLSUM=    QWLSUM+( QRADW+QCONW )*DTI ME 

IF     (TTIME     .GE.    MTIME)     GO    TO    210 

IF     (TGAS.LE.353..ANO.J.GE.10)     GO    TO    210 

IF     ( J.EQ.l )     GO    TO     150 

IF     ( JP.LT. JPRINT)     GO   TO    160 

JP=    0 
150    CALL    OUTPUT 
160    JP=    JP+1 

TTIME=    TTIME+DTIME 
170    CONTINUE 

END    TIME    STEP    DO-LOOP 
180    CONTINUE 
185    CALL    OUTPUT 

RETURN 

ERROR  IN  INPUT 
190  CONTINUE 

KNTRL=  2 

WRITE     (2.910) 
910    FORMAT     (///•     CRIB    ROUTINE    DOES    NOT    ACCEPT    STEADY-STATE    CAS£*) 

RETURN 

SQUARE    ROOT     ERROR 
195    CONTINUE 

IF(KTRACE.EQ.l)     WRITE(2.930)    TGAS.RC. RP. YPYR . ZW .RMA.MWOUT 
930    FORMAT     (/•     TGAS= • . F5.0. •     RC=«.E10.4.'     RP=».E10.4. 
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I        •     YPYR=« ,E10.A.«     ZW=«,F6.4.'     RMA=« ,E1 0.4, •     MWOUT=« , F6.1 ) 
C  FAIL    TO    CONVERGE.     ERROR    EXIT 

20  0    CONTINUE 
KNTRL=3 
RETURN 
C  FIRE    IS    OVER    (TRANSIENT    CASE) 

210    CONTINUE 

CALL    OUTPUT 

RETURN 

END 
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SUBROUTINE  DEQNS 

DIFFERENTIAL    EQUATION    SOLVER    BASED    ON       CRANK-NICOLSON    METHOD. 

COMMON    /GP/  AWDOW,BPFtCDiCNV,DTlME.EMS(2}»HWDOWtIX»IXC.IXLf 

1  J. JM.JPt JPRINTtKtKD,KH,KlTERtKNTRLiMWlNtMWOUT.RCtRP»SIGMA 

COMMON    /PRBLM/       AOIAt AFLOOR. AWALL t DENSWtFLO AD»I RUN, OPENF, 
1  PRNTtSTEADY.THlCKW 

COMMON    /TEMP/  OENF » DENU.TAMB.TGAS ,TI NPT tT I ( 20) ,T2( 20 ) ,TSF .TSU 

COMMON    /THERML/    CNDA (2» 1 0 > . CPW( 2, 1 0) , DX, EF, EMSA( 2 . 10) . 
1  NCNO.NCPWtNEMStNQGEN*NRP,QGEN(2* 10) *RPX( 2 « 50 ) 

DIMENSION    A (20) • B( 20) t C(20 ) • CND( 20) , D( 20 ) • HCP(20 ) 

ENTRY    RSTA 

ENTER    HERE    WHEN    READY    FOR    NEW    TIME    STEP     (FINISHED    ITERATING) 

DO     10    I^lflX 
TKI  )=    T2{I  ) 
10    CONTINUE 

TGOLD=    TGAS 

IF     (J.EQ.O)    TGOLD=TAMB 

RETURN 

ENTRY    DESOLV 

SOLVE    DIFFERENTIAL    EQUATION 

KD=     1 

DX1=    DX 

DO     20     1=1. IX 

CND(1)=   TLU(CNDA.NCN0.T1(I)  ) 
20    HCP(I)=    DEN£W«OX/OTlME*TLU(CPW,NCPW«Tl(I)) 

EMS(1)=     l./( 1 ./TLU(EMSA,NEMS.TSF)     +1 ./EF    -1.) 

EMS(2)=    TLU(EMSA,NEMS.TSU) 

DO    50     I=2tlXL 

CNL=     l./(DX/CND(I-l)+DX/CND( I) ) 

IF     (1.EQ.2     )     CN_=     l./(DXl/CND( 1)+0X/CND(2)) 

CNG=     1 ./(DX/CND(I )+DX/CND(I+l) ) 

A(I)=    -CNL 

B(I)=    HCP( I )+CNL+CNG 

C(I )=    -CNG 

0(  I  )=     (HCPd  )-CNL-CNG)*Tl  (I)+CNL*Tl(  I-l  )+CNG*Tl  (  I  +  l  ) 
50    IF     (NQGEN.GT.O)     D(I)=    D( I ) +DX*TLU C QGEN.NQGEN.Tl ( I ) ) 

CNG=     l./(DXl/CND( 1)+DX/CND( 2) ) 

C( 1 )=    -CNG 

CNL=     1./(0X/CND(IXL)+DX/CND(IX) ) 

A(IX)=    -CNL 

ENTER  HERE  WHEN  KD.GT.l  SINCE  PRIOR  EXPRESSIONS  DO  NOT  CHANGE 

30    TSFOLD=   TSF 

ZRF=    TGAS*CTGAS*(TGAS+TSF)+TSF«TSF)+TSF*TSF«TSF 
ZCF=    CNV*( (TGAS-TSF)*(TGAS-TSF))**0. 1666666  667 
HF=     ZCF+SIGMA*EMS( 1)*ZRF 
DENF=     HF»OXl/2./CND( 1) 
ZF=     HF/2./(DENF+l.) 
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B(l)=    HCP( 1 )+ZF+CNG 

D(  1)=     (HCP(  n-ZF-CNG)«Tl(  1  )+ZF*(  TGAS-I-TGOLD) +CNG*T  1  C  2> 
IF     (NQGEN.GT.O)     0(l)=    D( I ) +DX*TLU( QGEN. NQGEN .Tl ( 1 ) ) 
ZRU=    TAMB«{TAMB*(TAMB+TSU)+TSU*TSU)+TSU*TSU*TSU 
ZCU=     1 .87*( (TAMB-TSU)«(TAMB-TSU)>**0. 1666666667 
HU=    ZCU+SIGMA*EMS{2)*ZRU 
DENU=     HU*DX/2./CND(IX} 
ZU=    HU/2«/(DeNU+l«) 
B(IX)=    HCPCIX)+ZU+CNL 

D(IXJ=    {HCP(IX)-ZU-CNL)*T1C1X)+ZU*2.*TAMB+CNL«T1(IXL) 
IF     (NQGEN.GT.O)     D{IX)=    D( I  X >  +  DX*TLU( QGEN . NQGEN.T 1 ( I  X) ) 
CALL    TRIOGF     ( A* B» C. 0tT2 t IX) 
TSF=     (DENF«TGAS+T2< 1))/(DENF+1 .) 
TSU=     (0ENU*TAMB+T2(IX) )/(DENU+l •> 
KD=    KD+1 

IF     {ABS(TSF-TSF0LD).LT.4)     RETURN 
IF     (KD.LE.6)     GO    TO    30 
WRITE     (2»100)    TSF.TSFOLD 
100    FORMAT     CO  FAIL    TO    CONVERGE    D.E.       TSF=«,F7.2.  •  TSFOLO=«     F7.2) 

IF    (XD.LE.aO)     GO    TO    30 
RETURN 
END 


46 


II 


SUBROUTINE    ECHOlD 

SUBROUTINE    TO    ECHO    INPUT    DATA 

COMMON    /CP/  CPA,CPC0<2)tCPC02(2) .CPH2(2) .CPH20(2> tCPN2(2)t 

CP02(2} tCPPYR(2J 
COMMON    /FUEL/  C t CFLPC t C VGROS. CVNET tH .HFLPC . MWPYR.N* NFLPC* O. 

OFLPC.R.RO.REGREStSH, SHAPEtSlZE.W.WFLPC.WTFUEL 
COMMON    /GP/  AWDOW.BPFtCD,CNV,DTlME.EMS(  2).HWD0W,lX.lXCt iXLt 

J. JM, JPt JPRINTtK,KOtKH,KlTER,KNTRLiMWIN»MWOUT«RC.RP«SIGMA 
COMMON    /LOGIC/       FC.FLSPEC t KRI Tt NE WPLT t NEWPRP ,PLFUEL» PLOT » PNCH . 

RPSPEC.VTSPEC 
COMMON    /PLAST/    T BOlLC.OHPt STOlCHt SI ZEl t EIT A t EISCAN 
COMMON    /PRBLM/       AOIA. AFLOOR tA WALL . DENSW tFLO AOt I RUNt OPENFt 

PRNT.STEADYtTHICKW 
COMMON    /THERML/    CNDA( 2t 1 0 ) t CPW( 2t I  0) f DX . EF» EMSA(  2 1 1 0)  i 

NCN0*NCPW«NEMStNQGENtNRP.QGEN(2t 1 0 ) • RPX(  2 « 50 ) 
LOGICAL    AOIA. EI  SCAN. FC.FLSPECt KRl T« NEWPRPiPLFUELt 

PL OT.PNCH.RPSPEC. STEADY, ST OICH.VT SPEC 
REAL    MWIN.  MWOUT.MWPYRtN.NFLPC 
WRITE     (2.90)     IRUN 

WRITE     (2*91)     AWALL*     AFLOORt    HWOOW*     AWOOW.    OPENF.     CO 
WRITE     (2.92J     FLOAD 

IF     (PLFUEL)     WRITE     (2,93)     DHP,TBOILC 

WRITE     (2,94)        CFLPC,HFLPC,OFLPC,NFLPC,WFLPC,R,R0 
WRITE     (2,95)     CVGROS,CVNET»MWPYR,CPPYR.BPF,EF 

IF     (.NOT. (PLFUEL. OR. VTSPEC. OR. RPSPEO)     WRITE     (2,96)    REGRES.SIZE, 
I  SHAPE 

IF     (PLFUEL. AND. .NOT. STOICH)    WRITE    (2,97)     SIZE 

IF     (  .NOT. (PLFUEL. OR. VTSPEC. OR. RPSPEO. AND. SH.GT. 0.0)     WRITE    (2,908) 
1  SH 

IF     (RPSPEC.AND.NRP.EQ.l)     WRITE    (2,913)     RPX(2,1) 

IF     (RPSPEC.AND.NRP.NE.l)     WRITE    (2,91A)     ( ( RPX ( I , J ) , 1= 1 , 2 ) , J= 1 , NRP) 
IF(ADIA)     GOTO    200 
WRITE     (2,98)    THICKW,DENSW 

IF     (NCND.EQ.l)     WRITE    (2,900)     CNDA(2,1) 

IF     (NCND.GT.l)     WRITE    ( 2,90 1 ) ( (CNDA( I , J) . I >1 , 2 ) , J=l , NCNO ) 
IF     (NCPW.EQ.l)     WRITE     (2,904)     CPW(2,l) 

IF     (NCPW.GT.l)     WRITE    ( 2, 905) ( (CPW( I , J ) , I =1 . 2 ) . J-1 , NCPW ) 
IF    (NEMS.EQ.l)     WRITE    (2,902)     EMSA(2,1) 

IF     (NEMS.NE.l)     WRITE     ( 2, 903) (  ( EMSA( I , J ) , 1  =  1  ,2 ) . J=l , NEMS ) 
IF     (NQGEN.EQ.l)     WRITE    (2,906)     QGEN(2,1) 

IF     (NQGEN.GT.l)     WRITE    ( 2, 907) ( ( QGEN( I , J) , 1  =  1 ,2) . J  =  l  .  NQGEN) 
200    IF(AOIA)     WRITE(2,909) 
RETURN 

90  FORMAT     (1H+,     T86,«C0MPF2    VERSION     1.1    -    RUN    N0.»,I4) 

91  FORMAT     ('0 — • GEOMETRY    AND    VENTILATION----*// 

1  •  WALL    SURFACE    AREA    =     • , F8 .  1  , •     M2«/ 

2  •  FLOOR    AREA    =     •♦F8.2.«     M2 • / 

3  •  WINDOW     HEIGHT     =     •,F8.2,«     M«/ 

4  •  AREA    =     •,F8.2.«     M2«/ 

5  •  OPENING  FACTOR  =  •,F7.3, •  M2.5«/ 

6  •  DISCHARGE  COEFF.=  •,F4.2/) 

92  FORMAT  ('0 -FUEL  LOAD  PROPERTIES —*// 

1     •  FIRE  LOAD  PER  FLOOR  AREA  =«,F6.1,«  KG/M2M 

913  FORMAT  ('ORATE  OF  PYROLYSIS  =',F7.2,'  KG/S«) 

914  FORMAT  ( 'ORATE  OF  PYROLYSIS   (KG/S)'/'      TIME        RP' 
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I  50(/3XtF5.0.F9.3) ) 

93    FORMAT     (•    TOTAL    ENTHALPY     OF    PYROLYSIS=     • 

1  .lPE10.2t«     J/KG«/«     BOILING    TEMPERATURE= • , 0PF5 . 0 . •     DEG    C«) 

9A    FORMAT     COFUEL    COMPOSITION'/ 


PERCENT     BY    WEIGHT'/ 
PERCENT'/ 
PERCENT'/ 
PERCENT' / 
PERCENT' / 


1  •     CARBON  =     '.F4.1 

2  •     HYDROGEN    =     ' . FA . I 

3  •     OXYGEN  =     '.F4,i. 

4  '     NITROGEN    =     ' .FA. I 

5  '     WATER  =     ' ,FA.l 

6  •     R    =     '.FS.a/ 

7  •     R0=     '.F5.2) 

95  FORMAT        (     •     HEAT    OF    COMBUSTION    OF    DRY    FUEL    =     '.2PE10.2,'     J/KG'/ 

1  •     LOWER    ACTUAL    HEAT    OF    COMBUSTION    =     ',E10.2,'     J/KG'/ 

2  '     MOLECULAR    WEIGHT    OF    UNBURNT    PYROLYSATES    =     •.0PF6.2/ 

3  '     CP    OF    PYROLYSIS    GAS    =     ( • .F6. A. ' ♦TGAS    +     '.F6.0. 

4  •)     J/KG-K'/ 

5  •     MAXIMUM    FRACTION    OF    PYROLYSATES    BURNED    =    •.F5.2/ 

6  '     GREY-GAS    FLAME    EMISSIVITY    =     •F5.3) 

96  FORMAT     (•     RATE    OF    REGRESSION    =     ',2PE9.2,'     M/S'/ 

1  •     FUEL    DIMENSION    =     '.0PF5.3.'     M'/ 

2  '     SHAPE    FACTOR    =     '•F4.2    /) 

97  FORMAT     (•     FUEL    AREA= • . F 1 0 . 2» '     M2' ) 

908  FORMAT     ('     CRIB    SPACING/HEIGHT    RAT IO= ' . F6.3) 

98  FORMAT     ('0 --WALL    THERMAL    PROPERTIES — ■ — •// 

1  •     THICKNESS    =     '.F5.3,'     M'/ 

2  '     DENSITY    =     '.F6.0»'     KG/M3') 

909  FORMAT     CO — —WALL    THERMAL    PROPERTIES '// 

I  '     ADIABATIC    WALL'//J 

900  FORMAT     C'OTHERMAL    CONDUCTIVITY    =     '.F7,3.'     W/M-K'} 

901  FORMAT     COTHERMAL    CONDUCTIVITY    ARRAY     (W/M-K)'/ 

1  '     TEMPERATURE  CONDUCTIVITY'.     1 0 ( /3X . F7 . 1 .4X. F 1 0.3 > ) 

902  FORMAT     COEMISSIVITY    =     '.FA. 2) 

903  FORMAT     COEMISSIVITY    ARRAY'/'    TEMPERATURE       EMISSIVITY' 
1  10(/3X.F7.1,AX,F10.3)) 

90A    FORMAT     COHEAT    CAPACITY    =     '.F7.0t'     J/KG-K'J 

905  FORMAT        ('OHEAT    CAPACITY    ARRAY     (J/KG-K)'/ 

1  '     TEMPERATURE  HEAT    CAPACITY'.     I 0( /3X, F7. 1 . AX . FIO . 0 ) ) 

906  FORMAT     COWALL    HEAT    GENERATED    =     ',F9.3.'     W/M3') 

907  FORMAT     ('OWALL    HEAT    GENERATED    ARRAY     (W/M3)'/ 

1  '     TEMPERATURE  QGEN'.        10 ( /3X.F7. 1 . 4X.F i 0.3 ) ) 

END 


48 


■k 


SUBROUTINE  ICONDS 

SET  INITIAL  CONDITIONS  AND  CONSTANTS 

COMMON  /CNSTS/   AWALLN, BWDOW, DENSA , G , GASCNT 

COMMON  /CP/       CPA,CPC0(2) .CPC02(2) .CPH2(2 
1    CP02(2l .CPPYR(2) 

COMMON  /FUEL/     C .CFLPC , CVGR0S» CVNET.H.HFLP 
1     OFLPCR.RO.REGRES.SH.SHAPE.SIZE.Wt  WFLPC. 

COMMON  /GP/       AWDOW,BPF,CO,CNV,DTIME,EMS( 
1     Jt JM, JP. JPRINT.K.KD.KH.KITER.KNTRL.MWIN. 

COMMON  /LOGIC/   FC . FLSPEC. KRI T, NE WPLT. NEWPR 
1     RPSPEC. VTSPEC 

COMMON  /PLAST/  TBOILC .DHP , STOICH, SI ZE 1 , EI TA 

COMMON  /PRBLM/   AOI A, AFLOOR. AWALL , DENSW .FLO 
1     PRNT, STEADY. THICKW 

COMMON  /QS/       QC0NW»QFIRE.QFLOW.QRAD0.QRA 

COMMON  /TEMP/     DENF» DENU. TAMB. TGAS. T INPT ,T 

COMMON  /THERML/  CNDA ( 2. 1 0 ) . CPW( 2. 1 0 ) . DX,  EF. 
1     NCN0«NCPW.NEMS«NQ6EN,NRP.QGEN(2. 10) «RPX( 

COMMON  /WOUT/     BWORST. FLREM, HRAT I O. RMA, RMF 
1     VKA,WB.YC02.YH20fYN2.  Y02.YPYR 

LOGICAL  ADIA,EISCAN,FC,FLSPEC.KRIT.NEWPRP,P 
1     PLOT, PNCH, RPSPEC. STEADY. STOICH 

REAL  MWIN.MWOUT.MWPYR.MTI ME, N, NFL PC 

FC=  .FALSE. 

KNTRL=  1 

AWALLN=  AWALL-AWDOW 

BWOOW=  AWDOW/HWDOW 

WTFUEL=  FLOAD*AFLOOR 

OPENF=  AWOOW*SORT(HWOOW) 

MW1N=  28.97 

MWOUT=  MWIN 

TAMB=  298. 

TGAS=  1800. 

IF  (TINPT.GT.O. ) 

IXL=  IX-1 

IXC=  IX/2 

TSF=  TAMB 

TSU=  TAMB 

DENF=  0. 

DENU=  0. 

WA=  6H 

WB=  5H 

IF  < .NOT. FLSPEC) 

WA=  6HWIN00W 

WB=  5HWI0TH 
10  JP=  0 

IF  ( .NOT. STEADY) 

IF  (STEADY)  GOTO 

IF  (DTIME.GT. 0.00001)  GOTO  15 

WRITE  (2.95) 

KNTRL=  2 

RETURN 
95  FORMAT  (///•  FOR  NON-STEADY  PROBLEMS  MUST  S 

1     'GREATER  THAN  ZERO') 
15  IF  (MTIME.GT.DTIME)  GOTO  20 


.KTRACE.MTIME 

) .CPH20(2),CPN2(2) . 

C.MWPYR.N.NFLPC. O. 

WTFUEL 

2) .HWDOW, IX.  IXC.  I  XL. 

MWOUT.RC.RP, SIGMA 

P.  PL FUEL. PLOT, PNCH. 

,EISCAN 
AD.IRUN.OPENF, 

DW.QWLSUM 

1(20). T2(20) .TSF.TSU 

EMSA(2. 10) . 

2.50) 

.TTIME.VAVGIN, 

LFUEL. 


TGAS=  TINPT 


GO  TO  10 


JPRINT=  PRNT/DTIME  + 
20 


(l.O-l.E-6) 


PECIFY  DTIME  '. 
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MRITE  (2.92) 
92  FORMAT  ( /// •  FOR  NON-STEADY  PROBLEMS  MUST  SPECIFY  MTIME  •• 
1     'GREATER  THAN  DTIME* ) 

KNTRL=  2 

RETURN 
20  CONTINUE 

IF  < .NOT, STEADY)  JM=  MTIME/DTIME  +2 

IF  (STEADY)  JM=  2 

DX=  THICKW/IX 

FLREM=  WTFUEL  ^ 

TTIME=  0. 

EMS(1)=  1  ./( 1 ./TLU(EMSA,NEMS,  TAMB)  +1./EF-1.) 

QWLSUM=  0. 

DO  60  1=1, IX 
60  T2(I)=  TAMB 

J=  0 

CALL  RSTA 

TOTAL=  CFLPC+HFLPC+OFLPC+NFLPC+WFLPC 

IF  (TOTAL.LT. 101.1 .AND. TOTAL. GT. 98. 9)  GOTO  70 
CHECK  FOR  ERRORS  IN  FUEL  COMPOSITION 

KNTRL=  2 

WRITE  (2.90) 
90  FORMAT  (///•  SUM  OF  FUEL  COMPOSITION  INPUT  IS  INCORRECT') 

RETURN 
70  C=  CFLPC*( 10./12.) 

H=  HFLPC+10. 

0=  OFLPC*( 10./16.) 

W=  WFLPC*( lO./ia. ) 

N=  NFLPC*( 10./14. ) 

R0=  (C+H/4.-0/2. )*32./1000. 

R=  RO/0.232 

CVNET=  CVGROS*( 1 .-WFLPC/ 1 00 . )-( WFLPC+9. 0*HFLPC) /I 00.*2440 .E+3 
LATENT  HEAT  OF  H20  EVAPORATION=  2440E+3  J/KG  AT  25  C 

RMA=  0. 16*AWDOW*OENSA*SQRT(G*HWDOW) 

RMF=  RMA 

RC=  BPF*RMA/R 

Y02=  0.10 

IF  (STOICH)  EITA=  1. 

IF  (EISCAN)  SIZE=  SIZEl/EITA     ^ 

RETURN 

END 
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SUBROUTINE    INC 

ROUTINE    TO    INPUT    ALL    CONSTANTS    NOT    GIVEN    IN    NAMELIST    -VARS- 

COMMON    /THERML/    CNDA( 2t I  0 ) t CPW( 2 . 10)  .DX ♦ EFt EMSA( 2 . 1 0) . 
1  NCN0,NCPW.NEMStNQGEN,NRP«QGEN(2. 10) .RPX(2t50) 

NQGEN=    0 

READ  (lt91)  NCNtNCPt NEMtNR«NQG 

IF  (NCN.EQ.O)  GO  TO  10 

NCND=  NCN 

READ  (1,90)  (CNDA( l.I),CNDA{2.I ),I=1 ,NCND) 
10  IF  (NCP.EQ.O)  GO  TO  20 

NCPW=  NCP 

READ  (1,90)   (CPW( 1,1  ),CPW(2,I  ),I=1,NCPW) 
20  IF  (NEM.EQ.O)  GO  TO  30 

NEMS=  NEM 

READ  (1,90)  (EMSA(  1,I),EMSA(2,I ),I  =  1,NEMS) 
30  IF  (NR.EQ.O)  GO  TO  40 

NRP=  NR 

READ  (1,90)   (RPX( 1,I),RPX(2,I  ),I  =  1,NRP) 

NQGEN=  NQG 
40  IF  (NQG.EQ.O)  GO  TO  50 

READ  (1,90)  (QGEN( 1,1) ,QGEN(2*I ),I=1 ,NQGEN) 
50  RETURN 

90  FORMAT  (8F10.0) 

91  FORMAT  (1013) 
END 
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SUBROUTINE    OUTPUT 

PRINTS    OUTPUT    DATA 

COMMON    /FUEL/  C t CFLPC. C VGROS iCVNET . Ht HFLPC t MWPYRt N. NFLPC. Ot 

1  OFLPC.R,R0tREGRES»SH,SHAPE.SIZEtW,WFLPC.WTFUEL 

COMMON    /GP/  AWDOW.BPFtCOtCNV,DTlMEtEMS( 2) .HWDOWtlX.IXC. IXL. 

1  Jt JMtJPtJPRI NT .RtKOtKH, KIT ER»KNTRL.MWlN,MWOUT.RC.RP» SIGMA 

COMMON    /LOGIC/       FCt FLSPEC. KRITiNE WPLT» NEWPRP iPLFUEL.PLOT t PNCHt 
1  RPSPEC.VTSPEC 

COMMON    /PLAST/    TBOILC tOHP , STOI CH, SI 2E I .EIT A . EISCAN 

COMMON    /PRBLM/        ADI A, AFLOOR. AWALL . DENSW .FLO AO,I RUN» OPENF » 
I  PRNT.STEADYtTHICKW 

COMMON    /OS/  QCONW.QFlREtQFLOW.QRAOOtQRADWtQWLSUM 

COMMON    /TEMP/  DENF , DENUtTAMB.TGAStTlNPT .T I ( 20) .T2( 20 J »TSF .TSU 

COMMON    /WOUT/  BWORST t FLREMt HRAT 10. RMA.RMF «TTI ME , VAVGI N , 

1  WA,WB.YC02.YH20.YN2.V02.YPYR 

LOGICAL    AOIA.EISCAN.FC.FLSPEC.KRIT.NEWPRP.PLFUEL. 
1  PLOT. PNCH.RPSPEC. STEADY, ST O I CH.VT SPEC 

REAL    MWIN.MWOUT.MWPYR.N.NFLPC 

LOGICAL    DATPRT 

DIMENSION    T2C(3) 

DATA    DATPRT    /.FALSE./ 

IF     (KITER.EQ.I)     RETURN 

IF     (PLOT.OR.PNCH)     CALL     DSTO 


(TTIME.TGAS) 


IF     (ILINE.LE.47)     GO    TO    50 

DATPRT=    .TRUE. 

GO    TO    300 
50    TGASC=    TGAS-273, 

T2C(1)=    TSF-273. 

T2C(2)=    T2(lXC)-273. 

T2C{3)=    TSU-273. 

FUELPC=  FLREM/WTFUEL*100. 

QNORM=  QFIRE/100. 

IF  (QCONW.LT.Oi  QNORM=  ( QFLOW  +  QRADO) /1 00. 

ZFLOW=  QFLOW/QNORM 

ZRADO=  QRADO/QNORM 

ZCONW=  QCONW/QNORM 

ZRADW=  QRADW/QNQRM 

EXCE£S=  RP-RC 

ILINE  =  ILINE+1  \ 

WRITE  (2.90)  J. TTIME»TGASC.T2C.RP,RC. EXCESS. FUELPC.RMA.HRATIO, 
1     VAVGIN, MWOUT.FC 

IF  (FLSPEC)  WRITE  (2.91)  BWORST 

WRITE  (3.92)  TTIME.ZFLOW.ZRADO.ZCONW.ZRADW.QFIRE.QWLSUM. 
1     Y02. YN2. YC02.YH2O.YPYR 

IF  (STOICH)  WRITE  (2.901)  SIZEl 

RETURN 

ENTRY    HEADNG 

START    NEW    PAGE 
IF     (KITER.EQ.I)     RETURN 
IPG    =     1 
300    CONTINUE 
ILINE     =    0 

IF     { .NOT.RPSPEC)     GOTO    315 
WRITE     (2.94)     IPG.IRUN 
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325 


335 


(3«9A)  IPG*IRUN 
400 
(. NOT. FL SPEC)  GOTO  325 
(.NOT.PLFUEL)  WRITE  (2 
(.NOT.PLFUEL)  WRITE  (3 
(PLFUEL)  WRITE  (2,905J 
(PLFUEL)  WRITE  (3,905) 
TO    400 

(.NOT.VTSPEC)     GOTO    335 
WRITE     (2<96)     IPG.IRUN 
WRITE     (3t96)     IPG.IRUN 
GO    TO    400 

IF     (.NOT.PLFUEL)     GO    TO 
WRITE     (2»97)     IPG.IRUN 


WRITE 

GO    TO 

IF 

IF 

IF 

IF 

IF 

GO 

IF 


,95)     IPG.IRUN 
,95)      IPG.IRUN 

IPG.IRUN 

IPG.IRUN 


345 


345 


400 


WRITE 
GO  TO 
WRITE 
WRITE 
IPG  = 
WRITE 
WRITE 


IPG.IRUN 


IPG.IRUN 
IPG.IRUN 

WA,  WB 


90 


(3.97) 

400 

(2,96) 

(3,98) 

lPG+1 

(2,99) 

(3,900) 
IF  (  .NOT.DATPRT)  RETURN 
OATPRT=  .FALSE. 
GO  TO  50 

FORMAT  (IH  ,S4,T6,F8.0,T15,F7.0,T25.3Fe.0.T44.F6.3.T51,F6.3, 
I    T6  0,F6.3,T69,F5.l .T79,F4. 2 , T87 ,F7. 2 ,T96 . F9.2 .T 1 09. F6 .2 . 5X , L I ) 

91  FORMAT     (1H+,     T124,     Fe.2) 

92  FORMAT     (IH     , F7. 0, 4F1 1 .3, I PE 15 . 3, E 1 5. 3, IX , OP ( 5F8.3 ) ) 

94  FORMAT     (1H1,T14,«      INPUTTED    VALUES    OF    RP    ARE     USED* 
1  TIOO.'PAGE    NO. •,I3,Tll5.'RUN    NO.*, 14//) 

95  FORMAT     (1H1,T14,      •     FUEL    PYROLYSIS     (CRIB)     SPECIFIED,     VENTILATION     • 

1  •     ADJUSTED    FOR    WORST    CONDITIONS'.    TIOO.'PAGE    NO. • , 13 .T 1 1 5. 

2  'RUN    NO. ',14//) 

96  FORMAT  (1H1,T14.'  VENTILATION  SPECIFIED,  FUEL  PYROLYSIS  ADJUSTED' 

1  •  FOR  WORST  CONDITIONS' .TIOO.'PAGE  NO. ' . I3.T 1 1 5 , 'RUN  NO. '.14/ 

2  /) 

97  FORMAT  (1HI,T14,'  THERMOPLASTIC  POOL  FI RE • , T  100, ' PAGE  NO.' 
1     I3,T115,'RUN  NO.',  14//) 

98  FORMAT  ( IH  I  ,T  14, 'CRIB  Fl RE ' .T 1 00 . ' PAGE  NO. '.13. 
I    T115,'RUN  NO. ',14/) 

99  FORMAT  ( IH0,T10, 'TIME' .T17. 'TEMP' .T30. 'WALL  TEMPS ' .T47 .' RP' , 

1  T54, 'RC ,T59, 'EXC.PYR. ' .T70. 'FUEL' .T78. 'AIR  I N  '  .T90 . ' N.P . ' « 

2  T98, 'VELOCITY' ,T  1  09. ' MOL . WT '  .Tl 18, 'FUEL'  .T124, A6/ 

3  T11,'S',T17,  'GAS.C  .T34  .  '  C  '  .T46.  •  KG/S  '  .T53.'KG/S'  .T61  . 

4  'KG/S' .T71 ♦ 'PCT' .T79.'KG/S' .TlOl . 'M/S' .T108. 

5  Tl IS.'CNTRL' ,T124,A5/) 

900  FORMAT  (T5. 'TIME' .T24. 'HEAT  B ALANCE • ,T75, ' Q- WALL • , T88 , • Y02' ,T96» 
2     •YN2',T10A, 'YC02' ,T 1 12, ' YH 20 ' ,T 120 , ' YPYR • / 

2  T13,'GAS  FLOW' ,T25, 'WND  RAD'  ,T36,  '  WALL  CNV  ,T47  ,  '  WALL  RAD', 

3  T6  0, 'Q-FIRE',T76,  'SUM'  .T88.'PCT'  ,T96, 'PCT' .T I  0  4, 'PCT' .Tl 12. 

4  'PCT' .T 120. 'PCT'/ 

5  T16, 'PCT' ,T27,  'PCT' .T39.  'PCT' .T49, 'PCT' .T62.' W'  ,T77, 

6  ' J' ,T8a, 'MASS' ,T96. 'MASS' .T104. 'MASS' ,T112.'MASS'. 

7  T120,  'MASS'/) 

901  FORMAT     ( // '     STOICHIOMETRIC    FUEL     SIZE=     '.F8.3.'     M2 ' ) 

905    FORMAT     (1H1,T14,      '     FUEL    PYROLYSIS     (POOL)     SPECIFIED,     VENTILATION' 

1  •     ADJUSTED    FOR    WORST    CONDITIONS',     TIOC'PAGE    NO.  '  ,  13  .T  1  1  5 , 

2  'RUN    NO.  •  ,14//) 

END 
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SUBROUTINE    PFLFlX 


PESSIMIZATION    ROUTINE 

FIXED    FUEL    PYROLYSISt     WINDOW    WIDTH    VARIED    FOR    WORST 

BURNING    CONDITIONS. 

EITHER    CRIB    OR    POOL    FUELS    ACCEPTED 

FOR    POOL    FUELS t     MUST     ALSO    SET    PLFUEL=.TRUE. 

COMMON    /CNSTS/       AW ALLN t BW DOW» DENSA, G . GASCNT t KTRACE ,MTI ME 

COMMON    /CP/  CPA,CPC0(2}tCPCO2(2) tCPH2(2) tCPH20(2) .CPN2(2). 

1  CP02(2) tCPPY«(2) 

COMMON    /FUEL/  C . CFLPC tCVGROS t CVNET t H iHFLPC . MWPYR t Nt NFLPCt O. 

1  OFLPC.R. RO.REGRESiSH.SHAPEfSIZEtW.WFLPC. WTFUEL 

COMMON    /GP/  AWDOW.BPF,CD.CNV,DTIMEtEM£( 2) »HWDOW,IX.IXC.  IXL. 

1  Jt JMtJPtJPRI NT tKtKOiKHt KIT ERtKNTRLtMWiNtM WOUTtRCtRP • SIGMA 

COMMON    /LOGIC/       FC.FLSPECtKRIT t NEWPLT. NEWPRP . PLFUELtPLOT. PNCH t 
I  RPSPEC.VTSPEC 

COMMON    /PLAST/    T BOILCi DHP. STOICH, SI ZEl t EITA , EISCAN 

COMMON    /PRBLM/       ADI A, AFLOOR, A WALL , DENS W. FLOAD . I RUN» OPENF. 
1  PRNT.STEAOY.THICKW 

COMMON    /QS/  QCONW,QFIRE,QFLOW»QRADO,QRADW,QWLSUM 

COMMON    /TEMP/  DENF, DENUt T AMBt TGASt T INPT «T 1 ( 20 ) »T2{ 20 ) .TSF ,TSU 

COMMON    /THERML/    CNDA< 2 1 I  0 ) i CPW( 2. I  0 ) . DX ♦ EF. EMSA( 2 . 10) . 
1  NCNDf NCPWtNEMS«NQGENtNRP*QGEN(2» 10) «RPX(2«50} 

COMMON    /WOUT/  BWORSTt FLREMt HRATI Ot RMA.RMF .TTI ME . V AVGI Nt 

1  WA»WB.YC02fYH20.YN2» Y02tYPYR 

LOGICAL    ADIA.EISCAN.FCtFLSPEC.KRITtNEWPRP.PLFUEL, 
1  PLOT.PNCHtRPSPEC. SCANt  STEADY. STOI CH, VTSPEC 

REAL    MWIN.MWOUT. MWPYR. MTI ME.Ni NFLPC 

IF     (STEADY)     GOTO     190 

FC=     .FALSE. 

SCAN=     .FALSE. 

QRADW=0. 

QCONW=    0. 

F2=0. 

F1=0. 

DTGAS=10. 

TB0IL=TB0ILC+273. 

CALL    HEADNG  , 

START    TIME    LOOP 

DO     170     J=1.JM 

KH=    0 

DERIV1=     1. 

TGAS2=    0. 

TGASl=    0. 

TGASP=    2000. 

TGASN=    TAME 
20    CONTINUE 

K=    0 
30    CONTINUE 

KR=    0 

RMA=    0.666667*CD*0.5*AWDOW*DENSA«SQRT(G*HWDOW*C 1 .-TAMB/TGAS ) ) 

IF     (FLREM)     220,220.35 

AS    SOON    AS    FUEL    IS    EXHAUSTED    PROGRAM    MUST     STOP. 
SINCE    WINDOW    SIZE    WOULD    NOT    BE     WELL    DEFINED. 
35    IF     { .NOT.PLFUELI     GOTO    40 
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RP=    SIZE*EF*SI GMA*(TGAS**4.-TB0IL**4,>/0HP 

PLUME=     S1ZE*0.0014.*CVNET/DHP 

PROP=     1  .->(TGAS**4.-TBOIL**4.}/(1700.**4.=TBOIL«*4.) 

IF     CPROP.LT.O.)     PROP=    0. 

RP=    RP+PROP*PLUME 

GO    TO    50 
40    IF     (REGRES.LE.0.0)     GOTO    45 
C  USE    THIS     FORMULA     IF     INPUT     REGRES     IS    SPECIFIED 

RP=    REGRES*2.*SHAPE/S1ZE*FLREM**{  1 .-1 ./SHAPE )*WTFUEL**( 1. /SHAPE) 

GO    TO    50 
45    CONTINUE 
C  FUEL    SURFACE    CONTROL 

C  ASSUME    CRIB    STICK    DENSITY     RHOCR=    500    KG/M**3 

RHOCR=    500. 

REGREN=    1.24E-3/RHOCR*SI2E**-0.6 

RPl=     RE6REN*2.*SHAPE/SIZE*FLREM*«( 1 .- I ./SHAPE)* WT FUEL** ( 1 ./SHAPE) 
C  CRIB    POROSITY    CONTROL 

RP2=    0.22*WTFUEL/{RHOCR*SIZE)*SH 

RP=    AMINl     (RP1.RP2) 
50     RMF=    RMA+    RP 

YC02=    3.66e67*CFLPC*RC/100./RMF 

YH20=     (WFLPC*RP+9.0*HFLPC*RC)/100./RMF 

Y02=     (0.23*RMA-R0*RC)/RMF 

YN2=     0.77*RMA/RMF    +NFLPC*RP/l 00 ./RMF 

YPYR=     (RP-RC)/RMF 

IF(YPYR.LT..O)     YPYR=    0. 

MW0UT=    44.*YC02+18.*YH20+28.*YN2+32.*Y02+MWPYR*YPYR 

HRATIO=     l./( l.  +  ( (TGAS/TAMB)*( MWIN/MWOUT  ) ♦ ( 1 . +RP/RMA )**2 ) 
1  ♦♦0.3333333333) 

C  NOTE    HIN    IS    TAKEN    AS    POSITIVE 

HIN=  HWDOW*  HRATIO 

ZW=l .-MWOUT^TAMB/MWIN/TGAS 

IF(ZW)195t55.55 
55    VAVGIN=     0.ee6667^SQRT(2.^G+HlN^ZW) 

RMA=    CD^VAVGIN^HlN^BWDOW»DENSA 

IF  (RP-RMA/R)  60.60t65 
60  RC=  RP«BPF 

BWORST=  BWDOW^RC^R/RMA 

RMA=  RC*R/BPF 

GO    TO    70 
65     RC=     BPF*RMA/R 

BWORST=    BWDOW 
C  RECALCULATE    Y-     VALUES     SINCE    RP.RC    HAVE    BEEN    CHANGED 

70     KR=    KR+1 

IF     (KR-3)     50.75*75 
75    CONTINUE 

QFLOW=    RMF^(YCO2+(TGAS^(0.5«CPCO2( I ) ♦TGAS+CPC02( 2) )-TAMB*(0 .5* 

1  CPC02( 1 )«TAMB+CPC02(2) ) )    +YH20+(T GAS*( 0 . 5*CPH20( 1 ) *TGAS+ 

2  CPH20(2) )-TAMB*{0.5*CPH2O( 1)*TAMB+CPH20(2) ))     +Y02*(TGAS*( 

3  0.5*CP02( 1 )*TGAS  +  CP02(2) )-T AMB* ( 0 .5*CP02 (  I  )*TAMB+CP02 ( 2)  ) ) 

4  +YN2*(TGAS*<0.5*CPN2( 1 ) *TGAS+CPN2 ( 2 ) )-T AMB*( 0. 5*CPN2 ( I )* 

5  TAMB+CPN2(2)) )     +YPYR* (T GAS *( 0. S^CPPYR ( I ) *TGAS+CPPYR{ 2) ) 

6  -TAMB^(0.5*CPPYR( I ) ♦TA MB+CPPYR ( 2) ) )) 
QFIRE=     RC^CVNET 

IF     (ADIA)     GOTO    90 

CALL     DESOLV 

QRADW=     AWALLN*EMS( 1 )*SI GMA^( TGAS^^4.-TSF^*4 . ) 

QCONW=    AWALLN*(TGAS-TSF)*CNV*( { TGAS-TSF )*( T GAS-TSF ) )**0 . 1 66  66667 
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90    CONTINUE 

QRAOO=    HWOOW*BWORST*SIGMA*(TGAS4i4>4.-TAMB*4>4  •  ) 

K=     K-l-l 

F3=F2 

F2  =  F1 

Fl=    QFIRE-QFLOW-QRADO-QRADW-QCONW 

TGAS3=TGAS2 

TGAS2=TGAS1 

TGAS1=TGAS 

IF     (Fl.LT.O. .AND.TGAS.lt .TGASP)    TGASP=TGAS 

IF     (Fl  .GT«0..AND.TGAS.6T.TGASN>    TGASN=TGAS 

0ERIV2=  DERIVI 

IF  (TGAS1.EQ.TGAS2)  GOTO  130 

DERI V1=(F1-F2)/(T GAS X-TGAS2} 

IF     (KTRACE.GT.O)     WRITE    (4t99)     TGASl tTGAS2 t Fl ,F2 .DERI V 1 ,K. KD , 
1  KH,JtT2(l) tTSFfQFIRE,QFLOW,QRADW,RP,RC 

99    FORMAT (2F9.2t3( IPE9. 2) . 31 3 1 15 t 2( 0PF9. 2>  t3(  IPE10.3> .2C0PF7.3) ) 

IF     (.NOT. SCAN)     GOTO    95 

IF     {F1/F2.GE.0.0)     GOTO    93 

SCAN=    .FALSE. 

GOTO     100 
93    TGAS=    TGAS-DTGAS 

IF     ( (PLFUEL.AND.(TGAS.LT.TBOIL)) .OR.(TGAS.LT.TAMB) )     GOTO    200 

GOTO  120 
95  IF  (DERIV2. LT..O. AND. ABS(F2) .GT.. 0001 >GOTO  100 

IF(OERIV2.LT..O.ANO.J.GT.2)  GOTO  100 

TGAS=  TGASl  +  OTGAS 

GOTO  120 
100  DIF=  ABSCFl/QFLOW) 

IF  (DIF.LT.0.002.AND.ABS(TGAS2-TGAS1).LT.2. )  GOTO  130 

TGAS=(F1*TGAS2-F2*TGAS1)/(F1-F2) 

IF     (K. GT. 10. AND. Fl.LT.O.. AND. TGAS.GT. TGASP)     GOTO    105 

IF     (K. GT. 10. AND. Fl .GT.O.. AND.TGAS.lt .TGASN)     GOTO     105 

IF     (TGAS.GT. TGASP. OR.TGAS.lt. TGASN)    TGAS=(T GASP+TGASN) /2. 

IF     (K.EQ.l.AND.KH.EQ.O)    TGAS=    TGASl+10. 

IF     (TGAS.GT. 2000.)     GOTO     110 

IF     (TGAS.LT.(TAMB+30.) )     GOTO     110 

IF     (PLFUEL.AND.TGAS.LT.TBOIL)     GOTO    110 

GOTO    120 
105    TGAS=     (TGASN+TGASP)/2. 

GOTO    120 
110    SCANS     .TRUE. 

TGAS=    1900. 
120    CONTINUE 

IF     (K-200)     30.30.200 
130    CONTINUE 

CALL    RSTA 

FLREM=  FLREM-RP*OTlME 

IF(FLREM.LT.O)  FLREM=0. 

IF  (QCONW.GT.O. )  QWLSUM=  QWLSUM+( QRADW+QCON W )*DTI ME 

IF  (TTIME  .GE.  MTIME)  GO  TO  210 

IF  (TGAS.LE.353..AND.J.GE.10)  GO  TO  210 

IF  ( J.EQ.l )  GO  TO  150 

IF  < JP.LT.JPRINT)  GO  TO  160 

JP=  0 
150  CALL  OUTPUT 
160  JP=  JP+1 

TTIME=  TTIME+DTIME 
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170    CONTINUE 

END    TIME    STEP    OO-LOOP 
PET URN 

ERROR    IN     INPUT 
190    CONTINUE 
KNTRL=    2 
WRITE     (2.910) 
910    FORMAT     (///•     PFLFIX    ROUTINE    DOES    NOT     ACCEPT     STEADY-STATE    CASE') 
195    CONTINUE 

1F(KTRACE.EQ. 1)     WRITE(2«930)     TGAS .RC.RP. YPYR . ZW .RMA, MWOUT 
930    FORMAT     (/•     TGAS=' . F5 .0. •     RC=«.E10.4.«     RP=«.E10.4. 

I         •     YPYR=',E10.*. •     ZW=«.F6.4.«     RMA= • , El  0.4.  •     MWOUT  =  «  t F6 . 1  ) 
FAIL    TO    CONVERGE.     ERROR     EXIT 
200    CONTINUE 
KNTRL=3 
RETURN 

FIRE    IS    OVER    (TRANSIENT    CASE) 
210    CONTINUE 

CALL    OUTPUT 
220    CONTINUE 
RETURN 
END 
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SUBROUTINE    POOL 
C 

C  PLASTIC    FUEL     (POOL    FIRE)     ROUTINE 

C  FUEL    'SEES*     ONLY    COMPARTMENT    AND    NOT     ITSELF 

C 

COMMON    /CNSTS/       AWALLN , BWDOW , DENSA , G . GASCNT , KTRACE , MTI ME 

COMMON    /CP/  CPA,CPC0(2J,CPC02(2) .CPH2<2) »CPH20(2) .CPN2(2) t 

1  CP02(2> tCPPYR(2J 

COMMON    /FUEL/  C .CFLPC» CVGROS* CVNET tH ,HFLPC . MWPYR .N, NFLPC. O. 

1  OFLPC.RtR0tREGRES»SHtSHAPEtSlZE.WtWFLPC»  WTFUEL 

COMMON    /GP/  AWDOW»BPF,CDtCNVt DTlMEtEMS( 2) .HWDOW.IX. IXC. IXL» 

1  J« JMtJPtJPRI NT tK«KO«KH,K ITER (KNTRL«MWIN«MW0UT*RC*RP. SIGMA 

COMMON    /LOGIC/       FC.FLSPECf KRI T« NE WPLT. NEWPRP »PLFUEL. PLOT . PNCH . 
I  RPSPEC.VTSPEC 

COMMON    /PLAST/    TBOILC. DHP, STOI CH, SI ZEl .EIT A , EISCAN 

COMMON    /PRBLM/       ADI A, AFLOOR. AWALL t OENSW, FLO AD,IRUN .OPENF , 
1  PRNT.STEADYiTHICKW 

COMMON    /QS/  QCONtt,QFlRE»QFLOW.QRADO»QRADW.QWLSUM 

COMMON    /TEMP/  DENF. DENU.T AM B.TGAS. TINPT .T 1 ( 20 ) . T2( 20 ) tTSF ,T SU 

COMMON    /THERML/    CNDA( 2 1 10 ) . CPW( 2. I  0 ) .DX . EF, EMSA( 2 . I  0) » 
1  NCN0»NCPWtNEMS«NQGENtNRPtQGEN(2t 10) «RPX(  2*50) 

COMMON    /WOUT/  BWORST»FLREMt HRATI O t RMA.RMF ,TTI ME . VAVGI N, 

1  WAtWBt YC02t YH20.YN2t Y02tYPYR 

LOGICAL    AOIA, EI SCANt FCt FLSPEC » KRI T. NEWPRP. PLFUEL . 
1  PLOT .PNCH.RPSPECf SCAN, STEADY, STOICH,VT SPEC 

REAL    MWIN, MWOUT, MWPYR, MTI ME, N, NFLPC 
C  Y-    DENOTE    MASS    FRACTIONS    OF    OUTFLOW 

C  RMA=     MASS     INFLOW    RATE 

C  RMF=    MASS    OUTFLOW    RATE 

C  R=    STOICHIOMETRIC    AIR/FUEL    MASS    RATIO 

C  R0=    STOICHIOMETRIC    OXYGEN/FUEL    RATIO 

SCAN=     .FALSE. 
C  INITIALIZES    STOICHIOMETRIC    CASE 

IF     (STOICH)     EITA    =1. 

QRAOW=0. 

QCONW=    0. 

F2=0. 

F1=0. 

DTGAS  =  10. 

TBOIL=TBOlLC+273. 

FC=  .FALSE. 

CALL  HEAONG 

TF2=  TGAS 
C  START  TIME  STEP 

DO  170  J=1.JM 

KH=  0 

DER1V1=  I. 

TGAS2=  0. 

TGAS1=  0- 

TGASP=  2000. 

TGASN=  TAMB 
20  CONTINUE 

K=    0 
30    CONTINUE 

IF     (STOICH)     GOTO    34 
32    RP=    SIZE*EF«SIGMA*(TGAS**4.-TB0IL**A.)/DHP 

IF     (STEADY)     GOTO    33 


I 


58 


PLUME=    SIZE*0,0014*CVNET/DHP 

PROP=     1 .-(TGAS**A.-TB0IL**4.)/( I700.**4.-TBOIL«*4 . ) 

IF     (PROP.LT.O.)     PROP=    0. 

RP=    RP+PROP*PLUME 

33  CONTINUE 
RMF=     RMA+RP 
1F(FLREM.LE.0)RP=0. 

34  YC02=    3.66667*CFLPC«RC/100./RMF 
YH20=     (WFLPC*RP+9.0*HFLPC*RC)/100./RMF 
Y02=     C0.23*RMA=R0«RC)/RMF 

YN2=     0.77*RMA/RMF    +NFLPC*RP/1 00 ./RMF 
YPYR=     (RP-RC)/RMF 
IF(YPYR.LT. .0)     YPYR=    0. 
QFUEL=     (RP-RC)*DHP 

MWOUT=    44.«YC02+18.*YH20+28.*YN2+32«*Y02+MWPYR*YPYR 
HRAT2  0=    l./( l.+( (TGAS/TAMB)*(MWlN/MWOUT)*( 1 .+RP/RMA) ♦♦2) 
1  **0 .3333333333) 

NOTE    HlN    IS    TAKEN    AS    POSITIVE 
HIN=     HWDOW*     HRATIO 
ZW= 1 .-MWOUT*T  AM 8/ MW I N/T  GAS 
IF<ZW)195t35t35 

35  VAVGIN=    0.ee6667*SQRT(2.«G«HIN*ZW) 
RMA=    CO*VAVGlN*HIN*BWDOW*DENSA 
RMF=    RMA+RP 

RC=  BPF*RMA/R 

IF  (STOICH)  RP=  RC/BPF 

IF    < (EISCAN. AND. ElTA.LT.l.). OR. STOICH)     GOTO    37 

36  IF     (RC.GT.RP*BPF)     GOTO    40 

37  FC=     .FALSE. 
GO    TO    45 

40    FC=     .TRUE. 

RC=    RP*BPF 

RECALCULATE    VALUES     IF     IN    FUEL    CONTROL    REGIME. 

YC02=    3.66667*CFLPC*RC/100./RMF 

YH20=     (WFLPC*RP+9.0*HFLPC*RC)/100./RMF 

Y02=     {0.23*RMA=R0*RC)/RMF 

YN2=     0.77*RMA/RMF    +NFLPC*RP/l 00 ./RMF 

YPYR=  (RP-RC)/RMF 
45  CONTINUE 

QFLOW=  RMF»(YC02*(TGAS*(0.5*CPC02(1 )*TGAS+CPC02( 2 ) )-TAMB*{0.5* 

1  CPC02( 1 )*TAMB+CPC02<2) ) )     +YH20*( TGAS* ( 0 . 5«CPH20 ( 1 ) *T6AS+ 

2  CPH20(2) )-TAMB*( 0.5*CPH2O( l)*TAMB+CPH20(2)))     +Y02*(TGAS*( 

3  0.5*CPO2{l )*TGAS  +  CP02<2) )-T AMB*( 0.5*CP02 (  I  )«TAMB+CP02( 2)  ) ) 

4  +YN2*(TGAS*(0.5*CPN2( I ) ♦TGAS+CPN2 ( 2) )-TAMB*( 0.5*CPN2( 1 )* 

5  TAMB+CPN2(2) ) )     +YP YR* (T GAS* ( 0 . 5*CPPYR ( 1 ) ♦TGAS+CPPYR ( 2) ) 

6  -TAMB«(0 .5*CPPYR( I )«TAMB+CPPYR ( 2) ) ) > 
QFIRE=    RC*CVNET 

IF     (AOIA)     GOTO    90 
IF     ( .NOT. STEADY)     CALL    DESOLV 
IF     (STEADY)     CALL     STFLOW 

QRADW=    AWALLN*EMS(  1  )*SI  GMA*(T  GAS**4.-=TSF**4  .  ) 

QCONW=     AWALLN*(TGAS-TSF)*CNV*( ( TGAS-TSF) *(TGAS-TSF ) )**0 . 16666667 
90     CONTINUE 

QRADO=    AWOOW*SI GMA*(TGAS**4.-TAMB**4 . ) 

K=     K+1 

F3  =  F2 

F2=F1 

Fl=     QFIRE-QFLOW-QFUEL-QRAOO-QRADW-QCONW 
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TGAS3=TGAS2 

TGAS2=TGAS1 

TGAS1=TGAS 

IF     tFl.LT.0..ANO.TGAS.LT.TGASP)    TGA£P=TGAS 

IF     (Fl.GT.O. .AND.TGAS.GT.TGASN)    TGASN=TGAS 

DERIV2=  DERIVl 

IF  (TGAS1.EQ.TGAS2)  GOTO  130 

DER1VI=(F1-F2)/(TGASI-TGAS2) 

IF     (KTRACE.GT.O)     WRITE    (4,99)     T GAS! .TGAS2. F 1 ,F2.DERI V  I , K .KD , 
I  KH.JtT2(  1  J  tTSF,QFIRE,QFLOW,QRADW.RP,RC 

99    F0RMAT(2F9.2.3( 1PE9.2),3I3.15.2( 0PF9.2).3( lPEiO.3) .2(0PF7.3) ) 

IF     (.NOT. SCAN)     GOTO    95 

IF     (F1/F2.GE.0. 0)     GOTO    93 

SCAN=     .FALSE. 

GOTO     100 
93    TGAS=    TGAS-DTGAS 

IF     (TGAS.LT.TAMB)     GOTO    200 

IF     (TGAS.LT.TBOIL.ANO.(FLREM.GT.O.O) )     GOTO     190 

GOTO     120 
95     IF     (DERIVl .LT..0. AND. ABS(F2) .GT..0001)GOTO     100 

IF(DERIV2.LT..0.AND.J.GT.2)     GOTO    100 

TGAS=    TGASl+OTGAS 

GOTO     120 
100    DIF=    ABS(Fl/QFLOW) 

IF     (DIF.LT.0.002.ANO.ABS(TGAS2=TGASl ).LT.2. )     GOTO     130 

TGA£=(FI*TGAS2-F2*TGAS1)/(F1-F2) 

IF     (K.GT.10.AND.F1.LT.0..AND.TGAS.GT.TGASP)     GOTO     105 

IF     (K. GT. 10. AND. Fl.GT.O. .ANO.TGAS.lt .TGASN)     GOTO     105 

IF     (K.EQ.l .AND.KH.EQ.O)     TGAS=    TGASl+10. 

IF     (TGAS.GT.2000.  )     GOTO     110 

IF     (TGAS.LT.(TAMB+30.) )     GOTO    110 

IF    (TGAS.LT.TBOIL.AND.(FLREM.GT.O.O))     GOTO     110 

GOTO     120 
105    TGAS=     (TGASN+TGASP)/2. 

GOTO  120 
110  SCAN=  .TRUE. 

TGAS=  1900. 
120  CONTINUE 

IF  (STEADY. AND. .NOT. ADIA)  CALL  STFLOW 

IF  (K-200)  30.30.200 
130  IF  (STEADY)  GOTO  180 

CALL  RSTA 

FLREM=  FLREM-RP*DTIME 

IF(FLREM.LT.O)  FLREM=0. 

IF  (QCONW.GT.O. )  QWLSUM=  QWLSUM+( QRADW+QCONW )*DTI ME 

IF  (TTIME  .GE.  MTIME)  GO  TO  210 

IF  (TGAS.LE.353. .AND.J.GE. 10)  GO  TO  210 

IF  (J.EQ.l )  GO  TO  150 

IF  (JP.LT.JPRINT)  GO  TO  160 

JP=  0 
150  CALL  OUTPUT 
160  JP=  JP+1 

TTIME=  TTIME+DTIME 
170  CONTINUE 

END  TIME  STEP  DO-LOOP 
180  CONTINUE 

IF(  .NOT.STOICH)  GOTO  185 

FIND    STOICHIOMETRIC    FUEL    SIZE 
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SIZE1=    RP/(EF*SIGMA*(TGAS**4.-TB0IL**4.)/DHP) 
185    CALL    OUTPUT 

NORMAL    EXIT     WHEN    STEADY. EQ.T 
RETURN 

ERROR    EXIT 
190    CONTINUE 

IF     (KTRACE.EQ.l)     WRITE    (2t910>    TGAS  . 
910    FORMAT!///'     TGAS.LT.TBOIL  TGAS=« t F8. 1 . •  GOTO    NEXT    CASE'///) 

GOTO    200 

SQUARE    ROOT    ERROR 
195    CONTINUE 

IF(KTRACE.EQ.l )     WRITE(2t930)     TGAS t RC.RP. YPYR t ZW  ,RMA,  MWOUT 
930    FORMAT     (/'     TGAS=' t F5 .0 . '     RC=',E10.4.'     RP='.E10.4» 

1         '     YPYR=',E10.4.  '     ZW='.F6.4t'     RMA='  .El  0.4 t  •     MWOUT  =  '  . F6. I ) 
FAIL    TO    CONVERGE.     ERROR    EXIT 
200    CONTINUE 
KNTRL=3 
RETURN 

FIRE    IS     OVER     (TRANSIENT    CASE) 
210    CONTINUE 

CALL    OUTPUT 

RETURN 

END 
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SUBROUTINE  PP 

PLOTTING  SUBROUTINE 

THIS    ROUTINE    IS    LEFT     BLANK    SINCE     IT     IS    MACHINE-DEPENDENT 

ENTRY    PLTRST 

THIS    ENTRY    SETS    UP   THE     INIALIZATION    OF    PLOTTING 
ENTRY    DSTO 

THIS    ENTRY    IS    CALLED    EACH    TIME    TO    STORE    A    DATA    POINT 
ENTRY    DOUT 

THIS    IS    THE    LAST    ENTRY    FOR    A    GIVEN    RUN 
RETURN 
END 


I 


62 


SUBROUTINE  PVTFIX 
C 

C  PESSIMIZATION  ROUTINE 

C  FIXED  VENTILATIONt WORST  POSSIBLE  FUEL  PYROLYSIS  RATE. 

C 

COMMON  /CNSTS/   AWALLN, BW DOW, DENS A, G . GASCNT , KTRACE , MTI ME 

COMMON  /CP/       CPA,CPC0(2J.CPC02(2) .CPH2( 2) .CPH20(2>,CPN2(2) . 
1    CP02(2) ,CPPYR{2) 

COMMON    /FUEL/  C, CFLPCt CVGROS, CVNET t H,HFLPC. MWPYR , N. NFLPC,  O. 

1  OFLPC,R.R0tREGRES.SH,SHAPE,SIZE,W,WFLPC» WTFUEL 

COMMON    /GP/  AWDOW,BPF,CD»CNV. DTIMEtEMS(2>,HWD0W»IXtIXCtIXL. 

1  J.JMtJPt  JPRI NT t K. KD » KH » KIT ER,KNTRL,MWlN.MWOUT»RCtRP, SIGMA 

COMMON    /LOGIC/       FC, FLSPEC, KRITt NEWPLT, NEWPRP .PLFUEL.PLOT .PNCH , 
1  RPSPECt VTSPEC 

COMMON    /PLAST/    T BOlLC, OHP , STOlCH, SI ZEl , EITA , EI  SCAN 

COMMON    /PRBLM/        ADI A , AFLOOR , A  WALL , DENS W»FLOAD, I  RUN, OPENF , 
1  PRNT.STEADYtTHICKW 

COMMON    /QS/  QCONW,QFIRE,QFLOW, QRADO»QRADW,QWLSUM 

COMMON    /TEMP/  DENF, DENUtT AMB.TGAS, TINPT ,T 1 ( 20i .T2( 20 ) ,TSF ,T SU 

COMMON    /THERML/    CNDA( 2, I  0 ) , CPW( 2, 1 0)  , DX , EF, EMSA{ 2 , 1 0) , 
1  NCND,NCPW,NEMS,NQGEN,NRP,QGEN(2, 10) ,RPX(2,50) 

COMMON    /WOUT/  BWORSTtFLREM,HRATI0,RMAtRMF,TTlME,VAVGlN, 

1  WA,WB,YC02,YH20, YN2, Y02,YPYR 

LOGICAL    ADIA,EISCAN,FC,FLSPEC,KRIT,NEWPRP,PLFUEL, 
1  PLOT,PNCH,RPSPEC,SCAN,STEADY,STOICH,VTSPEC 

REAL    MWIN,MW0UT»MWPYR,MTIME,N,NFLPC 

IF     (STEADY)     GOTO     190 

FC=     .FALSE, 

SCAN=    .FALSE. 

QRADW=0. 

QCONW=    O. 

F2=0. 

F1=0. 

DTGAS=10. 

CALL    HEADNG 
C  START    TIME    LOOP 

DO    170     J=1,JM 

KH=     0 

DERIV1=     1. 

TGAS2=    0. 

TGAS1=    0. 

TGASP=    2000. 

T6ASN=    TAMB 
20    CONTINUE 

K=    0 
30    CONTINUE 

IF     (FLREM.GT.O.  )     GOTO    32 

RC=    0. 

RP=    0. 

FC=     .TRUE. 
32    RMF=    RMA+RP 

YC02=    3.6e667*CFLPC«RC/100./RMF 

YH20=     (WFLPC*RP+9.0«HFLPC«RC>/100./RMF 

Y02=     (0.23«RMA-R0«RC)/RMF 

YN2=  0.77«RMA/RMF  +NFLPC*RP/1 00 ./RMF 

YPYR=  (RP-RCJ/RMF 

1F{YPYR.LT..0)  YPYR=  0. 
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MWOUT=    44.*YC02+18.*YH20+28.*YN2+32.*Y02+MWPYR*YPYR 

HRATIO=     l./( I .+( (TGAS/TAMB)*(MWlN/MWOUT)*( I . ♦RP/RMA)**2 ) 
1  **0. 3333333333) 

NOTE    HIN    IS    TAKEN    AS    POSITIVE 

HIN=    HWDOW*     HRATIO 

ZW=1  .-MW0UT*TAMB/MWIN/TGAS 

IF( ZW)195t35.35 
35    VAVGIN=    0.666667*SQRT(2.*G*HIN*ZW) 

RMA=    CD*VAVGlN*HlN«BWOOW*DENSA 

RMF=     RMA+RP 

IF     (.NOT.FC)     RC=     BPF*RMA/R 

IF     (.NOT.FC)     RP=    RC/BPF 
45    CONTINUE 

QFLOW=    RMF«(YC02*(TGAS«(0.5«CPC02< 1 )*TGAS+CPC02 ( 2 ) )-TAMB*(0 .5* 

1  CPC02(1 J*TAMB+CPC02(2} ) )    +YH20*( T GAS*( 0 . 5*CPH20 ( 1 ) *TGAS+ 

2  CPH20(2i )-TAMB*(0«5*CPH2O( l)*TAMB+CPH20(2J ))     +Y02*(TGAS*( 

3  0.5*CPO2(l )*T6AS+CP02(2) >-T AMB*( 0 .5*CP02 ( 1 )*TAMB+CP02( 2) ) ) 

4  +YN2«(TGAS«(0.5*CPN2( I } *TGAS+CPN2( 2) )-T AMB*( 0 . 5«CPN2 ( 1 ) * 

5  TAM8+CPN2(2}> )     +YPYR* (T GAS* ( 0 . 5*CPPYR ( 1 > ♦TGAS+CPPYR( 2} ) 

6  -TAMB*{0.5*CPPYR( I > ♦TAMB+CPPYRt 2) ) )) 
QFIRE=    RC*CVNET 

IF  (ADIA)  GOTO  90 

CALL  DESOLV 

QRADW=  AWA'_LN*EMS(  1  )*SI  GMA*(T  6AS**4  .-TSF**4  .  ) 

QCONW=  AWALLN*(TGAS™TSF>*CNV*( (T6AS-TSF ) ♦(T GAS-TSF > }**0 . I  6666667 
90  CONTINUE 

QRAOO=  AWDOW*SI GMA* (TGAS**4 .-TAMB**4 . ) 

K=  K+1 

F3  =  F2 

F2=F1 

Fl=  QFIRE-QFLOW-QRADO-QRADW-QCONW 

TGAS3=TGAS2 

TGAS2=TGASl 

TGAS1=TGAS 

IF     (F1.LT.0..AND.TGAS.LT.TGASPJ     TGASP=TGAS 

IF     (Fl.GT.O. .AND.TGAS.GT.TGASN)    TGASN=TGAS 

DERIV2=  OERIVl 

IF     {TGAS1.EQ.TGAS2J     GOTO     130  i 

DERIV1=(FI-F2)/(TGAS1-TGAS2) 

IF  (KTRACE.GT.O)  WRITE  (4f99)  TGASl .TGAS2. F 1 .F2 .DERI  V 1 , K , KD , 
t     KHt J.T2( I )»TSF,QFIREtQFLOW.QRADW.RP.RC 
99  FORMAT (2F9. 2 t3{ 1PE9. 2) ,313 . 15 . 2( 0PF9 .2 ) .3< IPEIO.3) ,2(0PF7,3) } 

IF  (.NOT. SCAN)  GOTO  95 

IF  (F1/F2.GE.0.0}  GOTO  93 

SCAN=  .FALSE. 

GOTO     100 
93    TGAS=    TGAS-OTGAS 

IF     (TGAS.LT.TAMB)     GOTO    200 

GOTO  120 
95  IF  (OERIVl. LT..O. AND. ABS(F2)  .GT-.OOODGOTO  100 

IF(DERIV2.LT..0.AND.J.GT.2}  GOTO  100 

TGAS=    TGASl+OTGAS 

GOTO     120 
100    DIF=    ABS(F1/QFL0W> 

IF     {DIF.LT.0.002.AND.ABS(TGAS2-TGAS1) .LT.2.  )     GOTO     130 

TGAS=(F1*TGAS2-F2*TGAS1 )/(Fl-F2> 

IF     (K.GT.lO.AND.Fl .LT.O..AND.TGAS.GT.TGASP)     GOTO    105 

IF     (K.GT.lO.AND.Fl .GT.O-. ANO.TGAS.lt .TGASN)     GOTO     105 
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IF     (K.EQ.l.ANO.KH.EQ.O)    T GAS=    TGASl     +10. 

IF     (TGAS.GT.2000. )     GOTO     110 

IF     (TGAS.LT.(TAMB+30.J )     GOTO    110 

GOTO     120 
105    TGAS=     (TGASN+TGASP)/2. 

GOTO     120 
110    SCAN=     .TRUE. 

TGAS=     1900. 
120    CONTINUE 

IF     (K-200)     30t30»200 
130    CONTINUE 

CALL    RSTA 

FLREM=    FLREM-RP^DTIME 

IF(FLREM.LT.O)     FLREM=0. 

IF     (QCONW.GT.O. )     QWLSUM=     QWLSUM+( QRADW+QCCNW )*DTI ME 

IF     (TTIME     .GE.     MTIMEJ     GO    TO    210 

IF     (TGAS.LE.353..AND.J.GE.10)     GO    TO    210 

IF     (J.EQ.l )     GO    TO     150 

IF     { JP.LT. JPfiINT)     GO   TO    160 

JP=    0 
150    CALL    OUTPUT 
160    JP=    JP+1 

TTIME=    TTIME+DTIME 
170    CONTINUE 

END    TIME    STEP    DO-LOOP 

RETURN 

ERROR  IN  INPUT 
190  CONTINUE 

KNTRL=  2 

WRITE     (2.910) 
910    FORMAT     (///•     PVTFIX    ROUTINE    DOES    NOT    ACCEPT     STEADY-STATE    CASE'} 

RETURN 

SQUARE  ROOT  ERROR 
195  CONTINUE 

IF(KTRACE.EQ.l)  WRITE(2.930)  TGAS. RC. RP , YPYR , ZW ,RMA , MWOUT 
930  FORMAT  (/•  TG AS  =  « .F5. 0»  •  RC  =  «.E10.4.'  RP=«.E10.4, 

1    •  YPYR=« ,E10.4. •  ZW='.F6.4.*  RMA= • , E 10.4, •  MWOUT=« ,F 6. 1 ) 
FAIL  TO  CONVERGE,  ERROR  EXIT 
200  CONTINUE 

KNTRL  =  3 

RETURN 

FIRE  IS  OVER  (TRANSIENT  CASE) 
210  CONTINUE 

CALL  OUTPUT 

RETURN 

END 
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SUBROUTINE  RPFI X 


TABULAR  FUEL  PYROLYSIS  ROUTINE 

FUEL  PYROLYSIS  RATE  IS  AN  INPUT  VARIABLE. 

COMMON  /CNSTS/   A WALLN, BWOOW, DENSA .G « GASCNT » KTRACE .MTI ME 

COMMON  /CP/       CPA,CPC0(2).CPC02(2) .CPH2(2) .CPH20(2)  tCPN2(2)  • 
1     CP02{2).CPPYR(2) 

COMMON    /FUEL/  C. CFLPCt CVGROSt CVNET . H, HFLPC .MWPYR ,N , NFLPC . O, 

I  OFLPCt R.RO.REGRES.SH. SHAPE t SI ZEtW.WFLPC.WTFUEL 

COMMON    /GP/  AWOOWtBPF.CDtCNV. DTIME,EMS<2)»HWD0W,IX,IXC,IXL, 

I  J, JMt JP» JPRINT«KtKD«KH.KlTER.KNTRL»MWlN. MWOUT . RC. RP t SI GMA 

COMMON    /LOGIC/       FC» FLSPECt KRIT t NEWPLT. NEWPRP ,PLFUEL t PLOT .PNCH . 
1  RPSPEC. VTSPEC 

COMMON    /PLAST/    T BOlLCt DHP, STOICH, SI ZEl , EITA , EI  SCAN 

COMMON    /PRBLM/       ADI At AFLOOR t AWALL t DENSW.FLO AD .1  RUN, OPENFt 
1  PRNT.STEADY.THICKW 

COMMON    /QS/  QCONW,QFIREiQFLOW.QRADO»QRAOW,QWLSUM 

COMMON    /TEMP/  DENF, OENU.T AMB.TGAS.T INPT . T 1 ( 20 ) .T2( 20 ) ,TSF ,TSU 

COMMON    /THERML/    CNDA( 2,  1 0 ) , CPW( 2, I  0) .DX ,EF, EMSA( 2 . 1 0 ) , 
1  NCND.NCPW.NEMS.NQGENf NRPt  QGEN(2. 10)«RPX( 2.50} 

COMMON    /WOUT/  BWORST»FLREM.HRATIO.RMA.RMF,TTlME. VAVGiNt 

I  WAt WB»YC02tYH20.YN2» Y02,YPYR 

LOGICAL    ADIA.EISCAN.FC.FLSPEC.KRITtNEWPRP.PLFUEL* 
I  PLOT. PNCH, RPSPEC. SCAN, STEADY, STOICH, VTSPEC 

REAL    MWI N.MWOUT, MWPYR. MTI ME. N. NFLPC 

SCAN=     .FALSE. 

QRAOW=0. 

QCONW=    0. 

F2  =  0. 

Fl  =  0. 

DTGAS=10. 

CALL    HEADNG 

START    TIME    LOOP 

DO     170    J=1.JM 

KH=     0 

DERIVl=     1. 

TGAS2=    0.  ' 

TGASl=    0. 

TGASP=    2000. 

TGASN=    TAMB 
20    CONTINUE 

K=    0 
30     CONTINUE 

FC=     .FALSE. 

IF     (FLREM.GT.O. )     RP=    TLU ( RPX, NRP. TTI ME> 

IF     (FLREM.LE.O.i     RP=    0. 

RMF=    RMA+RP 

YC02=    3.66667*CFLPC«RC/100./RMF 

YH20=     (WFLPC*RP+9.0*HFLPC*RC)/100 ./RMF 

Y02=     (0.23*RMA-'R0*RC>/RMF 

YN2=     0.77*RMA/RMF    +NFLPC*RP/100 ./RMF 

YPYR=     (RP-RC)/RMF 

IF(YPYR.LT. .0)     YPYR=    0. 

MWOUT=    44.*YC02+18.*YH20+28.*YN2+32.*Y02+MWPYR*YPYR 

HRATIO=     l./( l.+( (TGAS/TAMB)*(MWIN/MW0UT)*(1.+PP/RMA)**2) 
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1     «*0 .3333333333) 

NOTE    HIN    IS    TAKEN    AS    POSITIVE 

HIN=  HWDOW*  HRATIO 

ZW=l.-MWOUT*TAMB/MWlN/TGAS 

IF{ZW>195.35.35 
35  VAVGIN=  0.6e6e67*SQRT(2.*G*HIN*ZW) 

RMA=    CD*VAVGIN*HIN*BW00W«DENSA 

RMF=    KMA+fiP 

IF     (RMA/R-RP)     40.*0t45 
40    RC=    BPF*RMA/'R 

GO    TO    50 
45    RC=    BPF*RP 

FC=     .TRUE. 
50    CONTINUE 

QFLOW=    RMF*(YC02*(TGAS*(0.5*CPC02( 1 ) *TGAS+CPCD2( 2 ) ) -TAME* (0.5* 

1  CPC02( 1 )*TAMB+CPC02(2)) )    +YH20*(T GAS*( 0 . 5«CPH20( 1 ) ♦TGAS+ 

2  CPH20(2) )-TAMB*{0.5*CPH2O{ l)*TAMB+CPH20(2i) )     +Y02*(TGAS*< 

3  0.5*CPO2(l )*TGAS+CP02(2} )-T AMB*{ 0  .5*CP02 (  1  )*TAMB+CP02 ( 2) ) ) 

4  +YN2*(TGAS*( 0.5*CPN2(1 }*TGAS+CPN2 ( 2  J )"TAMB*( 0 . 5*CPN2( 1 )* 

5  TAMB+CPN2C2)) i    +YPYR*( TGAS*( 0. 5*CPPYR( 1 ) ♦TGAS+CPPYRC 2 ) ) 

6  ■=TAMB«(0.5*CPPYR(  1  )  *T  AMB+CPPYR  (2)  )  )) 
QFIRE=    RC*CVNET 

IF     (AOIA)     GOTO    90 

IF  (  .NOT. STEADY)  CALL  OESOLV 

IF  (STEADY)  CALL  STFLOW 

QRADW=  AWALLN*EMS( 1 )*SI GMA*(TGAS**A.-TSF**4 . ) 

QCONW=  AWALLN*(TGAS-TSF)*CNV*( (TGAS-TSF) ♦( T GAS-TSF ) )**0. 1 6666667 
90  CONTINUE 

QRADO=  AWOOW*SI GMA*( TGAS**4 .-T AMB**4 . ) 

K=  K+l 

F3=F2 

F2  =  Fl 

F 1=    OFIRE-QFLOW-QRADO-QRADW-QCONW 

TGAS3=TGAS2 

TGAS2=TGAS1 

TGAS1=TGAS 

IF     (Fl.LT.O. .ANO.TGAS.lt .TGASP)    TGASP=TGAS 

IF     (Fl,GT.O..AND.TGAS.GT.TGASN)     T6ASN=TGAS 

DERIV2=  DERIVl 

IF  (TGASl .EQ.TGAS2)  GOTO  130 

DERIV1=(F1-=F2)/(TGAS1-TGAS2) 

IF    (KTRACE. GT.O)     WRITE    (4*99)     TGASl .TGAS2. F I . F2 . DERI V  I . K. KO . 
I  KH, J,T2( 1) .TSF.QFlRE.QFLOW.QRADW.RPtRC 

99    FORMAT ( 2F9.2.3( IPE9. 2) t 313. 15. 2( 0PF9.2 ) . 3 ( 1 PEl 0 .3 ) ,2(0PF7.3) ) 

IF     (.NOT. SCAN)     GOTO    95 

IF     (F1/F2.GE.0. 0)      GOTO    93 

£CAN=     .FALSE. 

GOTO     100 
93    TGAS=    TGAS-DTGAS 

IF     (TGAS.LT.TAMB)     GOTO    200 

GOTO     120 
95    IF     (DERIVl .LT..O. AND. ABS(F2) .GT..0001)GOTO     100 

IF(DERIV2.LT. .O.AND. J.GT.2)     GOTO     100 

TGAS=    TGASl+OTGAS 

GOTO     120 
100     DIF=     ABS(F1/QFL0W) 

IF     (DIF.LT.0.002.AND.ABS(TGAS2-TGAS1) .LT.2.  )     GOTO     130 

TGAS=(F1*TGAS2-F2*TGASI)/(F1-F2) 
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IF  (K.GT.lO.AND.Fl .LT.O..AND.TGAS.GT.TGASP)  GOTO  105 

IF  (K.GT.lO.AND.Fl .GT.O.. AND.TGAS.lt. TGASN)  GOTO  105 

IF  (K.EQ.l. AND.KH.EQ.O)  TGAS=  TGASl+lO. 

IF  (TGAS.GT.2000.>  GOTO  I  10 

IF  (TGAS.LT.(TAMB+30.) )  GOTO  110 

GOTO     120 
105    TGAS=     (TGASN+TGASP)/2. 

GOTO     120 
110    SCAN=     .TRUE. 

TGAS=     1900. 
120    CONTINUE 

IF     (STEADY. AND.. NOT. ADIA)     CALL     STFLOVK 

IF     (K-200)    30.30.200 
130    CONTINUE 

IF     (STEADY)     GOTO     180 

CALL    RSTA 

FLREM=    FLREM-RP^DTIME 

IF(FLREM.LT.O. J     FLREM=0. 

IF     (QCONW.GT.O.)     QWLSUMs:    QWLSUM-l-(  QRADW-I-QCON  W  )*DTI  ME 

IF     (TTIME     .GE.     MTlMEi     GO    TO    210 

IF     (TGAS.LE.353..ANO.J.GE.10)     GO    TO    210 

IF     (J.EQ.l )     GO    TO     150 

IF     ( JP.LT. JPRINT)     GO   TO    160 

JP=    0 
150    CALL    OUTPUT 
160    JP=    JP+l 

TTIME=    TTIME+DTIME 
170    CONTINUE 

END    TIME    STEP    00-LOOP 
180    CONTINUE 
185    CALL    OUTPUT 

NORMAL    EXIT     WHEN    STEADY. EQ.T 

RETURN  , 

SQUARE    ROOT    ERROR 
195    CONTINUE 

IF(KTRACE.EQ.l)     WRITE(2i930)     TGAS.RC .RP. YPYR» ZW ,RMAt MWOUT 
930    FORMAT     (/•     TGAS=« iF5 .0. •     RC=«.E10.4.»     RP=«,E10.4. 

1        •     YPYR=».E10.4. •     ZW=«.F6.4t«     RMA= • t El  0.4. •     MWOUT  =  « . F6 . 1 ) 
FAIL    TO    CONVERGE.     ERROR    EXIT 
200    CONTINUE 

KNTRL=3 

RETURN 

FIRE    IS    OVER    (TRANSIENT    CASE) 
210    CONTINUE 

CALL  OUTPUT 

RETURN 

END 
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II 


SUBROUTINE    STFLOW 
C 

C  CALCULATES     WALL    HEAT     CONDUCTION    WHEN    STEADY-STATE 

C  CONDITION    ONLY    IS    NEEDED, 

C 

COMMON    /GP/  AWDOWtBPF.CD,CNV,DTlME«EMS(2) .HWOOWtlX.IXCt IXLt 

1  J. JM. JP. JPRINTtK.KOtKH.KITERfKNTRLtMWlN.MWOUT.RC.RP.SIGMA 

COMMON    /PRBLM/       AOIAt AFLOORt AWALL, DENS W.FLO AD t I RUN.OPENF» 
1  PRNT.STEADY.THICKW 

COMMON    /TEMP/  DENF »DENU.TAMB,TGAS.TI NPT .T 1 ( 20) .T 2( 20 } .TSF t TSU 

COMMON    /THERML/     CNDAC2. 1 0 > . CPW( 2. 1 0) . DX . EF , EMSA ( 2 . 10 ) . 
1  NCNDtNCPWtNEMS*NQGEN«NRPtQGEN( 2« 10) «RPX( 2 t 50 ) 

C  BIOT=    BIOT     NUMBER 

KD=    0 

TSF=    TGAS    -30. 

TSU=    TAMB    +30. 
10    CONTINUE 

TSFOLD=    TSF 

TSUOLD=    TSU 

EMS(1)=     l./( l./TLU{EMSA.NEMS.TSF)     +l./EF-l.i 

EMS(2)=    TLU(EMSA. NEWS. TSU) 

TAVG=     (TGAS+TAMB)/2. 

CND=    TLU(CNDA.NCNO»TAVG) 

ZRF=    TGAS*(TGAS*(TGAS-i-TSF)+TSF*TSF)+TSF*TSF*TSF 

ZCF=    CNV*( (TGAS-TSF)*(TGAS-TSF) J**0. 16666667 

HF=    ZCF+EMS( 1 )*S1GMA«ZRF 

B10TF=    HF*THICKW/CND 

ZRU=    TAMB«(TAMB*{TAMB+TSU)+TSU«TSU)+TSU*TSU*TSU 

ZCU=     1.31*(  (TAMB-TSU)*(TAMB-TSU))** 0.166666  67 

HU=     ZCU+EMS(2)*SIGMA*ZRU 

BI0TU=    HU*THICKW/CND 

TSF=     ( (BlOTF+HF/HU)*TGAS+TAMB)/( l.+BIOTF+HF/HU) 

TSU=     ( (BI0TU+HU/HF)*TAMB+TGAS)/(1 . +BIOTU+HU/HF) 

T2(l)=    TSF-     {TSF-TSU)*DX/THICKW/2. 

T2(IXC)=     (TSF+TSU>/2. 

IF     ( (ABS<TSF-TSFOLD) .LT.3.) .AND. ( ABS{TSU-TSUOLD) . LT .3 . ) ) 
1  RETURN 

KD=    KD+1 

TSU=     <TSU+TSU0LD)/2. 

IF     <KD.LT.20)     GOTO     10 

WRITE     (2.90)    TSF.TSFOLD 
90    FORMAT     (//•     ♦♦UNSUCCESSFUL     ITERATION    IN    STFLOW'/ 
1  •  TSF=«,F15.2. •       TSFOLD=« ,F15.2) 

RETURN 

END 
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FUNCTION  TLU  ( AR«AY,NUM, VALIN ) 

TABULAR  LOOK-UP  INTERPOLATING  ROUTINE 

DIMENSION  ARRAY(2«NUM) 

ARRAY(1.I>=    INDEPENDENT     VARIABLE 
ARRAY(2»I>=    DEPENDENT     VARIABLE 

INTERPOLATES    LINEARLY    WITHIN    GIVEN    DOMAIN.       SETS    EQUAL    TO 
SMALLEST    OR    LARGEST    VALUE     IF    OUTSIDE    THE    DOMAIN. 
IF     (NUM.NE.l)     GO    TO     10 
TLU=    ARRAY (2. 1 ) 
RETURN 
10     IF     (NUM.NE.2)     GO    TO    20 
1=     2 

GO    TO    50 
20    IF     CVALIN.GT.ARRAY(i.l) )     GO    TO    30 
TLU=    ARRAY (2. 1 ) 
RETURN 
30    DO    AO     1=2. NUM 

IF     <VALIN.LE.ARRAY(1 , I))     GO    TO    50 
40    CONTINUE 

TLU=    ARRAY(2.NUM) 
RETURN 
50    TLU=    ARRAY(2.I-1)     +    (VALIN    -     ARRAY(  1  , I- I ) )* 

1  ((ARRAY(2tI}     -     ARRAY{2»I-1))     /     (ARRAYd.I)    -    ARRAY  (  1  , 1- 1  )  )  ) 

RETURN 
END 


SUBROUTINE    TRIDGF     (  A,  B«  C*  O.  E,  IX  ) 

TRIDIAGONAL    GAUSS     ELIMINATION    PROCEDURE  FOR    UNSYMMETRIC 

MATRictS. 

A=LEFT     OF     DIAGONAL.     B=DI  AGONAL.     C=RIGHT  OF    DIAGONAL. 

D=    CONSTANT    VECTOR.     E=    SOLUTION    VECTOR.  I X=    SIZE    OF    MATRIX. 

DIMENSION    A(20).B(20).C(20) .D(20) .E(20) .CP( 20) 

CP(l)=    C(l )/B(l ) 

Ed  )=    DdJ/BC  I) 

C(IX)=    0. 

IXL=     IX- 1 

00     10    I=29lX 

J=    I-l 

BX=    B(IJ-CP( J>*A(I) 

CP(I)=    C(I)/BX 

E(I)=     (D(I>-E(J}«ACIJ)/BX 
10    CONTINUE 

DO    20    1=1. IXL 

J=     IX-I 

E(J)=    E( J)-E( J+1)*CP(J) 
20    CONTINUE 

RETURN 

END 
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UFACTURERS advances  in  the  field  of  voluntary  and 
mandatory  standards.  The  new  DIMENSIONS/NBS  also 
carries  complete  listings  of  upcoming  conferences  to  be 
held  at  NBS  and  reports  on  all  the  latest  NBS  publications, 
with  information  on  how  to  order.  Finally,  each  issue  carries 
a  page  of  News  Briefs,  aimed  at  keeping  scientist  and  consum- 
alike  up  to  date  on  major  developments  at  the  Nation's  physi- 
cal sciences  and  measurement  laboratory. 
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Enter  my  Subscription  To  DIMENSIONS/NBS  at  S11.00.  Add  S2.75  for  foreign  mailing.  No  additional 
postage  is  required  for  mailing  within  the  United  States  or  its  possessions.  Domestic  remittances 
should  be  made  either  by  postal  money  order,  express  money  order,  or  check.  Foreign  remittances 
should  be  made  either  by  intefnational  money  order,  draft  on  an  American  bank,  or  by  UNESCO 
coupons. 
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NBS  TECHNICAL  PUBLICATIONS 


PERIODICALS 

JOURNAL  OF  RESEARCH— The  Journal  of  Research 
of  the  National  Bureau  of  Standards  reports  NBS  research 
and  development  in  those  disciplines  of  the  physical  and 
engineering  sciences  in  which  the  Bureau  is  active.  These 
include  physics,  chemistry,  engineering,  mathematics,  and 
computer  sciences.  Papers  cover  a  broad  range  of  subjects, 

j  with  major  emphasis  on  measurement  methodology,  and 
the  basic  technology  underlying  standardization.  Also  in- 
cluded from  time  to  time  are  survey  articles  on  topics  closely 
related  to  the  Bureau's  technical  and  scientific  programs.  As 
a  special  service  to  subscribers  each  issue  contains  complete 
citations  to  all  recent  NBS  publications  in  NBS  and  non- 
NBS  media.  Issued  six  times  a  year.  Annual  subscription: 
domestic  $17.00;  foreign  $21.25.  Single  copy,  $3.00  domestic; 
$3.75  foreign. 
Note:  The  Journal  was  formerly  published  in  two  sections: 

'  Section  A  "Physics  and  Chemistry"  and  Section  B  "Mathe- 
matical Sciences." 
DIMENSIONS/NBS 

This  monthly  magazine  is  published  to  inform  scientists, 
engineers,  businessmen,  industry,  teachers,  students,  and 
consumers  of  the  latest  advances  in  science  and  technology, 
with  primary  emphasis  on  the  work  at  NBS.  The  magazine 
highlights  and  reviews  such  issues  as  energy  research,  fire 
protection,  building  technology,  metric  conversion,  pollution 
abatement,  health  and  safety,  and  consumer  product  per- 
formance. In  addition,  it  reports  the  results  of  Bureau  pro- 
grams in  measurement  standards  and  techniques,  properties 
of  matter  and  materials,  engineering  standards  and  services, 
instrumentation,  and  automatic  data  processing. 

Annual  subscription:    Domestic,  $11.00;  Foreign  $13.75 

NONPERIODICALS 

Monographs — Major  contributions  to  the  technical  liter- 
ature on  various  subjects  related  to  the  Bureau's  scientific 
and  technical  activities. 

Handbooks — Recommended  codes  of  engineering  and  indus- 
trial practice  (including  safety  codes)  developed  in  coopera- 
tion with  interested  industries,  professional  organizations, 
and  regulatory  bodies. 

Special  Publications — Include  proceedings  of  conferences 
sponsored  by  NBS,  NBS  annual  reports,  and  other  special 
publications  appropriate  to  this  grouping  such  as  wall  charts, 
pocket  cards,  and  bibliographies. 

Applied  Mathematics  Series — Mathematical  tables,  man- 
uals, and  studies  of  special  interest  to  physicists,  engineers, 
chemists,  biologists,  mathematicians,  computer  programmers, 
and  others  engaged  in  scientific  and  technical  work. 
National  Standard  Reference  Data  Series — Provides  quanti- 
tative data  on  the  physical  and  chemical  properties  of 
materials,  compiled  from  the  world's  literature  and  critically 
evaluated.  Developed  under  a  world-wide  program  co- 
ordinated by  NBS.  Program  under  authority  of  National 
Standard  Data  Act  (Public  Law  90-396). 


NOTE:  At  present  the  principal  publication  outlet  for  these 
data  is  the  Journal  of  Physical  and  Chemical  Reference 
Data  (JPCRD)  published  quarterly  for  NBS  by  the  Ameri- 
can Chemical  Society  (ACS)  and  the  American  Institute  of 
Physics  (AIP).  Subscriptions,  reprints,  and  supplements 
available  from  ACS,  1155  Sixteenth  St.  N.W.,  Wash.,  D.C. 
20056. 

Building  Science  Series — Disseminates  technical  information 
developed  at  the  Bureau  on  building  materials,  components, 
systems,  and  whole  structures.  The  series  presents  research 
results,  test  methods,  and  performance  criteria  related  to  the 
structural  and  environmental  functions  and  the  durability 
and  safety  characteristics  of  building  elements  and  systems. 
Technical  Notes — Studies  or  reports  which  are  complete  in 
themselves  but  restrictive  in  their  treatment  of  a  subject. 
Analogous  to  monographs  but  not  so  comprehensive  in 
scope  or  definitive  in  treatment  of  the  subject  area.  Often 
serve  as  a  vehicle  for  final  reports  of  work  performed  at 
NBS  under  the  sponsorship  of  other  government  agencies. 
Voluntary  Product  Standards — Developed  under  procedures 
published  by  the  Department  of  Commerce  in  Part  10, 
Title  15,  of  the  Code  of  Federal  Regulations.  The  purpose 
of  the  standards  is  to  establish  nationally  recognized  require- 
ments for  products,  and  to  provide  all  concerned  interests 
with  a  basis  for  common  understanding  of  the  characteristics 
of  the  products.  NBS  administers  this  program  as  a  supple- 
ment to  the  activities  of  the  private  sector  standardizing 
organizations. 

Consumer  Information  Series — Practical  information,  based 
on  NBS  research  and  experience,  covering  areas  of  interest 
to  the  consumer.  Easily  understandable  language  and 
illustrations  provide  useful  background  knowledge  for  shop- 
ping in  today's  technological  marketplace. 
Order  above  NBS  publications  from:  Superintendent  of 
Documents,  Government  Printing  Office,  Washington,  D.C. 
20402. 

Order  following  NBS  publications — NBSIR's  and  FIPS  from 
the  National  Technical  Information  Services,  Springfield, 
Va.  22161. 

Federal  Information  Processing  Standards  Publications 
(FIPS  PUB) — Publications  in  this  series  collectively  consti- 
tute the  Federal  Information  Processing  Standards  Register. 
Register  serves  as  the  official  source  of  information  in  the 
Federal  Government  regarding  standards  issued  by  NBS 
pursuant  to  the  Federal  Property  and  Administrative  Serv- 
ices Act  of  1949  as  amended.  Public  Law  89-306  (79  Stat. 
1127),  and  as  implemented  by  Executive  Order  11717 
(38  FR  12315,  dated  May  11,  1973)  and  Part  6  of  Title  15 
CFR  (Code  of  Federal  Regulations). 

NBS  Interagency  Reports  (NBSIR) — ^A  special  series  of 
interim  or  final  reports  on  work  performed  by  NBS  for 
outside  sponsors  (both  government  and  non-government). 
In  general,  initial  distribution  is  handled  by  the  sponsor; 
public  distribution  is  by  the  National  Technical  Information 
Services  (Springfield,  Va.  22161)  in  paper  copy  or  microfiche 
form. 


BIBLIOGRAPHIC  SUBSCRIPTION  SERVICES 


The  following  current-awareness  and  literature-survey  bibli- 
ographies are  issued  periodically  by  the  Bureau: 
Cryogenic  Data  Center  Current  Awareness  Service.  A  litera- 
ture survey  issued  biweekly.  Annual  subscription:  Domes- 
tic, $25.00;  Foreign,  $30.00. 
Liquified  Natural  Gas.  A  literature  survey  issued  quarterly. 
Annual  subscription:  $20.00. 


Superconducting  Devices  and  Materials.  A  literature  survey 
issued  quarterly.  Annual  subscription:  $30.00.  Send  subscrip- 
tion orders  and  remittances  for  the  preceding  bibliographic 
services  to  National  Bureau  of  Standards,  Cryogenic  Data 
Center  (275.02)  Boulder,  Colorado  80302. 


U.S.  DEPARTMENT  OF  COMMERCE 
National  Bureau  of  Standards 

Washington,  D.C.  20234 


OFFICIAL  BUSINESS 

Penalty  for  Private  Use,  S300 


POSTAGE  AND  FEES  PAID 

U.S.  DEPARTMENT  OF  COMMERCE 

COM-215 


^ 


SPECIAL  FOURTH-CLASS  RATE 
BOOK 


